原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
數值分析在生活中的应用举例及Matlab实现
最小二乘法，用MATLAB实现
1. 数值实例
下面给定的是乌鲁木齐最近1个月早晨7：00左右（新疆时间）的天气预报所得到的温度，按照数据找出任意次曲线拟合方程和它的图像。下面用MATLAB编程对上述数据进行最小二乘拟合。
2008年10月26~11月26
天数
1
2
3
4
5
6
7
8
9
10
温度
9
10
11
12
13
14
13
12
11
9
天数
11
12
13
14
15
16
17
18
19
20
温度
10
11
12
13
14
12
11
10
9
8
天数
21
22
23
24
25
26
27
28
29
30
温度
7
8
9
11
9
7
6
5
3
1
下面用MATLAB编程对上述数据进行最小二乘拟合
2、程序代码
x=[1:1:30];
y=[9,10,11,12,13,14,13,12,11,9,10,11,12,13,14,12,11,10,9,8,7,8,9,11,9,7,6,5,3,1];
a1=polyfit(x,y,3) %三次多项式拟合%
a2= polyfit(x,y,9) %九次多项式拟合%
a3= polyfit(x,y,15) %十五次多项式拟合%
b1= polyval(a1,x)
b2= polyval(a2,x)
b3= polyval(a3,x)
r1= sum((y-b1).^2) %三次多项式误差平方和%
r2= sum((y-b2).^2) %九次次多项式误差平方和%
r3= sum((y-b3).^2) %十五次多项式误差平方和%
plot(x,y,*) %用*画出x,y图像%
hold on
plot(x,b1, r) %用红色线画出x,b1图像%
hold on
plot(x,b2, g) %用绿色线画出x,b2图像%
hold on
plot(x,b3, b:o) %用蓝色o线画出x,b3图像%
3、数值结果
不同次数多项式拟合误差平方和为：
r1=67.6659
r2=20.1060
r3=3.7952
r1、r2、r3分别表示三次、九次、十五次多项式误差平方和。
4、拟合曲线如下图
二、 线性方程组的求解（ 高斯-塞德尔迭代算法 ）
1、实例： 求解线性方程组（见书P233页）
记Ax=b, 其中
任取初始值，进行迭代。
2、用C语言实现，源程序如下：
#includestdio.h
#includestdlib.h
struct Node
{
double value; // 矩阵非零元素值
int index; // 元素的索引(序号或列标)
};
struct Matrix
{
struct Node *data; // 记录元素数据
int total_elem; // 矩阵非零元素个数
int total_ln; // 矩阵的总行数
int total_col; // 矩阵的总列数
};
// 采用MSR法实现系数矩阵的存储
void Create_Matrix( struct Matrix matrix )
{
int index = 0, i = 0;
double value = 0.0;
printf( 请输入线性方程组系数矩阵的非零元素个数及总行数、总列数: );
scanf( %d%d%d, matrix.total_elem, matrix.total_ln, matrix.total_col );
matrix.data = (struct Node *)malloc( sizeof(struct Node)*(matrix.total_elem+1) );
printf( 请依次输入矩阵的主对角线元素及每行的非零非对角线元素(行优先)\n );
for( i = 0; i = matrix.total_elem; i++ )
{
if( i == matrix.total_ln )
value = 0;
else
scanf( %lf, value );
matrix.data[i].value = value;
}
printf( \n系数矩阵的信息相关如下(Quote为添加项,表示序列号): );
printf( \nValue: );
for( i = 0; i matrix.total_elem + 1; i++ )
printf( %3.0lf, matrix.data[i].value );
printf( \nQuote
文档评论（0）