哈工大《结构力学》ch02 m文件与matlab程序设计知识课件.pptVIP

Ch02 M文件与MATLAB程序设计;教学过程;提问2;1. M文本编辑器的基本知识;% 新建M-file, area.m % 计算一个球的体积 r = input(Type radius 输入半径:); Area=pi*r^2; volume=(4/3)*pi*r^3; fprintf(半径 The radius is %12.5f\n,r) fprintf(面积 The area of a circle is %12.5f\n,Area) fprintf(体积 The volume of a sphere is %12.5f\n,volume) 在指令窗口运行指令 area;1. M文本编辑器的基本知识;（1）M文本编辑器;M文件的基本属性*;M文件的组成部分*;M文件的组成部分及其功能* ;（2）脚本式M文件 ;脚本式M文件举例;（3）函数式M文件;例题;function s=jiecheng(n) %此函数用来求非负整数n的阶乘 %参数n可以为任意的非负整数 %编写日期： 2007-5-2 if n0 %若用户将输入参数误写成负值，则报错 error(输入参数不能为负值！); return; else if n==0 %若n为0,则其阶乘为1 s=1; else s=1; for i=1:n s=s*i; end end end;函数的嵌套调用;【练习】;2. MATLAB程序设计;（1）顺序结构 ;（2）if-else-end 语句;【例】apples=10; cost=apples*25; if apples5 cost=(1-20/100)*cost; end cost cost = 200;if……end语句;if……end语句 ;【例】 if rand(1)0.5 disp(i love you) else disp(i donot love you) end;【例题】编写一个脚本，利用input()函数录入一个由百分数表示的等级数，然后根据下表，用if语句将其转化为一个字母等级： 90%或者更大：A 80%~90%：B 70%~80%：C 60%~70%：D 小于60%：F 请通过反复输入合法和不合法的等级百分数来测试你的程序。;(3) switch……case……end语句* ; (4) try……catch……end* ; (5) for循环语句 ;【例】 for n=1:10 x(n)=sin(n*pi/10); end ? x x = Columns 1 through 7 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns 8 through 10 0.5878 0.3090 0.0000; For循环可按需要嵌套 n=0:1:10; for i=1:11 for j=1:11 y(i)=sin(n(i)); n(j)=n(j)*10; end end;【例】一个简单的for循环示例。 for i=1:10; %i依次取1,2,…10,. x(i)=i; %对每个i值，重复执行由该指令构成的循环体 end; x %要求显示运行后数组x的值。 x = 1 2 3 4 5 6 7 8 9 10 ;例:[tang.m] a=2;b=2; x=-a:0.2:a;y=-b:0.2:b; for i=1:length(y) for j=1:length(x) if x(j)+y(i)1 z(i,j)=0.5457*exp(-0.75*y(i)^2-3.75*x(j)^2-1.5*x(j)); elseif x(j)+y(i)=-1 z(i,j)=0.5457*exp(-0.75*y(i)^2-3.75*x(j)^2+1.5*x(j)); else z(i,j)=0.7575*exp(-y(i)^2-6.*x(j)^2); end end en