| 查看: 535 | 回复: 2 | ||||
[求助]
matlab求解非线性方程组！！！
|
|
本人是matlab小菜，因科研需要需求非线性方程组，方程组如下： AT=0.5 K1=-14 K2=-9.24 K3=-4.97 K4=-15.26 K5=-8.55 K6=-11.33 K7=-2.80 K8=-5.04 K9=-6.77 K10=-7.96 K11=-8.71 K12=-8.74 K13=-14.8 K14=-2.11 K15=-3.74 K16=-4.79 K17=-5.55 K18=-5.73 K19=-5.11 K20=-4.25 K21=-9.2 K22=-10.5 K23=-22.9 K24=-3.9 K25=-12.72 K26=-8.3 K27=-0.8 K28=-1.3 O=10^(P-14) H=10^(-P) N=10^(12.74-2*P) C=10^(13.77-2*P) M=10^(15.28-2*P) NT=N*(1+A/K7+A^2/K8+A^3/K9+A^4/K10+A^5/K11+A^6/K12+O^2/K5+O/K3+O^3/K6) CT=C*(1+A/K14+A^2/K15+A^3/K16+A^4/K17+A^5/K18+A^6/K19+O/K20+O^2/K21+O^3/K22+O^4/K23) MT=M*(1+A/K27+A^2/K28+O/K24+O^3/K26) AT=A*(1+H/K2+N/K7+2*N*A/K8+3*N*A^2/K9+4*N*A^3/K10+5*N*A^4/K11+6*N*A^5/K12+C/K14+2*C*A/K15+3*C*A^2/K16+4*C*A^3/K17+5*C*A^4/K18+6*C*A^5/K19+M/K27+2*M*A/K28) 其中P的范围为[7,14]，求以P为横坐标，NT，CT，MT，AT对P做图的matlab语言编写 |
回复此楼» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有123人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
zmlmf50(金币+50): 2011-09-06 17:19:51
|
A=0.5; K1=-14; K2=-9.24; K3=-4.97; K4=-15.26; K5=-8.55; K6=-11.33; K7=-2.80; K8=-5.04; K9=-6.77; K10=-7.96; K11=-8.71; K12=-8.74; K13=-14.8; K14=-2.11; K15=-3.74; K16=-4.79; K17=-5.55; K18=-5.73; K19=-5.11; K20=-4.25; K21=-9.2; K22=-10.5; K23=-22.9; K24=-3.9; K25=-12.72; K26=-8.3; K27=-0.8; K28=-1.3; P=7:0.01:14; O=10.^(P-14); H=10.^(-P); N=10.^(12.74-2*P); C=10.^(13.77-2*P); M=10.^(15.28-2*P); NT=N.*(1+A/K7+A^2/K8+A^3/K9+A^4/K10+A^5/K11+A^6/K12+O.^2/K5+O./K3+O.^3/K6); CT=C.*(1+A/K14+A^2/K15+A^3/K16+A^4/K17+A^5/K18+A^6/K19+O/K20+O.^2/K21+O.^3/K22+O.^4/K23); MT=M.*(1+A/K27+A^2/K28+O/K24+O.^3/K26); AT=A*(1+H/K2+N/K7+2*N*A/K8+3*N*A^2/K9+4*N*A^3/K10+5*N*A^4/K11+6*N*A^5/K12+C/K14+2*C*A/K15+3*C*A^2/K16+4*C*A^3/K17+5*C*A^4/K18+6*C*A^5/K19+M/K27+2*M*A/K28); plot(P,NT,'r-',P,CT,'g-',P,MT,'-k',P,AT,'-y') legend('NT','CT','MT','AT') |
2楼2011-09-06 14:30:02
3楼2011-09-06 17:21:40
