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

Matlab与机电系统仿真 七 多项式与微积分 AlfredBernhardNobel 我看不出我应得到任何荣誉 我对此也没有兴趣 阿尔弗雷德 伯纳德 诺贝尔 1833 10 21 1896 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 x 3 21 x 2 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 x 2 x 5a x a 4 5 2 3 求y a 计算对象是数值 数值计算 计算对象是符号 符号计算 7 1符号变量和符号表达式 symsabc 建立符号变量x a 2 5 b c 用符号变量组成表达式 赋值给xy sym d 2 5 e f 直接赋值符号表达式z sym sqrt 2 3 4 还可以这样赋值哦 注意 用第一种方法 一种创建了abcx四个符号变量 而用第二种方法 只创建了y一个符号变量 7 1实变量与复变量 symsabreal 建立实变量a breal a ans asymscdunreal 建立复变量 默认 real c ans 1 2 c 1 2 conj c 7 1符号矩阵 symsabcdx ab cd x sym ef gh 7 2符号运算 普通算符 数组算符 关系算符 没有 函数 sincostanasinacosatansqrtexplogconjrealimagabs 7 1符号计算 问题 y x 2 4 sqrt 2 x 2 sqrt 2 求y 数值计算 符号计算 哪个精度高 哪个速度快 7 2运算精度 符号计算的精度可调 digits n 设置计算精度 默认值为32vpa sqrt 2 采用设置的精度计算 2v