数学物理建模作业讲评课件.pptVIP

对函数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);;解法2： function [Pd,Pp]=fp(P,Vp,V) Pv=P-Vp; L=sqrt(dot(V,V)), C=V*dot(Pv,V)/L*2 D=Pv-C; Pd=Vp+C; Pp=P-2*D; [Pd,Pp]=fp([2,-4,5],[2,-3,1],[3,7,-3]);解法3： k=[3,7,-3];p=[2,-4,5]; A=[7 -3 0;0 3 7;3 7 -3];b=[23 -2 -37]; x=inv(A)*b m=2*x-p;回顾一下泰勒级数展开 泰勒(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;