1、第二章3.设矩阵 A 为A=24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21;（1） B=A(2:5,1:2:5)B =65 24 2134 98 218 42 12143 45 21（2）A(7)=A =24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21（3）A+30（4）size(A);ndims(A)（5）题目有误（6）reshape(x,3,4)（7）abs(x)（8）char(x)4. L1 =0 0 0
2、0 1 0 0 0 0L2 =1 1 1 1 1 0 0 0 0L3 =0 0 0 1 1 1 0 0 0L4 =4 5 65.(1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D =65.0000 5.00000 32.000054.0000 3.1400E =1.0e+003
3、 *0.9141 -0.22391.2080 2.71231.1330 -0.2103(2)E=10endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*d+100*e+10*b+c;disp(加密后：,num2str(a)2.a=input(请输入a: );b=input(请输入b: );c=input(
4、请输入c: );x=0.5:1:5.5;x1=(x=0.5elseif b x=linspace(-10,10,200); y=; for x0=xif x00y=y,x0.2+(1+x0).(1/4)+5;elseif x0=0y=y,0;elseif x0 plot(x,y)fplot 函数：fplot(x0).*(x.2+(1+x).(1/4)+5),-10,10)第六章1.A=randn(10,5)（1）mean(A) ;均值 std(A) ;标准方差（2）max(max(A) ;最大元素 min(min(A) ;最小元素（3）B=sum(A,2) ;A 每行元素的和 sum(B) ;
5、A 全部元素之和（4）sort(A) ;A 的每列元素按升序排列sort(A,2,descend) ;A 的每行元素按将序排列2.（1）（2）X=1 4 9 16 25 36 49 64 81 100;Y=1:10;X1=1:100;Y1=interp1(X,Y,X1,cubic)3.x=165 123 150 123 141;y=187 126 172 125 148;P=polyfit(x,y,3)P =1.0e+003 *-0.0000 0.0013 -0.1779 8.4330所以它的线性拟合曲线为：p(x)=1.3x 2177.9x+84334.(1)P1=0 3 2;P2=5 -1
6、 2;P3=1 0 -0.5;P=conv(conv(P1,P2),P3)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)roots(P)ans =0.7071 0.1000 + 0.6245i0.1000 - 0.6245i-0.7071 -0.6667 (3)i=0:10;xi=0.2*i;polyval(P,xi)ans =-2.0000 -2.3920 -2.6112 -1.7024 2.7104 15.0000 42.1120 94.1408 184.9056
7、332.5264 560.00005.(1)建立函数文件：function f=fxy(u)x=u(1);y=u(2);f=3.*x.2+2*x.*y+y.2在命令窗口中输入以下命令：U,fmin=fminsearch(fxy,1,1)结果：U =1.0e-004 *-0.0675 0.1715fmin =1.9920e-010（2）f=inline(-sin(x)-cos(x.2);fmax=fminbnd(f,0,pi)fmax =0.73106.（1）x=pi/6 pi/4 pi/3;f=inline(sin(x).2+cos(x).2);dx=diff(f(x,5*pi/12)/(pi
8、/12) 可参见第 157 页例题 6.19dx =0 0 0x=pi/2 时单独计算：x=pi/2;f=inline(sin(x).2+cos(x).2);diff(f(x,pi)/(pi/2)ans =0（2）x=1:3;f=inline(sqrt(x.2+1);dx=diff(f(x,4)结果：dx =0.8219 0.9262 0.96087.（1）f=inline(sin(x).5.*sin(5*x);quad(f,0,pi)ans =0.0982（2）f=inline(1+x.2)./(1+x.4);quad(f,-1,1)ans =2.2214（3）f=inline(x.*sin
9、(x)./(1+cos(x).2);quad(f,0,pi)ans =2.4674（4）f=inline(abs(cos(x+y);dblquad(f,0,pi,0,pi)ans =6.28328.N=64; %采样点数T=5; %采样时间终点t=linspace(0,T,N); %给出 N 个采样时间 ti（i=1:N）y=exp(-t); %求各采样点样本值 ydt=t(2)-t(1); %采样周期f=1/dt; % 采样频率Y=fft(y); %计算 y 的快速傅里叶变换 YF=Y(1:N/2+1); %F(k)=Y(k)f=f*(0:N/2)/N; %使频率轴 f 从 0 开始plot
10、(f,abs(F) %绘制振幅-频率图9.（1）矩阵求逆法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;x=inv(A)*bx =-1.8060-0.53733.0448矩阵除法法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;x=Abx =-1.8060-0.53733.0448矩阵分解法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;L,U=lu(A);x=U(Lb)x =-1.8060-0.53733.0448（2）方法同（1）10.函数文件：line_solution(A,b)function x,y=line_solution(A,b)m,
11、n=size(A);y=;if norm(b)0if rank(A)=rank(A,b)if rank(A)=ndisp(原方程组有唯一解 x);x=Ab;elsedisp(原方程组有无穷个解，特解为 x，齐次方程组的基础解系为 y);x=Ab;y=null(A,r);enddisp(方程组无解);x=;endelsedisp(原方程组有零解 x);x=zeros(n,1);if rank(A) In line_solution at 11方程组无解x =y =-0.5000 0.50001.0000 00 1.00000 011.（1）f=inline(x-sin(x)./x);x=fzer
12、o(f,0.5)x =0.8767(2)f=inline(sin(x).2).*exp(-0.1.*x)-0.5.*abs(x);x=fzero(f,1.5)x =1.673812.函数文件：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)在命令窗口输入以下命令：x=fsolve(fxy,0.5,0.5,optimset(Display,off)结果：x =0.6354 0.373415.A=-1 2 0;-1 2 -1;-1 2 -1;-1 2 -1;0 2 -1;d=-
13、1;0;1;B=spdiags(A,d,5,5);b=1 0 0 0 0;x=(inv(B)*b)x =0.8333 0.6667 0.5000 0.3333 0.1667实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。(1) 012sin85ze(2) ，其中22l()x210.45ix(3) 0.3. 0.3sin(.)l,3.0,29,.3022aeaz(4) ，其中 t=0:0.5:2.524113ttztt解：M 文件 :z1=2*sin(85*pi/180)/(1+exp(2)x=2 1+2*i;-0.45 5;z2
14、=1/2*log(x+sqrt(1+x2)a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a)./2.*sin(a+0.3)+log(0.3+a)./2)t=0:0.5:2.5;z4=(t=0B2解： M 文件 ：A=12 34 -4;34 7 87;3 65 7;B=1 3 -1;2 0 3;3 -2 7;A+6.*BA-B+eye(3)A*B A.*B A3 A.3 A/BBA A,B A(1,3,:);B2运算结果：A=12 34 -4;34 7 87;3 65 7;B=1 3 -1;2 0 3;3 -2 7;A+6.*BA-B+eye(3)A*B A.*B A3 A.3 A/BBA A,B A(1,3,:);B2ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1
