matlab答案第六章

1、第六章 1. 利用MATLAB提供的randn函數(shù)生成符合正態(tài)分布的105隨機(jī)矩陣A，進(jìn)行如下操作：(1) A各列元素的均值和標(biāo)準(zhǔn)方差。(2) A的最大元素和最小元素。(3) 求A每行元素的和以及全部元素之和。(4) 分別對A的每列元素按升序、每行元素按降序排列。答：clear all; close all; clc;A=randn(10, 5);meanA=mean(A); %(1)A各列元素的均值stdA=std(A); %(1)A各列元素的標(biāo)準(zhǔn)方差maxA=max(max(A); %(2)A的最大元素minA=min(min(A); %(2)A的最小元素rowsumA=sum(A, 2)

2、; %(3)A每行元素的和sumA=sum(rowsumA); %(3)A全部元素之和sort1=sort(A); %(4)A的每列元素按升序排列sort2=sort(A, 2, descend); %(4)A的每行元素按降序排列2. 按要求對指定函數(shù)進(jìn)行插值和擬合。(1) 按表6.1用3次樣條方法插值計算090D范圍內(nèi)整數(shù)點(diǎn)的正弦值和075D范圍內(nèi)整數(shù)點(diǎn)的正切值，然后用5次多項(xiàng)式擬合方法計算相同的函數(shù)值，并將兩種計算結(jié)果進(jìn)行比較。表6.1 特殊角的正弦和正切值表度0 15 30 45 60 75 90sin0 0.2588 0.5000 0.7071 0.8660 0.96591.0000

3、tan0 0.2679 0.5774 1.0000 1.7320 3.7320(2) 按表6.2用3次多項(xiàng)式方法插值計算1100之間整數(shù)的平方根。表6.2 1100內(nèi)特殊值的平方根表N 1 4 9 16 25 36 49 64 81 100N的平方根 1 2 3 4 5 6 7 8 9 10 答:(1) 程序設(shè)計：clear all; close all; clc;alpha1=0:15:90;sin_alpha1=sin(alpha1*pi/180); %精確正弦值plot(alpha1, sin_alpha1, k:p); hold on; %繪精確正弦曲線alpha2=0:90;sin_

4、Y1=interp1(alpha1, sin_alpha1, alpha2, spline); %3次樣條正弦插值plot(alpha2, sin_Y1, r-*); hold on; %繪3次樣條插值正弦曲線P1=polyfit(alpha1, sin_alpha1, 5); %5次多項(xiàng)式擬合sin_Y2= polyval(P1, alpha2); %5次多項(xiàng)式求值plot(alpha2, sin_Y2, b-o); %繪5次多項(xiàng)式插值正弦曲線legend(精確正弦值, 3次樣條正弦插值, 5次多項(xiàng)式正弦插值); title(正弦值比較); alpha3=0:15:75;tan_alpha3

5、=tan(alpha3*pi/180); %精確正切值figure, plot(alpha3, tan_alpha3, k:p); hold on; %繪精確正切曲線alpha4=0:75;tan_Y1=interp1(alpha3, tan_alpha3, alpha4, spline); %3次樣條正切插值plot(alpha4, tan_Y1,r-*); hold on; %繪3次樣條正切曲線P2=polyfit(alpha3, tan_alpha3, 5); %5次多項(xiàng)式擬合tan_Y2= polyval(P2, alpha4); %5次多項(xiàng)式求值plot(alpha4, tan_Y2

6、, b-o); %繪5次多項(xiàng)式插值正弦曲線legend(精確正切值, 3次樣條正切插值, 5次多項(xiàng)式正切插值); title(正切值比較); (2)程序設(shè)計：clear all; close all; clc;X=1, 4, 9, 16, 25, 36, 49, 64, 81, 100; Y=1:10;X1=1:100; Y1=interp1(X, Y, X1, cubic);plot(X, Y, r:o); hold on; %繪精確平方根曲線plot(X1, Y1, k-x); %繪3次多項(xiàng)式插值平方根曲線legend(精確平方根, 3次多項(xiàng)式插值);3. 已知一組實(shí)驗(yàn)數(shù)據(jù)如表6.3所示。

7、表6.3 一組實(shí)驗(yàn)數(shù)據(jù)i12345Xi165 123 150 123 141 Yi187 126 172 125 148 求它的線性擬合曲線。 答：程序設(shè)計：clear all; close all; clc;x=165, 123, 150, 123, 141;y=187, 126, 172, 125, 148;P=polyfit(x, y, 3)運(yùn)行結(jié)果：P =1.0e+003 *-0.0000 0.0013 -0.1779 8.4330所以它的線性擬合曲線為：P(x)=1.3x2-177.9x+84334. 已知多項(xiàng)式P1(x)=3x+2，P2(x)=5x2x+2，P3(x)=x20.5，

8、求 ： (1) P(x)=P1(x)P2(x)P3(x)。(2) P(x)=0的全部根。(3) 計算xi=0.2i (i=0,1,2,10)各點(diǎn)上的P(xi)。答：(1) 程序設(shè)計：clear all; close all; clc;P1=0, 3, 2; P2=5, -1, 2; P3=1, 0, -0.5;P=conv(conv(P1, P2), P3)運(yùn)行結(jié)果：P =0 15.0000 7.0000 -3.5000 0.5000 -2.0000 -2.0000所以P(x)=15x5+7x4-3.5x3+0.5x2-2x-2(2) 程序設(shè)計：r=roots(P)運(yùn)行結(jié)果：r =0.7071

9、0.1000 + 0.6245i0.1000 - 0.6245i-0.7071-0.6667(3) 程序設(shè)計：i=0:10;xi=0.2*i;Pxi=polyval(P, xi)運(yùn)行結(jié)果：Pxi =-2.0000 -2.3920 -2.6112 -1.7024 2.7104 15.0000 42.1120 94.1408 184.9056 332.5264 560.00005. 求函數(shù)在指定點(diǎn)的數(shù)值導(dǎo)數(shù)。答：(1) 程序設(shè)計：clear all; close all; clc;x=pi/6, pi/4, pi/3;f=inline(sin(x).2+cos(x).2);dx=diff(f(x,

10、 5*pi/12)/(pi/12) %可參見教材第157頁例6.19運(yùn)行結(jié)果：dx =0 0 0%當(dāng)x=pi/2時單獨(dú)計算x=pi/2;f=inline(sin(x).2+cos(x).2);dx=diff(f(x, pi)/(pi/2)dx =(2) 程序設(shè)計：clear all; close all; clc;x=1:3;f=inline(sqrt(x.2+1);dx=diff(f(x, 4)運(yùn)行結(jié)果：dx =0.8219 0.9262 0.96086. 求數(shù)值積分。答：(1) 程序設(shè)計：clear all; close all; clc;f=inline(sin(x).5.*sin(5*

11、x);qf=quad(f, 0, pi)運(yùn)行結(jié)果：qf =0.0982(2) 程序設(shè)計：clear all; close all; clc;f=inline(1+x.2)./(1+x.4);qf=quad(f, -1, 1)運(yùn)行結(jié)果：qf =2.2214(3) 程序設(shè)計：clear all; close all; clc;f=inline(x.*sin(x)./(1+cos(x).2);qf=quad(f, 0, pi)運(yùn)行結(jié)果：qf =2.46747. 分別用矩陣求逆、矩陣除法以及矩陣分解求線性方程組的解。答：(1) 程序設(shè)計：方法1：矩陣求逆法clear all; close all; c

12、lc;A=2, 3, 5; 3, 7, 4; 1, -7, 1;b=10; 3; 5;x=inv(A)*b運(yùn)行結(jié)果：x =-1.8060-0.53733.0448方法2：左除運(yùn)算符法clear all; close all; clc;A=2, 3, 5; 3, 7, 4; 1, -7, 1;b=10; 3; 5;x=Ab運(yùn)行結(jié)果：x =-1.8060-0.53733.0448方法3：矩陣LU分解法clear all; close all; clc;A=2, 3, 5; 3, 7, 4; 1, -7, 1;b=10; 3; 5;L, U=lu(A);x=U(Lb)運(yùn)行結(jié)果：x =-1.8060-

13、0.53733.0448 (2) 程序設(shè)計：方法1：矩陣求逆法clear all; close all; clc;A=5, 1, -1, 0; 1, 0, 3, -1;-1,-1,0, 5; 0, 0, 2, 4;b=1; 2; 3; -1;x=inv(A)*b運(yùn)行結(jié)果：x =1.4000-5.90000.1000-0.3000方法2：左除運(yùn)算符法clear all; close all; clc;A=5, 1, -1, 0; 1, 0, 3, -1;-1,-1,0, 5; 0, 0, 2, 4;b=1; 2; 3; -1;x=Ab運(yùn)行結(jié)果：x =1.4000-5.90000.1000-0.3

14、000方法3：矩陣LU分解法clear all; close all; clc;A=5, 1, -1, 0; 1, 0, 3, -1;-1,-1,0, 5; 0, 0, 2, 4;b=1; 2; 3; -1;L, U=lu(A);x=U(Lb)運(yùn)行結(jié)果：x =1.4000-5.90000.1000-0.30008. 求下列方程的根。答：(1) 程序設(shè)計：clear all; close all; clc;f=inline(x-sin(x)./x);x=fzero(f, 0.5)運(yùn)行結(jié)果：x =0.8767(2) 程序設(shè)計：clear all; close all; clc;f=inline(sin(x).2).*exp(-0.1.*x)-0.5.*abs(x);x=fzero(f, 1.5)運(yùn)行結(jié)果：x =1.67389. 求非線性方程組在(0.5,0.5)附近的數(shù)值解。答：函數(shù)文件fxy.m：function f=fxy(u)x=u(1); y=u(2);f(1)=x-0.6*sin(x)-0.3*cos(y);f(2)=y-0.6*cos(x)+0.3*sin(y);在命令窗口輸入以下命令：clear all; close all; clc;x=f