| 查看: 1632 | 回复: 4 | ||
[求助]
求用牛顿迭代法解非线性方程组的matlab程序 已有2人参与
|
方程组形式如图所示，希望高手帮忙，谢谢~ D~(P4F`R%ZU2}EJ}BZ}Z38A.jpg |
回复此楼» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有103人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限差分法化之后化简为线性方程组
已经有50人回复
- 求教 如何用matlab解带约束的非线性方程组
已经有4人回复
- matlab 非线性微分方程求解
已经有3人回复
- 用matlab求解一个非线性方程组的解
已经有4人回复
- 大家能给看看这个方程属于什么类型，应该怎么求解？
已经有10人回复
- 用fortran程序遗传算法解非线性方程组
已经有7人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 对流扩散方程迭代求解发散
已经有11人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助matlab解非线性规划。求代码。
已经有13人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 请问牛顿法和高斯牛顿法有什么不同？非常感谢
已经有6人回复
- matlab求解非线性方程组
已经有16人回复
- 求一个mathmetica的程序
已经有10人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
2楼2014-05-08 21:11:35
895701777
木虫 (著名写手)
进城农民工
- 应助: 36 (小学生)
- 金币: 1384.8
- 散金: 706
- 红花: 12
- 帖子: 1547
- 在线: 287.9小时
- 虫号: 2052402
- 注册: 2012-10-10
- 性别: GG
- 专业: 新闻理论
3楼2014-05-09 09:23:57
reko34
木虫 (正式写手)
- 应助: 117 (高中生)
- 金币: 4794.2
- 散金: 22
- 红花: 12
- 帖子: 780
- 在线: 417.4小时
- 虫号: 3055206
- 注册: 2014-03-15
- 性别: GG
- 专业: 机械测试理论与技术
【答案】应助回帖
感谢参与，应助指数 +1
|
白天没机会碰MATLAB，急的话先试，有问题的话等我晚上试，但注释的2个东西要LZ根据实际情况自己定。 MATLAB解N维牛顿法有优势的，我之前用EXCEL编的，可能会有语法搞起来。 a=?; b=?; x=?; %迭代初始向量 Jacob=zeros(n,n); y=zeros(n,1); while (stopcriterion) %迭代停止条件自己定 pdx=5*x.^4+3*x.^2+1-b; for i=1:n Jacob(i,i)=pdx(i); end y(1)=x(1)^5+x(1)^3+x(1)+a-b*(x(1)+x(2)); for i=2:n-1 y(i)=fx(x(i),x(i-1),x(i+1),a,b); Jacob(i-1,i)=-b; Jacob(i,i-1)=-b; end y(n)=x(n)^5+x(n)^3+x(n)+a-b*(x(n)+x(n-1)); Jacob(n-1,n)=-b; Jacob(n,n-1)=-b; dx=y/Jacob; x=x+dx; end %向量x是结果 function fx(x,x0,x1,a,b) fx=x^5+x^3+x+a-b*(x+x0+x1); end |
4楼2014-05-09 12:49:16
reko34
木虫 (正式写手)
- 应助: 117 (高中生)
- 金币: 4794.2
- 散金: 22
- 红花: 12
- 帖子: 780
- 在线: 417.4小时
- 虫号: 3055206
- 注册: 2014-03-15
- 性别: GG
- 专业: 机械测试理论与技术
5楼2014-05-09 20:24:51
