| 查看: 440 | 回复: 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语言编写 |
回复此楼» 猜你喜欢
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有282人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求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
10