http://imageshack.us/photo/my-images/210/schapery.png/
想請問一下這個積分式用matlab該怎麼做呢
a(x)為一分段函數
我的問題是
1.分段函數有辦法用int算出解析解嗎
2.裡面的積分式上限為變數t 和 tau
所以可以先用int做不定積分 然後最後再用數值積分運算嗎??
-------------------------------
程式碼
a=@(x)1.*(x<540)+(237.4*(0.033*x).^-1.891).*(x>=540&x<900)+(237.4*(30).^-1.891).*(x>=900&x<1800)+(237.4*(30-0.033*(x-1800)).^-1.891).*(x>=1800&x<2160)+1.*(x>=2160&x<3600);
rt = int(a,x,0,t);
rt2= int(a,x,0,tau);
k=(rt-rt2).^1.8;
s=quad(k,0,3600);
可是就是一堆錯誤QQ
有人知道怎麼改嗎......
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.155.119