matlab求解物理函数
需要用matlab计算,求T下的值,谢谢大佬。
clear;
q=1;
k=8.6174e-5;
BN=0.45;
A=110;
VR=45;
JF=IF/S;
IF=20;
s=112;
S=(112*25.4/1000)^2;
f=JF*(BN+k*T/q*ln(JF/(A*T^2)))-A*T^2*e^(-q*BN/(k*T)*VR);
ans=solve(f==0,T)
返回小木虫查看更多
今日热帖
需要用matlab计算,求T下的值,谢谢大佬。
clear;
q=1;
k=8.6174e-5;
BN=0.45;
A=110;
VR=45;
JF=IF/S;
IF=20;
s=112;
S=(112*25.4/1000)^2;
f=JF*(BN+k*T/q*ln(JF/(A*T^2)))-A*T^2*e^(-q*BN/(k*T)*VR);
ans=solve(f==0,T)
返回小木虫查看更多
直接数值解 用fzero函数
题主的问题是没有按软件函数的要求编程。完整的代码为

clear;
syms T
q=1;
k=8.6174e-5;
BN=0.45;
A=110;
VR=45;
IF=20;
S=(112*25.4/1000)^2;
JF=IF/S;
f=JF*(BN+k*T/q*log(JF/(A*T^2)))-A*T^2*exp(-q*BN/(k*T)*VR);
ans=vpasolve(f==0,[T])
这里,用vpasolve计算比solve更精确。
计算结果如下。
123.jpg
,