原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
山东大学数值计算实验报告3
山东大学 计算机科学与技术 学院
数值计算 课程实验报告
学号： 姓名： 鸿武 班级：计机14.1
实验题目：
1、Computer Problems 3.1(P152)
1.1 画出散点图 （plot 函数）
1.2 用一次多项式拟合。给出拟合函数，确定矛盾方程组的系数矩阵、右端
项。正规方程的方法得到结果，并画图显示该拟合函数
1.3 用二次多项式拟合，画图
1.4 用三次多项式拟合，画图
1.5 前面几个曲线拟合，假设对应的矛盾方程组为Ax ≌b，分别计算A 和
T
A A 的条件数，并比较
2、例题3.2 （图3.1 ）的21 个点
1.1 画出散点图
1.2 用一次多项式拟合。给出拟合函数，确定矛盾方程组的系数矩阵、右端
项。正规方程的方法得到结果，并画图显示该拟合函数
1.3 用二次多项式拟合，画图
1.4 用双曲线拟合 （选做）
实验学时：4 实验日期：2016 年11 月12 日
实验目的：最小二乘法数据拟合
硬件环境： PC 机
软件环境： MATLAB、JAVA、C++
实验步骤与内容：
1、
1.1 画出数据散点图如下：
代码如下所示：
t = [0.0 1.0 2.0 3.0 4.0 5.0];
y = [1.0 2.7 5.8 6.6 7.5 9.9];
scatter(t,y)
1.2 一次多项式拟合，拟合函数：
(, ) = 1 + 2
可以写成如下形式：
1 0 1
1 1 2.7
= 1 2 [1] ≅ 5.8 =
1 3 2 6.6
1 4 7.5
[1 5] [9.9]
用正规方程的方法来编写MATLAB 函数如下所示。
test12.m
clc,clear
t = [0.0 1.0 2.0 3.0 4.0 5.0];
b = [1.0 2.7 5.8 6.6 7.5 9.9];
t = t;
b = b;
A = ones(6,2);
A(:,2) = t;
AA = A*A;
Ab = A*b;
L = Cholesky(AA);
y = inv(L)*A*b;
x = inv(L)*y
scatter(t,b)
hold on;
syms m yy
ezplot(yy==x(1)+x(2)*m,[-1,10])
hold on;
B = A*A;
fprintf(%s\n,1 范数下矩阵AA 的条件数)
Acond(B)
fprintf(%s\n,2 范数下矩阵AA 的条件数)
Acond2(B)
fprintf(%s\n,无穷范数下矩阵AA 的条件数)
Acond3(B)
Cholesky.m
function [ B ] = Cholesky(A )
%UNTITLED2 此处显示有关此函数的摘要
% 楚列斯基分解
B = A;
n = length(A);
for k = 1:n
B(k,k) = sqrt(B(k,k));
for i = k+1:n
B(i,k) = B(i,k)/B(k,k);
end
for j = k+1:n
for i = k+1:n
B(i,j)
文档评论（0）