Integration by matlab

2011 5 30
1 个回答
回答已采纳
4 次查看（30 天）

0 个投票

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

1 个投票

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 个评论
显示 1更早的评论 隐藏 1更早的评论

reem
reem 2011-5-30
Thanks a lot Matt
when I do that
This problem appear
??? Error using ==> plot
Conversion to double from function_handle is not possible.
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 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by