资源描述:
《matlab典型例题.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、例2.1>>muw0=1.785;>>a=0.03368;>>b=0.000221;>>t=0:20:80;>>muw=muw0./(1+a*t+b*t.^2)例2.2数值数组和字符串的转换>>a=[1:5];>>b=num2str(a);>>a*2ans=246810>>b*2ans=986464100646410264641046464106例2.9比较左除和右除求解恰定方程>>rand('seed',12);>>a=rand(100)+1.e8;>>x=ones(100,1);>>b=a*x;>>cond(a)ans=5.0482e+011>>
2、tic;x1=b'/a;t1=toct1=0.4711>>er1=norm(x-x1')er1=139.8326>>re1=norm(a*x1'-b)/norm(b)re1=4.3095e-009>>tic;x1=ab;t1=toct1=0.0231>>tic;x1=ab;t1=toct1=0.0011>>er2=norm(x-x1)er2=1.5893e-004>>re1=norm(a*x1-b)/norm(b)re1=4.5257e-016例2.14:计算矩阵的指数>>b=magic(3);>>expm(b)ans=1.0e+006*1.08
3、981.08961.08971.08961.08971.08971.08961.08971.0897例2.18:特征值条件数>>a=[-149-50-154;537180546;-27-9-25]a=-149-50-154537180546-27-9-25>>[V,D,s]=condeig(a)V=0.3162-0.4041-0.1391-0.94870.90910.9740-0.00000.1010-0.1789D=1.00000002.00000003.0000例2.415阶多项式在【0,2pi】最小二乘拟合>>x=0:pi/20:pi/2;>>
4、y=sin(x);>>a=polyfit(x,y,5);>>x1=0:pi/30:pi*2;>>y1=sin(x1);>>y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^3+a(4)*x1.^2+a(5)*x1+a(6);>>plot(x1,y1,'b-',x1,y2,'r*')>>legend('原曲线','拟合曲线')>>axis([0,7,-1.2,4])例3.7gradient绘制矢量图>>x=0:pi/20:pi/2;>>y=sin(x);>>a=polyfit(x,y,5);>>x1=0:pi/30:pi*2;>>y1
5、=sin(x1);>>y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^3+a(4)*x1.^2+a(5)*x1+a(6);>>plot(x1,y1,'b-',x1,y2,'r*')>>legend('原曲线','拟合曲线')>>axis([0,7,-1.2,4])>>>>[x,y]=meshgrid(-2:.2:2,-2:.2:2);>>z=x.*exp(-x.^2-y.^2);>>[px,py]=gradient(z,.2,.2);>>contour(z),>>holdon>>quiver(px,py)>>holdoff例基本绘
6、图命令rand(100,1);plot(y)例4.1绘制如图>>x=1:0.1*pi:2*pi;>>y=sin(x);>>z=cos(x);>>plot(x,y,'--k',x,z,'-.rd')例4.5绘制如图>>x=1:10;>>y=rand(10,1);>>bar(x,y);>>x=0:0.1*pi:2*pi;>>y=x.*sin(x);>>feather(x,y)例4.6绘制如图>>lim=[0,2*pi,-1,1];>>fplot('[sin(x),cos(x)]',lim)例4.7绘图如下>>x=[2,4,6,8];>>pie(x,{'m
7、ath','english','chinese','music'})例4.9绘图如下三维螺旋线>>x=0:pi/50:10*pi;>>y=sin(x);>>x=0:pi/50:10*pi;>>y=sin(x);>>z=cos(x);>>plot3(x,y,z);例4.10绘图如下。矩阵三维图>>[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);>>z=x.*exp(-x.^2-y.^2);>>plot3(x,y,z)例4.13绘图如下>>[X,Y]=meshgrid([-4:0.5:4]);>>Z=sqrt(X.^2+Y.^2);>
8、>meshc(Z)例4.19绘制柱面图>>x=0:pi/20:pi*3;>>r=5+cos(x);>>[a,