quadl('t.^x*sin(x)',1,2,[],[],2)ans =2.1217>> quadl('q.^p*sin(p)',1,2,[],[],2)??? Error using ==> inlineevalError in inline expression ==> q.^p*sin(p)??? Error using ==> mtimesInner matrix dimensions must agree.Error in ==> inline.feval at 34INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);Error in ==> quadl at 64y = feval(f,x,varargin{:}); y = y(:).';>>
3 回答
哆啦的時光機
TA貢獻1779條經驗 獲得超6個贊
@（x）是用來定義無名函數的。以你的函數為例
y=@(x)4*a*sqrt(1-(c/a)^2*sin(x).^2);
這就定義了這樣的一個以x為自變量的函數y，但是與用function來定義的函數不同，它在matlab里沒有函數名，并且只存在于內存中（關掉matlab后就被刪除了）。
之后直接積分就可以
true=quad(y,0,2*pi,1e-8);
你用這個y='4*a*sqrt(1-(c/a)^2*sin(x).^2)';也可以，只是不能用a和c，要把a和c改成對應的數值。比如當a=2和c=1時，你要輸入y='4*2*sqrt(1-(1/2)^2*sin(x).^2)';之后再用true=quad(y,0,2*pi,1e-8);
搖曳的薔薇
TA貢獻1793條經驗 獲得超6個贊
quadl('q.^p*sin(p)',1,2,[],[],2)
這個函數的變量是q而不是p，符號計算時，Matlab自動調用函數確定變量，默認為最靠近x的作為變量
呼喚遠方
TA貢獻1856條經驗 獲得超11個贊
y=@(x)4*3*sqrt(1-(2/3)^2.*sin(x).^2);true=quad(y,0,2*pi)
按上面語句試試
@（x）表示函數以x為自變量
- 3 回答
- 0 關注
- 178 瀏覽
添加回答
舉報
0/150
提交
取消