MATLAB上机实验报告086092班魏旭阳20091000567一、实验内容实验六循环结构程序设计（1）编写程序实现求两矩阵A*B,A.*B.A*BA=input('输入矩阵A:');B=input('输入矩阵B:');C=zeros(size(A,1),size(B,2));for ii=1:size(A,1)for jj=1:size(B,2)for k=1:size(A,2)C(ii,jj)=C(ii,jj)+A(ii,k)*B(k,jj);endendendC输入矩阵A:[1 2;3 4]输入矩阵B:[4 3;2 1]C =8 520 13A.*BA=input('请输入矩阵A:');B=input('请输入矩阵B:');if size(A)==size(B)[m,n]=size(A);for i=1:mfor j=1:nC(i,j)=A(i,j)*B(i,j);endenddisp('输出A.*B:')Celsedisp('无效')end请输入矩阵A:[1 2;3 4]请输入矩阵B:[4 3;2 1]输出A.*B:C =4 6 6 4（2）编写程序实现求矩阵的转置。

A=input('请输入矩阵A:'); [m,n]=size(A); B=zeros(n,m) for i=1:m for j=1:nB(i,j)=A(j,i) end end disp(B)请输入矩阵A:[5 6;7 8] B =5 7 6 8实验八 二维曲线的绘制 （1）绘制下列曲线：①33xx y -=x=-10:1:10;y=x-(x.^3)/3;plot(x,y)-10-8-6-4-20246810-400-300-200-1000100200300400②2221xey π=x=-1:0.001:1;y=1/2*pi*exp(1/2*x.^2);plot(x,y)-1-0.8-0.6-0.4-0.200.20.40.60.811.41.61.822.22.42.62.8③64222=+y x ezplot('x^2+2*y^2-64',[-8,8])-8-6-4-202468-8-6-4-202468xyx 2+2 y 2-64 = 0④⎩⎨⎧==tr y t r x cos sint=0:2*pi/100:2*pi;x=sin(t);y=cos(t);plot(x,y)（r 视为常量）-1-0.8-0.6-0.4-0.20.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81（2）通过用plot 和fplot 函数绘制xy 1sin的曲线，并分析其区别。

x=-1:0.01:1;y=sin(1./x);plot(x,y)-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81fplot('sin(1/x)',[-1,1],1e-4)-1-0.8-0.6-0.4-0.20.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81实验九 三维曲线和曲面的绘制 （1）绘制下列三维曲线1) ⎪⎩⎪⎨⎧=≤≤==--t z t t ey t e x t t π20,sin cos 20/20/ t=0:2*pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t),z=t;plot3(x,y,z);title('Line in3-D Space');xlabel('X');ylabel('Y');zlabel('Z'); grid on;-1-0.50.51-1-0.500.5102468XLine in 3-D SpaceYZ2) ⎪⎩⎪⎨⎧=≤≤==3210,tz t t y t xt=0:0.001:1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);title('Line in 3-D Space'); xlabel('X');ylabel('Y');zlabel('Z'); grid on;0.20.40.60.810.5100.20.40.60.81XLine in 3-D SpaceYZ（2）绘制下列曲面图1) z=5,|x|≤5,|y|≤5[x,y]=meshgrid(-5:1:5);z=0*x+0*y+5;mesh(x,y,z)-55-5544.555.562)3,3,15),(22≤≤++-=y x yx y x f[x,y]=meshgrid(-3:1:3);z=-5./(1+x.^2+y.^2);mesh(x,y,z)-4-224-4-224-5-4-3-2-10实验十 MATLAB 数据统计处理利用MATLAB 提供的rand 函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：均值和标准方差；最大元素和最小元素；大于0.5的随机数个数占总数的百分比A=rand(30000,1); mean(A)= 0.4977 std(A)=0.2874 max(A)= 1.0000 min(A)= 9.4306e-006 n=0for i=1:30000 if A(i)>0.5 n=n+1 end end p=n/30000 n =14818 p =0.4939实验十一 多项式四则运算有3个多项式32)(,2)(,542)(2322341++=+=+++=x x x P x x P x x x x P ， 1） 求P(x)= P 1(x)+ P 2(x)P 3(x) p4=conv([1 2],[1 2 3])syms x;p1=x^4+2*x^3+4*x^2+5 syms x;p4=x^3+4*x^2+7*x+6 p=p1+p4p = x^4+3*x^3+8*x^2+11+7*x 2）求P(x)的根 P=[1,3,8,7,11];x=roots(P) x =-1.3840 + 1.8317i -1.3840 - 1.8317i -0.1160 + 1.4400i -0.1160 - 1.4400i3) 当x 取矩阵A 的每一元素时，求P(x)的值，其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=5.255.3275.04.12.11AP = [1 3 8 7 11];A=[-1 1.2 -1.4;0.75 2 3.5;0 5 2.5]; Y=polyval(P,A) Y =1.0e+003 *0.0100 0.0382 0.0125 0.0223 0.0970 0.4122 0.0110 1.2460 0.16444)当以矩阵A 为自变量时，求P(x)的值 P = [1 3 8 7 11];A=[-1 1.2 -1.4;0.75 2 3.5;0 5 2.5]; Z=polyvalm(P,A) Z =1.0e+003 *0.0076 -0.1281 -0.0775 0.1328 1.3900 1.1644 0.1824 1.7364 1.5198实验十二 线性方程组求解及函数求极值（1） 求下列方程组的解1） ⎪⎩⎪⎨⎧=+-=++=++57347310532z y x z y x z y xsyms x y z;[x,y,z]=solve('2*x+3*y+5*z=10','3*x+7*y+4*z=3','x-7*y+z=5','x','y','z') x = -121/67 y = -36/67 z = 204/672） ⎪⎪⎩⎪⎪⎨⎧=+-=-++=-+--=+-+1129312243134945256421432143214321x x x x x x x x x x x x x x xsyms x1 x2 x3 x4;[x1,x2,x3,x4]=solve('6*x1+5*x2-2*x3+5*x4=-4','9*x1-x2+4*x3-x4=13','3*x1+4*x2+2*x3-2*x4=1','3*x1-9*x2+2*x4=11','x1','x2','x3','x4')x1 =2/3 x2 = -1 x3 =3/2 x4 = 0（2）求下列函数在指定区间的最大值1）)2,0(,11)(42∈++=x xx x ffunction fx=mymin(x) fx=-(1+x.^2)/(1+x.^4); x=fmin('mymin',0,2) x =0.6436 f =-1.2071fmax=-f=1.20712）),0(,cos sin )(2π∈+=x x x x f fun=inline('-sin(x)-cos(x.^2)','x'); [x,f]=fminsearch(fun,0,pi)x =0.7310 f =-1.5282fmax=-f=1.5282。