网站大量收购闲置独家精品文档,联系QQ:2885784924

线性代数科学应用举例.ppt

  1. 1、本文档共51页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
x 0 1 2 3 4 y -27 0 21 0 -75 例 7 多项式插值与拟合 求: (1) 过这五个点作一个四次多项式函数 (2) 请根据这五个点,拟合一个二次多项式函数 下表给出了平面坐标系中五个点的坐标。 并求 x=5 时的函数值 p4(5)。用 MATLAB 绘制多项式函 数 p4(x) 的曲线、已知点及插值点 (5, p4(5))。 并用 MATLAB 绘制 p2(x) 的曲线及已知的五个点。 其中矩阵: 解:(1) 根据已知条件,把五个点的坐标值分别代入 四次多项式函数,可以得到如下线性方程组: 系数矩阵 A 的行列式为范德蒙 (Vandermonde) 行列式, 且五个坐标点的横坐标各不相同,则该行列式不等于零, 所以方程组有唯一解。 MATLAB 程序:(ma8) x=[0;1;2;3;4]; % 输入已知点坐标 y=[-27;0;21;0;-75]; A=[x.^0, x.^1, x.^2, x.^3, x.^4]; % 构造 vandermonde 矩阵 a=A\y; % 得到适定方程组的唯一解 a 运行程序,得到 a(1)=-27, a(2)=12, a(3)=26, a(4)=-12, a(5)= 1. 把五个点的坐标值分别代入二次多项式函数,可以得 到如下线性方程组: 其中, (2) 多项式拟合要解一个超定方程 该方程组有三个未知数,但有五个方程,进一步分析 可以得到该方程组无解,即不存在一个二次多项式曲 线刚好能过已知的五个点。MATLAB 软件提供了一个 利用最小二乘法解决超定方程组解的方法。求系数的 公式也是 a = A\y,以找到一条二次曲线来近似地描述 已知 5 个点的变化情况。 对比插值和拟合的曲线如下图 用平面坐标系中的一个闭合图形来描述刚体,用一个矩 阵 X 来表示它。X 的一列表示刚体一个顶点的坐标。为 了使图形闭合,X 的最后一列和第一列相同;为了实现 刚体的平移运算,给矩阵 X 添加元素值都为 1 的一行, 使 X 为 3×n 矩阵。 若有矩阵: 则可以证明,矩阵 Y1 是刚体 X 沿 x 轴正方向平移 c1, 沿 y 轴正方向平移 c2 后的结果;矩阵 Y2 是刚体 X 以坐 标原点为中心逆时针转动 t 弧度的结果。 例 8 刚体的平面运动 x 0 4 6 10 8 5 3.5 6.1 6.5 3.2 2 0 y 0 14 14 0 0 11 6 6 4.5 4.5 0 0 实 例 用下列数据表示字母 A: 对 A 进行以下平面运动, 并绘制移动前后的图形。 (1) 向上移动 15, 向左移动 30; (2) 逆时针转动π /3; (3) 先逆时针转动3π /4, 然后向上平移 30, 向右平移 20。 解 构造刚体矩阵 X,平移矩阵及转动矩阵。 M1= 1 0 0 0 1 0 -30 20 1 M2= 1 0 0 0 1 0 20 30 1 R1= 0 0 1 0 0 R2= 0 0 1 0 0 MATLAB 程序(ma9) ? X=[0,4,6,10,8,5,3.5,6.1,6.5,3.2,2,0; 0,14,14,0,0,11,6,6,4.5,4.5,0,0; ones(1,12)]; % 构造刚体矩阵 X M1=[1,0,-30; 0,1,15; 0,0,1]; % 构造平移矩阵 M1 Y1=M1*X; % 计算平移结果 ? fill(Y1(1, :), Y1(2, :), ‘red’); % 绘制平移后刚体 Y1 ? plot(X(1, :), X(2, :)); % 绘制原来刚体 X hold on axis equal ? R1=[cos(pi/3), -sin(pi/3), 0; sin(pi/3), cos(pi/3), 0; 0,0,1]; % 构造转动矩阵 R1 Y2=R1*X;

您可能关注的文档

文档评论(0)

sandaolingcrh + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档