mallab程序示例.docxVIP

第三章MATLAB编程与仿真重点：1、M文件M 脚本文件和M函数文件2、程序结构2．1顺序结构2．2循环结构for循环while循环3、选择结构3．1、switch…case结构3．2、if…else…end结构4、函数4．1、子函数4．2、私有函数4．3、嵌套函数……例3-1a=input(a=?);b=input(b=?);c=input(c=?);d=b*b-4*a*c;x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp([x1=,num2str(x(1)),,x2=,num2str(x(2))]);例3-4%输出0-50之间能被3整除的数disp(能整除3的数);for n =0:50if mod(n,3)~=0continueenddisp(n);end%输出0-50之间3和5的公倍数disp(输出0-50之间3和5的公倍数);for n =0:50if mod(n,3)==0if mod(n,5)~=0continueenddisp(n);endend选择结构例3-7clear alldisp(switch 语句)c=input(请输入成绩;);switch c case num2cell(90:100),disp(A);case num2cell(80:89),disp(B);case num2cell(70:79),disp(C);case num2cell(60:69),disp(D);otherwisedisp(E);end例3-8cleardisp(if…..else语句!)x=input(请输入分数：);if(x=100x=90)disp(A);elseif(x=89x=80)disp(B);elseif(x=79x=70)disp(C);elseif(x=69x=60)disp(D);elseif(x60)disp(E);end文件名：EX38.m 文件类型:函数文件例3-11clear allload mandrillimage(X)colormap(map)axis offaxis image文件名：picture.m 文件类型:M脚本文件例%脚本程序%花瓣图theta=-pi:0.01:pi;rho(1,:)=2*sin(5*theta).^2;rho(2,:)=cos(10*theta).^3;rho(3,:)=sin(theta).^2;rho(4,:)=5*cos(3.5*theta).^3;for k=1:4 polar(theta,rho(k,:))pauseend文件名:petals.m 文件类型:M脚本文件例3-12ax(1)=subplot(1,2,1);rgb=imread(ngc6543a.jpg);image(rgb);title(真彩色图像)ax(2)=subplot(122);im=mean(rgb,3);image(im); title(图像增强)colormap(hot(256))linkaxes(ax,xy)axis(ax,image)运行失败例3-14h2=@magicfunctions(h2) %显示句柄函数feval(h2,5) %调用句柄函数例3-17function[avg,med]=newstats(u) %主函数%在主函数中寻找子函数mean和mediann=length(u);avg=mean(u,n);med=madian(u,n);function a=mean(v,n) %子函数%计算平均水平a=sum(v)/n;function m =median(v,n) %子函数%计算中位数w=sort(v);if rem(n,2)==1 m=w((n+1)/2);elseif m=(w(n/2)+w(n/2+1))/2;end分段函数演示x=linspace(0,4);y=pw2(x); %自变量取样点集plot(x,y);%计算相应函数值hold on;%绘制曲线plot(1*ones(1,2),ylim,k:);%绘制曲线间隔线plot(2*ones(1,2),ylim,k:);%绘制曲线间隔线set(gca,fontsize,14);xlable(\itx,fontname,time new roman,fontsize,16);function y=pw2(x)if sum(x=0x4)0.5; %检查定义域是否在规定范围内 error(x must be in the interval [0,4]);endy=x.^2; %计算函数值y(x=1x2)=cos([x(x=1x2)+1]*pi);%计算函数值y(x=2)=-x(x=2).^2./[x(x=2)-2];%计算函数值保存为pw2.m文