1樓：匿名使用者
function [x,n]=sor(a,b,x0,w,eps,m)if nargin==4
eps= 1.0e-6;
m = 200;
elseif nargin<4
error
return
elseif nargin ==5
m = 200;
endif(w<=0 || w>=2)
error;
return;
endd=diag(diag(a)); %求a的對角矩陣l=-tril(a,-1); %求a的下三角陣u=-triu(a,1); %求a的上三角陣b=inv(d-l*w)*((1-w)*d+w*u);
f=w*inv((d-l*w))*b;
x=b*x0+f;
n=1; %迭代次數while norm(x-x0)>=epsx0=x;
x =b*x0+f;
n=n+1;
if(n>=m)
disp('warning: 迭代次數太多，可能不收斂！');
return;
endend
這個程式無誤。
w=1時，它們倆是等價的。它們得到的結果如果精度要求低可能一樣，但是精度高了後就會不一樣的。
2樓：匿名使用者
哈哈，我剛才也弄了半天。
dw=w*(d+w*l)\b;
最好改為：
dw=(d+w*l)\b*w;
可能效果會好一些。
關於自考求助，求助 關於自考
1。凡事想要好結果都需認真對待 不一定說非要每天都認真學 但至少學校交給你的東西應該掌握吸收 大學都是開放式教學 教學的時間佔的比例不大 你可以根據自己的情況來安排時間 另外值得提醒的是 大學的活動也是很培養人的一個形式 可以多多參加 2。自考分大自考和小自考 全日制助學班舉辦的都是大自考 寬進嚴出...
求助matlab如何利用小矩陣生成大矩陣
使用repmat 函式 來即可實現上述要求，具體 源過程如bai下 1 建立自定義函式repeat matrix a,m,n 保du存為repeat matrix.m function y repeat matrix a,m,n 將給定方zhi陣a擴充為m n的矩陣 n size a,1 rm fl...
求助利用MATLAB進行復雜函式的定積分方法！詳情見補充，多謝了
只要不是不定積分要求得到積分後的表示式，能實際計算出結果的積分就是定積分了。定積分不就是求函式在某區間或領域下的面積或體積麼，所以可以用梯形法之類的方法來近似求解這個面積。比如說步長為h，那麼求出f 0 f h f 2h f nh 則f x 在0 nh的定積分得到的面積就是h f 0 f h f n...