数学物理建模作业讲评.ppt

对函数sin(x)+cos(x)取任意多个数据点，求拟合6阶多项式，并图示拟合情况，并对图形进行修饰。 用符号运算应该怎么做？ %计算位移和速度 function [s,v]=velocity(X0,V0,A,T); syms x x0 v0 a t y1=dsolve(D2x=a,Dx(0)=v0,x(0)=x0,t); %求解满足的微分方程，得到位移 y2=diff(y1,t); %对位移作一次导数得到速度 s=subs(y1, {x0,v0,a,t},{X0,V0,A,T}); %变量替换得到位移和速度的值 s=double(s); v=subs(y2, {x0,v0,a,t},{X0,V0,A,T}); v=double(v); 求点P(2,-4,5)关于直线的垂足和对称点。要求编M函数文件来实现，并绘制这些直线和点。 %解法1：求p点垂足和对称点 syms x y z p=[2,-4,5]; q=[x,y,z]; g1=(x-2)/3=(y+3)/7; g2=(y+3)/7=(z-1)/(-3); g3= dot(q-p,[3 7 -3]); f=solve(g1,g2,g3) C=double([f.x,f.y,f.z]) S=2.*C-p %求中点和对称点 %入口参数 P0 已知点坐标 % x0 直线上点坐标 % k 直线方向矢量 function [C,S]=points(P0,x0,k); %x0=[2,-3,1]; %k=[3 7 -3]; %P0=[2,-4,5]; P1=sym([x,y,z]); %符号表达式，待求对称点坐标 P2=(P0+P1)/2; %符号表达式，垂足点坐标 Online=(P2-x0)./k; %符号表达式，垂足点代入直线方程 Perpen=sum((P2-x0).*k); %符号表达式，垂直条件代入 f1=sym(strcat(char(Online(1)),=,char(Online(2)))); %垂足点在线条件1 f2=sym(strcat(char(Online(1)),=,char(Online(3)))); %垂足点在线条件2 f3=sym(strcat(char(Online(2)),=,char(Online(3)))); %垂足点在线条件3 f4=sym(strcat(char(Perpen),=0)); %连线垂直于直线条件 s=solve(f1,f2,f4); %求条件方程，注意这里4个条件只有3个是独立的，求解3个未知数 S=double([s.x,s.y,s.z]); C=double((S+P0)/2); 回顾一下泰勒级数展开 泰勒(Taylor)展开定理 设f(z)在区域D:|z-z0|R内解析，则在D内f(z)可展开成泰勒级数 其中 一些常见的泰勒级数展开公式 回顾一下罗朗级数展开 罗朗级数展开定理 设函数f(z)在D:R1|z-z0|R2内解析，则在此圆环内f(z)必可展开成罗朗级数 其中 应用两种级数展开复变函数的时候首先确定： 展开所在的区域属于哪一个环域 若内部区域解析，则采用泰勒展开 若内部有奇点或不解析，则采用罗朗展开 绘制复变函数 的图像，并验证 其泰勒展开和罗朗展开 泰勒展开 罗朗展开 罗朗展开 图形的绘制 (1)直接绘制 z=4*cplxgrid(40); cplxmap(z,1./(z+eps*(abs(z)==1)-1)./ (z+eps*(abs(z)==2)-2),10*pi); 图形的绘制 (2)泰勒展开 z=2*cplxgrid(30); z1=z; z1(abs(z1)=1)=NaN; w1=1; u1=1; p1=1; q1=1; for k=1:100 u1=u1.*z1; w1=w1+u1; p1=p1.*z1/2; q1=q1+p1; end cplxmap(z1,w1+q1/2) colorbar 图形的绘制 (3)罗朗展开 z=3*cplxgrid(30); z2=z; z2(abs(z2)=1|abs(z2)=2)=NaN; w2=1./z2; u2=1./z2; p2=1; q2=1; for k=1:100 u2=u2./z2; w2=w2+u2; p2=p2.*z2/2; q2=q2+p2; end cplxmap(z2,-w2-q2/2)