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

>> 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 s>10000, break; end end >> i i = 141 4 MATLAB 语言的函数的基本结构 例：求m,使 function [m,s]=findsum(k) s=0; m=1; while (s> [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;