• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

matlab:追赶法(crout分解)工具箱

原作者: [db:作者] 来自: [db:来源] 收藏 邀请
 1 function x=chase (a,b,c,f) 
 2 % the method of chaase*******************************
 3 %  a, b, c,分别是是方程组的下对角,主对角,上对角的三条
 4 %  alpha,beta, gama分别是Crout分解的主和次对角线
 5 %    /beta(1)    0    ...             ...        0  0     \
 6 %    |alpha(2)   beta(2)  ...          ...       0  0     | 
 7 %  L=|0          alpha(3) ...beta(3)                      |
 8 %    |                                                    |
 9 %    |                                                    |
10 %    \                                   alpha(3) beta(n) /
11 %********************************************************
12 beta(1)=b(1);gama(1)=c(1)/beta(1);y(1)=f(1)/beta(1);
13 for i=2:length(b)
14 alpha(i-1)=a(i-1);beta(i)=b(i)-alpha(i-1)*gama(i-1);
15 if i<=length(b)-1
16     gama(i)=c(i)/beta(i);
17 end
18 y(i)=(f(i)-alpha(i-1)*y(i-1))/beta(i);
19 end
20 x(length(b))=y(length(b));
21 for i=length(b)-1:-1:1
22 x(i)=y(i)-gama(i)*x(i+1);
23 end

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
Matlab强行终止程序执行发布时间:2022-07-22
下一篇:
MATLABGUI图片添加背景发布时间:2022-07-22
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap