哪位高人能够用C语言求解下列方程组??50分相送!

来源:百度知道 编辑:UC知道 时间:2024/05/19 19:35:52
0.70965*a=9.5*(T-c)+K*((t-0.5)/log((T-5)/(T-t)))+0.1032935(1-0.0045(T-25))*a;
K*0.937732*((t-0.5)/log((T-5)/(T-t)))=b(0.6973t+130.9435);
K=1/(0.001/237+1/α);
α=0.002345*pow(b,1.4)+1087.536*pow(K*((t-0.5)/log((T-5)/(T-t))),0.7)/pow(b,0.7);
其中a,b,c是自变量,T,t是因变量,K与α可以消元。
要到时候只要我输入a,b,c的数值能求出T,t的值就可以,谢谢各位!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!用matlab也行!a=1000,b=0.045,c=25

你能给出a,b,c具体值吗?我试了一下用a=1,b=2,c=3代入,结果根本解不出方程,所以我认为无法实现给出任意a,b,c都能求出解,所以请你给出a,b,c具体值,这样用遗传算法或者GLP算法都能一下子求出解来。

按你的新方程的要求,a=1000,b=25计算结果如下:

T:92.1106
b:3.162e-010
k:0
误差:3.14932e-005

请按任意键继续. . .

我的邮箱:hu_hu605@163.com

晕呀...看了.都晕了........不解.不解.不解...

如果你告诉a b c 的值,一下子就算出来了。

function hhh
global a b c
a=input('a=')
b=input('b=')
c=input('c=')
a=1;b=2;c=3;
fsolve(@myfun,[100 10 100 10])
function f=myfun(x)
global a b c
T=x(1);t=x(2);alpha=x(3);K=x(4);
f1=9.5*(T-c)+K*((t-0.5)/log((T-5)/(T-t)))+0.1032935*(1-0.0045*(T-25))*a-0.70965*a;
f2=K*0.937732*((t-0.5)/log((T-5)/(T-t)))-b*(0.6973*t+130.9435);
f3=K-1/(0.001/237+1/alpha);
f4=0.002345*power(b,1.4)+1087.536*power(K*((t-0.5)/log((T-5)/(T-t))),0.7)/power(b,0.7)-alpha;
f=[f1;f2;f3;f4]

看吐了,楼主牛人,牛方程

晕,犯不着为了点分这么伤脑细胞啊。你这问题首先还得分析这里面的变量关系,得得到数学模型。太烦了!

你可以