24小时热门版块排行榜    

CyRhmU.jpeg
查看: 173  |  回复: 2
当前主题已经存档。

spring180

银虫 (小有名气)

[交流] 【求助】请教高手帮忙解决问题所在,重金酬谢!

E=linspace(4000,200000,1600);
S=int('exp(-E/(8.314*T)','T',507,512);
S
       我本意是想求活化能E在4000-200000范围内,所对应的各个不同的积分值S,但是却出现问题,请问如何解决?非常感谢!
回复此楼
科研更多的是代表一种精神
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)

★ ★ ★ ★ ★
spring180(金币+5,VIP+0):谢谢 我先看看 9-9 10:40
可以用下面的修改版,但是速度很慢
syms E T
Ek=linspace(4000,200000,1600);
S=int(exp(-E/(8.314*T)),T,507,512)
Sn=double(subs(S,5000))
for k=1:length(Ek);
    Sn(k)=double(subs(S,Ek(k)));
    k
end
plot(Ek,Sn)
2楼2009-09-09 10:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kuhailangyu(金币+2,VIP+0):欢迎常来 9-9 11:04
spring180(金币+25,VIP+0): 10-16 09:05
下面的程序计算速度快一些:
Ek=linspace(4000,200000,1600);
F=inline('exp(-E./(8.314*T))','T','E');
for k=1:length(Ek);
    Sd(k)=quad(F,507,512,1e-6,[],Ek(k));
    k
end
plot(Ek,Sd)
3楼2009-09-09 10:40:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 spring180 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见