当前位置：首页 > news >正文

Matlab学习记录25

news 2026/5/11 19:06:24

书籍：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')

>>Ex0501c

5.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=111111

5.2.2 while…end循环结构

%Ex0505 sum=0;n=1;whilen<=100sum=sum+n n=n+2;end

>>Ex0505 sum=2500

5.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=65

2、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=1060

3、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