上限可变积分的求解

wxwChina

       目的：求解随时间和频率变化的暂态能量趋势。公式如图所示，程序如附件。       思路是对一被积函数进行积分，其中有两个自变量，只是对其中时间t自变量进行积分，频率f自变量按照t=0：50e-6：0.02s进行变化，然后输出三维图。
       求解上限可变积分的时候，用int积分总是出错，不知道怎么改。结果就改成利用for循环进行，但是很慢。看看各位吧友能不能给指点如何修改，利用[f,t] = meshgrid(0:500,0:50e-6:0.02); mesh(f,b,result),xlabel('f'),ylabel('t'),zlabel('E');之类的，不用for循环。

wxwChina

                               
登录/注册后可看大图