求高手用matlab求数值解求高手帮忙求如下方程组的数值解c2*(T-t)*[a+1/2*b(T-t)+1/3c*(T^2+Tt+t^2)]=0(c*t^2+b*t+a)*[(c1+c2+h/m+k/m^2)**exp(m*t)-(h+kt)/m-k/m^2-c2*(T-t)]=0其中:c1、c2、a、b、c、h、k、m均已知,需求T,t的
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/06 01:48:22
求高手用matlab求数值解求高手帮忙求如下方程组的数值解c2*(T-t)*[a+1/2*b(T-t)+1/3c*(T^2+Tt+t^2)]=0(c*t^2+b*t+a)*[(c1+c2+h/m+k/m^2)**exp(m*t)-(h+kt)/m-k/m^2-c2*(T-t)]=0其中:c1、c2、a、b、c、h、k、m均已知,需求T,t的
求高手用matlab求数值解
求高手帮忙求如下方程组的数值解
c2*(T-t)*[a+1/2*b(T-t)+1/3c*(T^2+Tt+t^2)]=0
(c*t^2+b*t+a)*[(c1+c2+h/m+k/m^2)**exp(m*t)-(h+kt)/m-k/m^2-c2*(T-t)]=0
其中:c1、c2、a、b、c、h、k、m均已知,需求T,t的数值解.
请求高手,万分感谢!
求高手用matlab求数值解求高手帮忙求如下方程组的数值解c2*(T-t)*[a+1/2*b(T-t)+1/3c*(T^2+Tt+t^2)]=0(c*t^2+b*t+a)*[(c1+c2+h/m+k/m^2)**exp(m*t)-(h+kt)/m-k/m^2-c2*(T-t)]=0其中:c1、c2、a、b、c、h、k、m均已知,需求T,t的
你没有给出参数的具体值,我只好编了.
function hh
fsolve(@myfun,[100,10])
function f=myfun(x)
T=x(1);t=x(2);
c1=1;c2=9;a=8;b=3;c=4;h=5;k=6;m=7;
f(1)=c2*(T-t)*(a+1/2*b*(T-t)+1/3*c*(T^2+T*t+t^2));
f(2)=(c*t^2+b*t+a)*((c1+c2+h/m+k/m^2)*exp(m*t)-(h+k*t)/m-k/m^2-c2*(T-t));
结果:
ans =
66.7804 1.0423
或者这样,更简单:
c1=1;c2=9;a=8;b=3;c=4;h=5;k=6;m=7;
f=@(x) [c2*(x(1)-x(2))*(a+1/2*b*(x(1)-x(2))+1/3*c*(x(1)^2+x(1)*x(2)+x(2)^2));
(c*x(2)^2+b*x(2)+a)*((c1+c2+h/m+k/m^2)*exp(m*x(2))-(h+k*x(2))/m-k/m^2-c2*(x(1)-x(2)))]
x=fsolve(f,[100,10])
T=x(1),t=x(2)
结果:
x =
66.7804 1.0423
T =
66.7804
t =
1.0423