黑马求学
2010-07-18, 16:57
我的问题:已知f(x,T)=a*c1/(x^5*(exp(c2/(x*T))-1));
其中,a=1,c1=3.743e-16,c2=1.439e-2;
x的范围:[m,n]; m=820e-9 ; 820e-9
T=[900,1000,1200,1400,1600,1800,2000,2200,2400,2600,2800,3000],
F(T)=int(f(x,y),x,m,n),
求F(T)-T的图形
我的思路:
思路一:先积分,得到F(T)关于T的表达式,在做图,
思路二:代入T值,积分得到对应的F(T)值,再经曲线拟合作图
程序代码:
syms x
c1=3.743e-16;
c2=1.439e-2;
m=820e-9 ;
n=880e-9 ;
T=[900,1000,1200,1400,1600,1800,2000,2200,2400,2600,2800,3000];
f=c1./(x^5*(exp(c2./(x*T))-1));
F=int(f,x,m,n);
T=double(T);
F=double(F);
plot(F,T)
怎么也得不到正确的结果,只用MATLAB能解决吗?该怎么解决呢?请帮我看看,诚心求教,谢谢。
我在附件中附上了正确的图形结果
其中,a=1,c1=3.743e-16,c2=1.439e-2;
x的范围:[m,n]; m=820e-9 ; 820e-9
T=[900,1000,1200,1400,1600,1800,2000,2200,2400,2600,2800,3000],
F(T)=int(f(x,y),x,m,n),
求F(T)-T的图形
我的思路:
思路一:先积分,得到F(T)关于T的表达式,在做图,
思路二:代入T值,积分得到对应的F(T)值,再经曲线拟合作图
程序代码:
syms x
c1=3.743e-16;
c2=1.439e-2;
m=820e-9 ;
n=880e-9 ;
T=[900,1000,1200,1400,1600,1800,2000,2200,2400,2600,2800,3000];
f=c1./(x^5*(exp(c2./(x*T))-1));
F=int(f,x,m,n);
T=double(T);
F=double(F);
plot(F,T)
怎么也得不到正确的结果,只用MATLAB能解决吗?该怎么解决呢?请帮我看看,诚心求教,谢谢。
我在附件中附上了正确的图形结果