maker系列－Jacobi迭代法求线性方程组

MVH

哈哈，可能本例有点简单，不过这可能就是我的考试题目，我就练一下手，呵呵，好久没有看书了。希望大家不要见笑。

1. function [x1]=yakebi(A,b)
   
2. %Jacobi mehtod to solve the linear equation
   
3. %chinamaker@dytrol.com,12.18
   
4. num=length(A);
   
5. M=zeros(num,num);
   
6. x=zeros(num,1);
   
7. N=-A;
   
8. for i=1:num
   
9. M(i,i)=A(i,i);
   
10. N(i,i)=0;
    
11. x(i,1)=1;
    
12. end
    
13. err=1;
    
14. B=inv(M)*N;
    
15. g=inv(M)*b;
    
16. while abs(err)>0.00001
    
17. xx=B*x+g;
    
18. err=norm(xx-x,1);
    
19. x=xx;
    
20. end
    
21. x1=x;%Jacobi method result

复制代码

实例演示
>> A=[8,1,3,2;2,9,2,3;1,2,7,1;3,2,1,9];
>> b=[1,3,4,6]';
>> x=yakebi(A,b)
x =
-0.2363
0.0496
0.4939
0.6795

>> inv(A)*b
ans =
-0.2363
0.0496
0.4939
0.6795

girl-058

请问jacobi旋转法matlab编程怎么编啊

girl-058

请问jacobi旋转法matlab编程怎么编啊

happy

请问jacobi旋转法matlab编程怎么编啊

没听说过有jacobi旋转法这种方法
Jacobi方法是一种用平面旋转矩阵所构成的正交相似变换将对称矩阵化为对角型的方法，所以又叫做平面旋转法

girl-058

就是平面旋转法啊
那怎么编程呢
用它来求矩阵的特征值及特征向量啊

happy

上面的程序不就是了吗？

girl-058

不懂啊
一个是特征值一个是求线性方程一样吗？
学得不好所以不明白

happy

找本算法书看看Jacobi法算特征值是怎么回事吧，搞懂了算法写程序不难

girl-058

在鱼塘中投放N0尾鱼苗，随着时间的增长，尾数将减少而每尾的重量将增加。
（1）设尾数n(t)的相对减少率为常数；由于喂养引起的每尾鱼重量的增加率与鱼表面积成正比，由于消耗引起的每尾鱼重量的减少率一重量本身成正比，分别建立尾数和每尾鱼重量的微分方程，并求解
（2）用控制网眼的办法不捕小鱼，到时刻才开始捕捞，捕捞能力用尾数的相对减少量|N/n|表示，记作E，即单位时间捕获量是En(t)，问如何选择T和E使从T开始的捕捞量最大。
帮忙一下谁会数学建模啊

girl-058

1. w(t)=[(k1/k2)*(1-exp(-(k2*t)/3))].^(3);
   
2. 已知t=[1 2 3 4 5 6 7 8];
   
3. w=[150 200 300 370 450 507 600 756];

复制代码

怎么用MATLAB求出k1,k2啊
帮忙一下啊

happy

试一下lsqcurvefit

girl-058

居然没有人知道哦