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

MVH · 发表于 2005-7-28 17:48

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

function [x1]=yakebi(A,b)
%Jacobi mehtod to solve the linear equation
%chinamaker@dytrol.com,12.18
num=length(A);
M=zeros(num,num);
x=zeros(num,1);
N=-A;
for i=1:num
M(i,i)=A(i,i);
N(i,i)=0;
x(i,1)=1;
end
err=1;
B=inv(M)*N;
g=inv(M)*b;
while abs(err)>0.00001
xx=B*x+g;
err=norm(xx-x,1);
x=xx;
end
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 · 发表于 2006-3-28 13:53

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

girl-058 · 发表于 2006-3-28 13:58

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

happy · 发表于 2006-3-28 14:57

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

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

girl-058 · 发表于 2006-3-29 16:41

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

happy · 发表于 2006-3-29 19:27

上面的程序不就是了吗？

girl-058 · 发表于 2006-3-30 17:29

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

happy · 发表于 2006-3-30 17:41

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

girl-058 · 发表于 2006-5-18 16:33

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

girl-058 · 发表于 2006-5-19 16:56

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

复制代码

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

happy · 发表于 2006-5-19 18:10

试一下lsqcurvefit

girl-058 · 发表于 2006-5-21 12:36

居然没有人知道哦

账号		自动登录	找回密码
密码			我要加入

[编程技巧] maker系列－Jacobi迭代法求线性方程组

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。