Integration by matlab

1 回表示 (過去 30 日間)
reem
reem 2011 年 5 月 30 日
Good morning,
Please,I have a problem in integration
please anyone help me and I will appreciate it so much.
My program
N=400;
f1=1; % Frequency of the cosine wave
FS=200; % Sampling Frequency
t=1:400;
n=0:N-1;
X= @(Q)cos((2*pi*f1*n/FS)+Q);
M=quadgk(X,0,(2*pi));
disp(M)
plot(t,M)
grid;
title('cosinewave ');
xlabel('Time index');
ylabel('cosine value');
Kind regards

採用された回答

Matt Fig
Matt Fig 2011 年 5 月 30 日
Use this instead:
M=quadv(X,0,(2*pi));
%
%
%
In the future, you need to be more specific than just saying that you have an error. You should have said this:
This code produces an error, and here is the error message:
??? Error using ==> plus
Matrix dimensions must agree.
Error in ==> @(Q)cos((2*pi*f1*n/FS)+Q)
Error in ==> quadgk>evalFun at 363
fx = FUN(x);
Error in ==> quadgk>f1 at 381
[y,too_close] = evalFun(tt);
Error in ==> quadgk>vadapt at 275
[fx,too_close] = f(x);
Error in ==> quadgk at 214
[q,errbnd] = vadapt(@f1,interval);
>>
  3 件のコメント
Matt Fig
Matt Fig 2011 年 5 月 30 日
Not when I ran the code. Did you change more than the line I showed you to change? The only line I changed from what you originally pasted was the one line I showed, and I got a plot with no problems.
reem
reem 2011 年 5 月 30 日
Sorry,there is no any problem
In the future, I will be more specific
Thank you so so much, you helped me a lot now and before, and I will not forget your help forever
You are the best teacher for me
I benefited from you a lot

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeMathematics についてさらに検索

タグ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by