数学建模社区-数学中国

标题: [求助]数值分析课程设计，求求各位大人啦！！ [打印本页]

---

作者: hmily777    时间: 2007-6-7 15:15
标题: [求助]数值分析课程设计，求求各位大人啦！！

一、 要求

1．来源参考：可以是数学分析和线性代数课程中的相关题目；可以通过查找文献和利用网上资源等自主选定一个项目；

2．内容要求：以下三类方法至少应包含其中两类

第一类、求解线性代数方程组的列主元Gauss消去法和Gauss-Seidel迭代法以及改进的平方根法和追赶法；

第二类、求解非线性方程的二分法、牛顿法及其改进方法；

第三类、三次样条插值法和最小二乘曲线拟合法； 6 u ]2 O8 \8 [1 u9 g5 v' j. F G; G

3. 要求写出完整的实验报告纸质版，算法程序和演示所用的ppt电子版。

真的是很迫切需要，小女子无以为报，只能跪地求饶啦

---

作者: hmily777    时间: 2007-6-7 15:19
啊，忘记说了，请用MATLAB做 ，再次衷心的谢谢啦！！！！！

---

作者: yunizi    时间: 2008-4-22 22:10

这个问题不是很难，我可以告诉你其中一部分，其余的你可以靠自己

Gauss-Seidel迭代

function[x,k]=Gauss(A,b,x0,N.emg)

n=length(A);

x1=zeros(n,1);x2=zeros(n,1);

x1=x0；

r=max(abs(b-A*x1));

k=0;

while r>emg

for i=1:n

sum=0;

for j=1:n

if j>i

sum=sum+A(i,j)*x1(j);

eldeif j<i

sum=sum+A(i,j)*x2(j);

end

end

x2(i)=(b(i)-sum)/A(i,i);

end

r=max(abs(x2-x1));

x1=x2;

k=k+1;

if k>N

disp("迭代失败")；

return;

end

end

x=x1;

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5