Matlab与机电系统仿真（七）+(2).ppt

Matlab与机电系统仿真（七）-多项式与微积分,AlfredBernhardNobel,我看不出我应得到任何荣誉，我对此也没有兴趣。,阿尔弗雷德伯纳德诺贝尔（1833.10.211896.12.10）是瑞典化学家、工程师、发明家、军工装备制造商和炸药的发明者。他曾拥有Bofors军工厂，主要生产军火；还曾拥有一座钢铁厂。在他的遗嘱中，他利用他的巨大财富创立了诺贝尔奖，各种诺贝尔奖项均以他的名字命名。人造元素锘（Nobelium）就是以诺贝尔命名的。,猜一猜，这里有多少个诺贝尔奖获得者？,JohnBardeen,约翰巴丁（JohnBardeen，1908年5月23日1991年1月30日），美国物理学家，因晶体管效应和超导的BCS理论两次获得诺贝尔物理学奖。,电气工程方面的诺贝尔奖,铅笔尖上的诺贝尔奖,7.1多项式计算,已知：ai，x，求y-求多项式的值已知：ai，y，求x-求多项式的根已知：x，y，求ai-求多项式的系数,7.1多项式求值,%已知：x=3,求yx=3;%传统方法y=x3+21*x2+20*xa=121200;%新方法y=polyval(a,2),7.1多项式求根,%已知：y=0,求x%没有传统方法a=121200;%新方法x=roots(a),7.1多项式求系数,%已知：y=0,x=12510,求系数an%没有传统方法x=12510;%新方法a=poly(x),7.1多项式,将(x-7)(x-6)(x-5)(x-4)(x-3)展开？,7.1多项式,将化成乘积的形式？,将乘(x-7)(x-6)(x-5)(x-4)(x-3)结果？,7.1多项式,7.1多项式,将乘(x-7)(x-6)(x-5)(x-4)(x-3)结果？%传统方法a=roots(345678);b=76543;c=a,b;poly(c),7.1多项式的乘法,将乘(x-7)(x-6)(x-5)(x-4)(x-3)结果？%新方法a=345678;b1=1,-7;b2=1,-6;b3=1,-5;b4=1,-4;b5=1,-3;p1=conv(a,b1);p2=conv(p1,b2);c=a,b;poly(c),7.1多项式的除法,将除以(x-7)结果？%传统方法对不起，没想出来！,7.1多项式的除法,将除以(x-7)结果？%新方法y,r=deconv(34567,1,-7)%y是商，r是余数,7.1多项式的微分,求,7.1多项式的微分,求,%传统方法b=an,an-1,an-2,a1;c=length(b)-1:-1:1d=b.*c,7.1多项式的微分,求,%新方法b=an,an-1,an-2,a1,a0;d=polyder(b),7.1多项式的积分,求,%传统方法b=an,an-1,an-2,a1,a0;c=length(b):-1:1;d=b./c,k,7.1多项式的积分,求,%sorry!没有新方法,7.1多项式的拟合,%guess!F(x)=?,-0.330.99,-0.03-0.120.79,0.1-0.91.9-0.2,0.0040.05-0.721.70-0.14,-0.0060.09-0.410.290.860.014,7.1多项式的拟合,%对于任意一组数据x,y%对x,y进行4阶多项式拟合，多项式系数赋值给aa=polyfit(x,y,4);,7.1多项式的拟合,%建立一组数据x=0:0.1:6;y=sin(x);%多项式拟合a=polyfit(x,y,3);z=polyval(a,x);plot(x,y,x,z,r*),7.2数据统计,A=1234214643;%a是一个向量max(a)%最大值min(a)%最小值std(a)%标准差var(a)%均方差,7.2数据差分（微分）,a=123；421；464;%a是一个矩阵diff(a,1,2)%求a的第2维的1阶差分；x,y=gradient(a)%对第2维,第1维求微分x,y=gradient(a,0.1,0.2)%dx=0.1,dy=0.2,7.2数据积分,a=123；421；464;%a是一个矩阵consum(a,2)%对第2维求累计和；trapz(a)%梯形法求积分trapz(x,y)%梯形法求积分cumtrapz(x,y,n)%梯形法求积分，,请大家帮忙！,Matlab与机电系统仿真（八）-符号运算,7.1符号计算,问题：y=x2-x+5a,x=(a-4.5)2+3,求y(a)。计算对象是数值-数值计算；计算对象是符号-符号计算；,7.1符号变量和符号表达式,symsabc%建立符号变量x=a2+5*b-c;%用符号变量组成表达式，赋值给xy=sym(d2+5*e-f);%直接赋值符号表达式z=sym(sqrt(2)+34);%还可以这样赋值哦%注意：用第一种方法，一种创建了abcx四个符号变量；而用第二种方法，只创建了y一个符号变量。,7.1实变量与复变量,symsabreal%建立实变量a、breal(a)ans=asymscdunreal%建立复变量(默认)real(c)ans=1/2*c+1/2*conj(c),7.1符号矩阵,symsabcdx=ab;cdx=sym(ef;gh),7.2符号运算,普通算符：+-*/数组算符：.*./.关系算符：=（没有=）函数：sincostanasinacosatansqrtexplogconjrealimagabs,7.1符号计算,问题：y=x2-4*sqrt(2),x=2+sqrt(2),求y。数值计算？符号计算？哪个精度高？哪个速度快？,7.2运算精度,符号计算的精度可调：digits(n)%设置计算精度，默认值为32vpa(sqrt(2)%采用设置的精度计算2vp