matlab解一个复杂的一元非线性方程解不出来,求助代码 5
function[yu]=pvf(x,u3)%%一大堆常数参数s=1000;t=25;Sref=1000;Sref2=2500;Tref=25;Uoc0=335/n;Um...
function [yu] =pvf(x,u3)
%%一大堆常数参数
s=1000;t=25;
Sref=1000;Sref2=2500;Tref=25;Uoc0=335/n;Um0=240/n;Isc0=48;Im0=41.67;
a=0.0025;b=0.0005;c=0.00288;
Di=s/Sref*(1+a*(t-Tref));
Du=(1-c*(t-Tref))*log(exp(1)+b*(s-Sref));
Uoc=Uoc0*Du;Um=Um0*Du;Isc=Isc0*Di;Im=Im0*Di;
C2=(Um/Uoc-1)/log(1-Im/Isc);
C1=(1-Im/Isc)*exp(-Um/(C2*Uoc));
Di2=s/Sref2*(1+a*(t-Tref));
Du2=(1-c*(t-Tref))*log(exp(1)+b*(s-Sref2));
Uoc2=Uoc0*Du2;Isc2=Isc0*Di2;
%%%方程
yu=4*C2*Uoc*log(1+(1-x/Isc)/C1)+C2*Uoc2*log(1+(1-x/Isc2)/C1)-u3;
end
主程序:
x0=5;
x=fzero(@pvf,x0,[],250);
fuction上面都是常数参数设定,方程就一个,但是求解出了问题。
求解方程的代码啊啊啊啊 展开
%%一大堆常数参数
s=1000;t=25;
Sref=1000;Sref2=2500;Tref=25;Uoc0=335/n;Um0=240/n;Isc0=48;Im0=41.67;
a=0.0025;b=0.0005;c=0.00288;
Di=s/Sref*(1+a*(t-Tref));
Du=(1-c*(t-Tref))*log(exp(1)+b*(s-Sref));
Uoc=Uoc0*Du;Um=Um0*Du;Isc=Isc0*Di;Im=Im0*Di;
C2=(Um/Uoc-1)/log(1-Im/Isc);
C1=(1-Im/Isc)*exp(-Um/(C2*Uoc));
Di2=s/Sref2*(1+a*(t-Tref));
Du2=(1-c*(t-Tref))*log(exp(1)+b*(s-Sref2));
Uoc2=Uoc0*Du2;Isc2=Isc0*Di2;
%%%方程
yu=4*C2*Uoc*log(1+(1-x/Isc)/C1)+C2*Uoc2*log(1+(1-x/Isc2)/C1)-u3;
end
主程序:
x0=5;
x=fzero(@pvf,x0,[],250);
fuction上面都是常数参数设定,方程就一个,但是求解出了问题。
求解方程的代码啊啊啊啊 展开
展开全部
Matlab解非线性方程组很慢,有时甚至解不出来,建议用 Mathmatic来解
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
北京埃德思远电气技术咨询有限公司
2023-07-25 广告
2023-07-25 广告
整定计算怎么样?这个问题不能一概而论,要根据您的具体情况进行分析。这里简单介绍一下整定计算的特点,供您参考。①整定计算要决定保护的配置与使用,它直接关系到确保系统安全和对重要用户连续供电的问题,同时又和电网的经济指标,运行调度,调试维护等多...
点击进入详情页
本回答由北京埃德思远电气技术咨询有限公司提供
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询