MATLAB程序及设计与应用-5.ppt

第五讲 MATLAB数值计算 符号计算 授课教师：田 鹏 数值计算+符号计算 5.1 函数、极限和导数 5.2 空间解析几何 5.3 数列和级数 5.4 数值方法和符号积分 5.5 线形代数 复习 分子物理学绘图 例4.6：利用气体分子运动的麦克斯韦速度分布率，求27°C下氮分子运动的速度分布曲线，并求速度在300-500m/s范围内的分子所占的比例，讨论温度T及分子量?对速度分布曲线的影响。 解： 建模 1.麦克斯韦速度分布率为： 2.考虑到该公式较复杂，建立.m文件。 程序： 1）子程序（mksw.m）： function f=mksw(T,mu,v) R=8.31; %气体常数 k=1.381*10^(-23); %玻尔茨曼常数 NA=6.022*10^23; %阿伏伽德罗数 m=mu/NA; %分子质量 f=4*pi*(m/(2*pi*k*T)).^(3/2) ... .*exp(-m*v.^2./(2*k*T)).*V.^2; %速度分布率 2）主程序： T=input(绝对温度T=); mu=input(气体分子量mu=); vmin=input(速度下限vmin=); vmax=input(速度上限vmax=); v=0:1500; y=mksw(T,mu,v); plot(v,y); hold on; v1=vmin:vmax; %速度分布率 y1=mksw(T,mu,v1); fill([v1,500,300],[y1,0,0],r); trapz(y1); 结果： 当输入: 绝对温度T=300 气体分子量mu=0.028 速度下限vmin=300 速度上限vmax=500 ans = 0.3763 复习．我国古代数学家张丘在“算经”里提出一个世界数学史上有名的百鸡问题：鸡翁一，值钱五，鸡母一，值钱三，幼鸡三，值钱一，百钱买百鸡，问各几何？ 解： 建模 (怎么建？？？) 方法1-1： for x=0:25 y=(100-7*x)/4; if mod(100-x-y,3)==0y=0 fprintf( %g %g %g\n,x,y,100-x-y); end end 方法1-1：运用符号函数 syms+solve syms x y z p=x;q=z; for y=0:33 f1=x+y+z; x=solve(f1-100); f2=5*x+3*y+z/3; z=solve(f2-100); if mod(eval(z),3)==0y+eval(z)=100 fprintf(%g %g %g\n,100-y-eval(z),y,eval(z)); end x=p;z=q; end 方法1-2： disp(鸡翁 鸡母 幼鸡); for x=0:20 for y=0:33 z=100-x-y; if 5*x+3*y+z/3==100 fprintf(%g %g %g\n,x,y,z); end end end 方法1-3： disp(鸡翁 鸡母 幼鸡); for x=0:20 for y=0:33 for z=3:3:99 if x+y+z==1005*x+3*y+z/3==100 fprintf( %g %g %g\n,x,y,z); end end end end 方法2-1： [x,y]=meshgrid(0:20,0:33); t=(find(5*x+3*y+(100-x-y)/3==100)); x(t) y(t) 100-x(t)-y(t) 5.1 函数、极限和导数 一.单变量函数的计算和绘图 例5.1：已知 要求以0.01秒为间隔，绘出y及其导数的曲线. 分析：用diff(y,n)求Dy ， 每求导一次，y的维数减一。 Dy=diff(y)结果为Dy=y1-y2, 故