书籍：Matlab实用教程
工具：Matlab2021a
电脑信息：Intel® Xeon® CPU E5-2603 v3 @ 1.60GHz
系统类型：64位操作系统，基于X64的处理器 windows10 专业版
第5章 Matlab程序设计
5.1 脚本文件和函数文件
5.1.1 M文本编辑器
x=0:0.1:20;zeta=0.3y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));plot(x,y,'r')>>Ex0501 zeta=0.3000
5.1.2 M文件的基本格式
>>Ex0501a>>Ex0501b(0.505)5.1.3 M脚本文件
Ex0501c.m
%Ex0501c x=0:0.1:20;y1=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+acos(0.3));plot(x,y1,'r')hold on y2=1-1/sqrt(1-0.707^2)*exp(-0.707*x).*sin(sqrt(1-0.707^2)*x+acos(0.707));plot(x,y2,'g')y3=1-exp(-x).*(1+x)plot(x,y3,'b')>>Ex0501c5.1.4 M函数文件
function y=Ex0501d(zeta)%UNTITLED9 此处显示有关此函数的摘要%此处显示详细说明 x=0:0.1:20;y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));plot(x,y)end>>f=Ex0501d(0)>>hold on>>f=Ex0501d(0.1)>>hold on>>f=Ex0501d(0.3)
5.2 程序流程控制
5.2.1 for…end循环结构
%Ex0503 sum=0;forn=1:2:100sum=sum+n end>>Ex0503 sum=2500%Ex0504 sum=zeros(6,1);forn=eye(6,6)sum=sum+n end>>Ex0504 sum=100000sum=110000sum=111000sum=111100sum=111110sum=1111115.2.2 while…end循环结构
%Ex0505 sum=0;n=1;whilen<=100sum=sum+n n=n+2;end>>Ex0505 sum=25005.2.3 if…else…end条件转移结构
function y=Ex0506(zeta)%UNTITLED9 此处显示有关此函数的摘要%此处显示详细说明 x=0:0.1:20;if(zeta>0)&&(zeta<1)y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));elseif zeta==1y=1-exp(-x).*(1+x);endplot(x,y)>>Ex0506(0.5)5.2.4 switch…case开关结构
%Ex0507formonth=1:12switchmonthcase{3,4,5}season='spring'case{6,7,8}season='summer'case{9,10,11}season='autumn'otherwise season='winter'end end>>Ex0507 season='winter'season='winter'season='spring'season='spring'season='spring'season='summer'season='summer'season='summer'season='autumn'season='autumn'season='autumn'season='winter'5.2.5 try…catch…end试探结构
%Ex0508 n=4;a=magic(n);m=3;b=eye(3);tryc=a*bcatchc=a(1:m,1:m)*b end lasterr>>Ex0508 c=162351110976ans='错误使用*用于矩阵乘法的维度不正确。请检查并确保第一个矩阵中的列数与第二个矩阵中的行数匹配。要执行按元素相乘，请使用'.*'。'5.2.6 流程控制语句
1、break
%Ex0509 sum=0;n=1;whilen<=100ifsum<1000sum=sum+n n=n+2elsebreakend end>>Ex0509 sum=1024n=652、continue
%Ex0510 sum=2;ss=0;forn=3:100form=2:fix(sqrt(n))ifmod(n,m)==0ss=1;break;elsess=0;end endifss==1continue;end sum=sum+n;end sum>>Ex0510 sum=10603、return
4、pause
5、keyboard
6、input
>>a=input('input a number:')input a number:45a=45>>b=input('input a number:','s')input a number:45b='45'>>input('input a number:')input a number:2+6ans=8
)
)
)
)
