QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4727|回复: 2
打印 上一主题 下一主题

[求助]数值分析课程设计,求求各位大人啦!!

[复制链接]
字体大小: 正常 放大
hmily777        

1

主题

0

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2007-6-7 15:15 |只看该作者 |倒序浏览
|招呼Ta 关注Ta

一、& p, D: z P" g( R. B# V 要求

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

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

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

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

第三类、三次样条插值法和最小二乘曲线拟合法;1 W0 o7 D! f. W0 s2 W

3.. P) M4 C+ t, z( D: s) u 要求写出完整的实验报告纸质版,算法程序和演示所用的ppt电子版。

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

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
hmily777        

1

主题

0

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

yunizi        

1

主题

3

听众

16

积分

升级  11.58%

该用户从未签到

新人进步奖

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

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;

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2025-8-14 06:44 , Processed in 0.537256 second(s), 67 queries .

回顶部