原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE 6
1. LU分解：
[L U]=lu(A);
2. 追赶法
function x=zhuiganfa(A,b)
[n,n]=size(A);
for i=1:n
if(i==1)
l(i)=A(i,i);
y(i)=b(i)/l(i);
else
l(i)=A(i,i)-A(i,i-1)*u(i-1);
y(i)=(b(i)-y(i-1)*A(i,i-1))/l(i);
end
if(in)
u(i)=A(i,i+1)/l(i);
end
end
x(n)=y(n)
for j=n-1:-1:1
x(j)=y(j)-u(j)*x(j+1);
end
数值试验：
n=101;
a= 12 1 0………………………………..….0
1 12 1 ………………………………..…0
0 1 12 1 ………………………..…0
0 0 1 12 1 ……………………0
1
1 12
关键是如何定义上述矩阵：
n=101;
c1=ones(1,n-1);
a1=diag(c1,-1); 这个-1说明行位置-1
c2=12*ones(1,n);
a2=diag(c2);
c3=ones(1,n-1);
a3=diag(c3,1);
a=a1+a2+a3;
3. 拉格朗日插值
function yh=lage(x,y,xh)
n=length(x);
m=length(xh);
yh=zeros(1,m);
c1=ones(n-1,1);
c2=ones(1,m);
for i=1:n
xp=x([1:i-1 i+1:n]);
yh=yh+y(i)*prod((c1*xh-xp*c2)./(x(i)-xp*c2));
end
end
x=[11,12];
y=[2,4];
xh=[11.75];
lage(x,y,xh)
ans =
3.5000
4 最小二乘法
最小二乘的xi和yi为：
xi
1953
1964
1982
1990
2000
yi
5.82
6.95
10.08
11.34
12.66
要拟合的函数为：y=a+bx-cxy 注意不是多项式
编程函数为：
function z = erchen(x,y)
x1=ones(5,1);
A=[x1,x,-x.*y]; 注意点乘
z=A\y; 注意左除
a=z(1);
b=z(2);
c=z(3);
end
输入：
?? x=[1953 1964 1982 1990 2000];
?? y=[5.82 6.95 10.08 11.34 12.66];
?? erchen(x,y)
ans =
2.9456 =a
-0.0014 =b
-0.0005 =c
最小二乘的xi和yi为：
xi
0
0.25
0.5
0.75
1
yi
1
1.284
1.6487
2.1170
2.7183
要拟合的函数为：y=a+bx+cx2 是多项式
编程函数为：
function z = erchen2(x,y)
x1=ones(5,1);
A=[x1,x,x.^2];
z=A\y;
a=z(1);
b=z(2);
c=z(3);
end
输入：
?? x=[0 0.25 0.5 0.75 1.00];
?? y=[1.00 1.284 1.6487 2.1170 2.7183];
?? erchen2(x,y)
ans =
1.0051
0.8642
0.8437
最小二乘多项式拟合的简单函数方法：
?? x=[0 0.25 0.5 0.75 1.00];
?? y=[1.00 1.284 1.6487 2.1170 2.7183];
?? P=polyfit(x,y,2) 要拟合成4次，则2改成4就可以了
P =
0.8437 0.8642 1.0051 注意此内置函数输出的结果c,b,a是反的
5 复合辛普森公式求解积分
先定义函数：
function v=f(x)
v=sin(x); “若定义有除数要点除，分母有0时要特殊定义”
end
定义程序：
function I=fsps(f,a,b,
文档评论（0）