rogerliu 发表于 2005-12-14 12:04

求matlab解线性方程组的源代码（5-10阶）急用 谢谢

那位高手帮帮忙，我是新手，参考一下就行 谢谢了

happy 发表于 2005-12-14 13:26

回复：（rogerliu）求matlab解线性方程组的源代码（5-...

方程：AX=b，解法：X=A\b<BR><BR>

linqus 发表于 2005-12-14 16:28

对照一般的数值计算书的线性方程部分的算法，用matlab该是很容易实现的。<BR>手头的fortran版的倒是可以找到、matlab版的没有，呵呵

rogerliu 发表于 2005-12-17 10:24

<P>谢谢楼上两位，我想要用高斯列主消元的解法的代码 </P>

happy 发表于 2005-12-17 13:23

回复：（rogerliu）求matlab解线性方程组的源代码（5-...

<FONT style="FONT-SIZE: 14px" color=#555576><FONT style="FONT-SIZE: 14px; LINE-HEIGHT: 17pt" color=#555576>function x=nagauss2(a,b,flag)<BR>% 用途：选列主元Gauss消去法解线性方程组ax=b<BR>% 格式：x=nagauss2(a,b,flag) a为系数矩阵，b为右端列向量，flag若为0，则显示中间过程<BR>%      否则不显示，默认值为0，x为解向量<BR>if nargin&lt;3,flag=0;end<BR>n=length(b);a=;<BR>for k=1:(n-1)<BR>% 选主元<BR>=max(abs(a(k:n,k)));p=p+k-1;<BR>if p&gt;k,<BR>t=a(k,:);a(k,:)=a(p,:);a(p,:)=t;<BR>end<BR>% 消元<BR>a((k+1:n,(k+1):(k+1):(n+1))=a((k+1):n,(k+1):n,(k+1):(n+1))-a((k+1):n,k)/a(k,k)*a(k,(k+1):(n+1));<BR>a((k+1):n,k)=zeros(n-k,1);<BR>if flag==0,a,end<BR>end <BR>% 回代<BR>x=zeros(n,1);<BR>x(n)=a(n,n+1)/a(n,n);<BR>for k=n-1:-1:1<BR>x(k,:)=(a(k,n+1)-a(k,(k+1):n)*x((k+1:n))/a(k,k);<BR>end <BR></FONT></FONT>

青涩果子 发表于 2006-6-21 09:00

用雅可比迭代法求解线性方程组

<P>求VB高手</P>
页: [1]
查看完整版本: 求matlab解线性方程组的源代码（5-10阶）急用 谢谢