[2018年最新整理]5-MATLAB(M).pptVIP

f([153.06,2.0079,0.1855],12) ans = 18.2063 f([153.06,2.0079,0.1855],11) ans = 21.9171 f([153.06,2.0079,0.1855],11.5) ans = 19.9757 即在很短时间内喝完3瓶酒，且喝完酒后的11.5小时之内驾车就会违反国家标准。 问题 二 由于酒是在较长一段时间内喝的，即不能忽略饮酒的时间，对于此问我们通过假设对喝酒的过程进行平均分割，m（小时）时间内共饮了n次，分割成每m/n（小时）喝一次，，相当于在相同时间段内连续饮酒，我们可以根据 m： 饮酒总时间 n : 饮酒次数 D=75000; m=2; n=40; dd=D/n; k1=2.0715; k2=0.1810; k=k1/(k2-k1); T=m/n; a=dd*exp(-k1*T); b=dd*k*(exp(-k1*T)-exp(-k2*T)); l=2; for l=2:1:n-1 b=b*exp(-k2*T)+(a+dd)*k*(exp(-k1*T)-exp(-k2*T)); a=(a+dd)*exp(-k1*T); l=l+1; End a0=(a+dd)/490; b0=b/490; C2(n)(t)=b0*exp(-k2*t)+a0*k1*(exp(-k1*t)-exp(-k2*t))/(k2-k1) 用迭代的方法算出在第n次喝酒的c2(n)(t)，然后算出对应的t 值， 求得： t = 10.81小时,酒精含量会下降到19.9781； 即在连续喝两个小时酒的情况下，且喝完后的10.81小时之内驾车就会违反国家标准。 * 例：用循环求解 s=0;for i=1:100 s=s+i;end s=0; i=1;while (i=100) s=s+i; i=i+1;end sum(1:100) ans = 5050 例：用循环求解求最小的 m s=0; m=0; while (s=10000), m=m+1; s=s+m; end, [s,m] % 求出的 m 即是所求 ans = 10011 141 转移结构 例：用循环求解求最大的 m s=0; for i=1:10000 s=s+i; if s10000, break; end end i i = 141 4 MATLAB 语言的函数的基本结构 例：求m,使 function [m,s]=findsum(k) s=0; m=1; while (sk) s=s+m; m=m+1; end [m1,s1]=findsum(145323) m1 = 539 s1 = 145530 inline 函数和匿名函数 inline 函数，可以免去文件 f=inline(‘sin(x.^2+y.^2)’,’x’,’y’) MATLAB 7.0 5 二维图形绘制 二维图形绘制基本语句 构造向量: 例：选项为红色点划线且每个转折点用五角星表示 ‘r-.pentagram’ 例： x=[-pi : 0.05: pi]; % 以 0.05 为步距构造自变量向量 y=sin(tan(x))-tan(sin(x)); % 求出各个点上的函数值 plot(x,y) plot(x,y,r-.pentagram) 6 三维图形绘制 三维曲线绘制 stem3, fill3, bar3 等 例：参数方程 t=0:.1:2*pi; % 构造 t 向量，注意下面的点运算 x=t.^3.*sin(3*t).*exp(-t); y=t.^3.*cos(3*t).*exp(-t); z=t.^2; plot3(x,y,z), grid % 三维曲线绘制 stem3(x,y,z); hold on; plot3(x,y,z), grid 三维曲面绘制 一般曲面绘制 mesh( )绘制网格图,surf( )绘制表面图。 其他函数，光照下 surfl( ), surfc( )，瀑布型waterfall( ) 等高线绘制 contour( ), contour3( ) 例:Butterworth 滤波器 [x,y]=meshgrid(0:31); n=2; D0=200; D=sqrt((x-16).^2+(y-16).^2); % 求距