MATLAB解非线性方程组,请帮忙看下代码哪里错了-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 软件培训>>

Matlab软件培训

>>

MATLAB解非线性方程组,请帮忙看下代码哪里错了

MATLAB解非线性方程组,请帮忙看下代码哪里错了

发布:萱草椰子 | 分类:Matlab软件培训

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

functionf=fun19(Y)a1=0.670115;a2=0.038793;a3=0.000398;a4=0.290694;sigma=1.25;x1=970;x2=23.06452;x3=6.666667;x4=231.5385;pper=30;E=1025.13;PE=1.290114;a=0.169354;co2=2466.477;allco2=52034.07;coef1=2.76 ...
免费学术公开课,扫码加入


function f=fun19(Y)
a1=0.670115;a2=0.038793;a3=0.000398;a4=0.290694;
sigma=1.25;
x1=970;x2=23.06452;x3=6.666667;x4=231.5385;
pper=30;
E=1025.13;
PE=1.290114;
a=0.169354;
co2=2466.477;
allco2=52034.07;
coef1=2.763;coef2=2.145;coef3=1.642;coef4=0;
P1=1;P2=6.2;P3=0.3;P4=2.6;
b=1561022;
c1=83.89;c2=70.55;c3=49.56;c4=2.6;
g=0.785379;
F=8078.219;
f(1)=a1 * Y(8)^(1 / sigma - 1) + a2 * Y(9)^(1 / sigma - 1) + a3 * Y(10)^(1 / sigma - 1) + a4 * Y(11)^(1 / sigma - 1)-Y(1);
f(2)=E + sigma / (1 - sigma) * (Y(1) / a - 1) * E-Y(2);
f(3)=Y(14) * Y(2) + pper * allco2-Y(3);
f(4)=Y(16) + pper * coef1-Y(4);
f(5)=Y(17)+ pper * coef2-Y(5);
f(6)= Y(18)+ pper * coef3-Y(6);
f(7)=Y(19) + pper * coef4-Y(7);
f(8)=x1 * (sigma - 1) * (Y(2) / E-1) + X1 * sigma * (Y(3) / b - Y(4)/ c1) + x1-Y(8);
f(9)=x2 * (sigma - 1) * (Y(2)/ E-1) + X2 * sigma * (Y(3)/ b - Y(5)/ c2) + x2-Y(9);
f(10)=x3 * (sigma - 1) * (Y(2)/ E-1) + X3 * sigma * (Y(3)/ b - Y(6)/ c3) + x2-Y(10);
f(11)=x4* (sigma - 1) * (Y(2)/ E-1) + X4 * sigma * (Y(3)/ b - Y(7)/ c4) + x2-Y(11);
f(12)= a1^sigma * Y(16)^(1 - sigma) + a2^sigma * Y(17)^(1 - sigma) + a3^sigma * Y(18)^(1 - sigma) + a4^sigma * Y(19)^(1 - sigma)-Y(12);
f(13)=coef1 * Y(8)+ coef2 * Y(9) + coef3 * Y(10) + coef4 * Y(11)-Y(13);
f(14)= pE + pE / (1 - sigma) * ((Y(12) - g) / g )-Y(14);
f(15)=Y(16) * Y(8)^sigma + Y(17) * Y(9)^sigma + Y(18) * Y(10)^sigma + Y(19) * Y(11)^sigma-Y(15);
f(16)=p1 * Y(15) / f + sigma * p1 + sigma * p1 * Y(8) / x1-Y(16);
f(17)=p2 * Y(15) / f + sigma * p2 + sigma * p2 * Y(9) / x2-Y(17);
f(18)=p3 * Y(15) / f + sigma * p3 + sigma * p3 * Y(10) / x3-Y(18);
f(19)= (Y(2) * Y(14) - Y(8) * Y(16) - Y(9) * Y(17) - Y(10) * Y(18) - pper * (Y(13) - co2)) / Y(11)-Y(19);
Y0=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
Y=fsolve(@fun19,Y0,1)
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
以上是我的代码,解19个方程,19个未知数的方程组,错误提示是
Error in ==> fun19 at 17
f(1)=a1 * Y(8)^(1 / sigma - 1) + a2 * Y(9)^(1 / sigma - 1) + a3 * Y(10)^(1 / sigma - 1) + a4 * Y(11)^(1 / sigma - 1)-Y(1);
我才开始接触MATLAB,不知道为什么错误,请各位好心人帮忙看下,非常感谢!
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-2374201-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
经管之家 人大经济论坛 大学 专业 手机版