设为首页收藏本站联系我们

声振论坛 › 展示 ›科学计算› 算法编程 › 查看内容

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

2015-11-3 01:19| 发布者: aspen| 查看: 1637| 评论: 11|原作者: MVH|来自: 声振论坛

摘要: 哈哈，可能本例有点简单，不过这可能就是我的考试题目，我就练一下手，呵呵，好久没有看书了。希望大家不要见笑。 function =yakebi(A,b) %Jacobi mehtod to solve the linear equation %chinamaker@dytrol.com,1 ...

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

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

收藏分享邀请

上一篇：常用算法设计方法下一篇：matlab lms 算法实例

最新评论

引用 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: 居然没有人知道哦

查看全部评论(11)

相关分类

|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-10 04:17 , Processed in 0.088412 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部