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

Matlab与机电系统仿真（七） ------多项式与微积分 Alfred Bernhard Nobel 我看不出我应得到任何荣誉，我对此也没有兴趣。 阿尔弗雷德·伯纳德·诺贝尔 （1833.10.21－1896.12.10）是瑞典化学家、工程师、发明家、军工装备制造商和炸药的发明者。他曾拥有Bofors军工厂，主要生产军火；还曾拥有一座钢铁厂。在他的遗嘱中，他利用他的巨大财富创立了诺贝尔奖，各种诺贝尔奖项均以他的名字命名。人造元素锘（Nobelium）就是以诺贝尔命名的。 猜一猜，这里有多少个诺贝尔奖获得者？ John Bardeen 约翰·巴丁（John Bardeen，1908年5月23日－1991年1月30日），美国物理学家，因晶体管效应和超导的BCS理论两次获得诺贝尔物理学奖。 电气工程方面的诺贝尔奖 铅笔尖上的诺贝尔奖 7.1 多项式计算 已知：ai，x，求y ---求多项式的值 已知：ai，y，求x ---求多项式的根 已知：x，y，求ai ---求多项式的系数 7.1 多项式求值 %已知：x=3, 求 y x=3; %传统方法 y=x^3+21*x^2+20*x a=[1 21 20 0]; %新方法 y=polyval(a,2) 7.1 多项式求根 %已知：y=0, 求 x %没有传统方法 a=[1 21 20 0]; %新方法 x=roots(a) 7.1 多项式求系数 %已知：y=0, x=[1 2 5 10], 求系数an %没有传统方法 x=[1 2 5 10]; %新方法 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([3 4 5 6 7 8]); b=[7 6 5 4 3]; c=[a,b]; poly(c) 7.1 多项式的乘法 将 乘 (x-7)(x-6)(x-5)(x-4)(x-3) 结果？ %新方法 a=[3 4 5 6 7 8]; 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([3 4 5 6 7],[1,-7]) % y是商，r是余数 7.1 多项式的微分 求 7.1 多项式的微分 求 %传统方法 b=[an,an-1,an-2,…,a1]; c=length(b)-1:-1:1 d=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.33 0.99] [-0.03 -0.12 0.79] [0.1 -0.9 1.9 -0.2] [0.004 0.05 -0.72 1.70 -0.14] [-0.006 0.09 -0.41 0.29 0.86 0.014] 7.1 多项式的拟合 %对于任意一组数据x,y %对x,y进行4阶多项式拟合，多项式系数赋值给a a=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,’