清空筛选条件
清空筛选条件

Numerical analysis please help

1 次查看（过去 30 天）
francis
francis 2013-7-26
if true
i=1;
for q=0:0.05:4;
x(i)=q;
s(i)=quad('sin(x.^2)',0,x(i));
c(i)=quad('cos(x.^2)',0,x(i));
i=i+1;
end
fprintf('\n The value of s(x)=%5.3f',s(i-1))
fprintf('\n The value of c(x)=%5.3f\n',c(i-1))
plot(x,s,x,c)
xlabel('x')
ylabel(' C(x) and S(x)')
grid
figure
plot(c,s)
xlabel('c'),ylabel('s')
end
  4 个评论
显示 2更早的评论
Mike Hosea
Mike Hosea 2013-7-31
I got no error. What version of MATLAB are you using?
Andrew Reibold
Andrew Reibold 2013-7-31
I am using 2012b I believe, and I also received no errors. I did get two sexy looking graphs though.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

dpb
dpb 2013-7-26
Replace the loop structure
i=1;
for q=0:0.05:4;
x(i)=q;
s(i)=quad('sin(x.^2)',0,x(i));
c(i)=quad('cos(x.^2)',0,x(i));
i=i+1;
end
with
q=[0:0.05:4]';
s=accumarray([1:length(q)]',q,[],@(x) quad('sin(x.^2)',0,x));
c=accumarray([1:length(q)]',q,[],@(x) quad('cos(x.^2)',0,x));
If you want to get even more clever, define the functions as function handles... :)
I didn't pursue the initial too much; you should at a minimum preallocate the vectors s and c

类别

在 Help Center 和 File Exchange 中查找有关 Electromechanical 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by