- 1、本文档共51页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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;
您可能关注的文档
- 物理选修教案-磁场及电场基础知识汇编.ppt
- 线性变换及其矩阵表示.ppt
- 线性代数教案.ppt
- 物理选修教案恒定电流系统知识合编.ppt
- 物流和物流信息系统.ppt
- 线性代数总结.ppt
- 物流优化技术绪论.ppt
- 物质结构第二章分子结构和性质.ppt
- 误差理论-创新设计和实践.ppt
- 宪法是国家根本大法说课.ppt
- 2020版 沪科技版 高中生物学 必修2 遗传与进化《第4章 生物的进化》大单元整体教学设计[2020课标].docx
- 情绪价值系列报告:春节消费抢先看-国证国际证券.docx
- 精品解析:北京市东直门中学2023-2024学年高二下学期3月阶段性考试(选考)物理试题(解析版).docx
- 2020版 沪科技版 高中生物学 必修2 遗传与进化《第4章 生物的进化》大单元整体教学设计[2020课标].pdf
- 2020版 沪科技版 高中生物学 选择性必修1 稳态与调节《第1章 人体的内环境和稳态》大单元整体教学设计[2020课标].pdf
- 2020版 沪科技版 高中生物学 选择性必修1 稳态与调节《第1章 人体的内环境和稳态》大单元整体教学设计[2020课标].docx
- 液冷盲插快接头发展研究报告-全球计算联盟.docx
- 精品解析:北京市东直门中学2023-2024学年高二下学期3月阶段性考试(选考)物理试题(原卷版).docx
- 精品解析:北京市东直门中学2024届高三考前练习数学试卷(解析版).docx
- 2020版 沪科技版 高中生物学 选择性必修1 稳态与调节《第2章 人体的神经调节》大单元整体教学设计[2020课标].docx
文档评论(0)