how to create 1*101 of some integer array?
显示 更早的评论
Hi,
Very basic que: I would like to have 1*101 of -4 array. I have to solve an integral numerically
int_0^x -4 dx;
My code is like:
x= linspace(0,1,0.01);
A = zeros(1,101);
..
cumtrapz(x,-4*A);
But the command prompt shows all zeros for cumtrapz which is certainly not.
If I just use :
cumtrapz(x,-4);
it says matrix dimensions must agree. How can I create this array?
采纳的回答
Mohammad Abouali
2015-4-5
编辑：Mohammad Abouali
2015-4-5
Yo have:
A = zeros(1,101);
and then later you have
cumtrapz(x,-4*A);
all elements of A are zero so -4*A is going to be still zero.
change A=zeros(1,101) to A=ones(1,101)
Also as imageAnalyst said, the third argument for linspace is not the spacing so you need to change the code to something like this
x=linspace(0,1,101);
A=ones(1,101);
cumtrapz(x,-4*A);
By the way, are you sure you want to use cumtrapz()? don't you want the output of trapz()?
2 个评论
Mohammad Abouali
2015-4-6
Both trapz(x,y) and cumtrapz(x,y) integrate between x and y for the entire domain.
Trapz gives you one number. But cumtrapz gives all the values from start to end.
So if x=xStart:dx:xEnd then trapz gives int_xStart^xEnd but cumtrapz will give an array the first one is for int_xStart^{xStart+dx) the second one is for int_xStart^{xStart+2dx} and the n-th number is for int_xStart^{xStart+n*dx}.
the last number of cumtrapz and trapz should be the same.
Both can be useful. just make sure which one is it that you want.
更多回答（1 个）
Image Analyst
2015-4-5
The third argument of linspace is the number of elements, not the increment between elements. Try this:
numberOfElements = 101; % Whatever you want
x= linspace(0, 1, numberOfElements);
A = zeros(1, numberOfElements);
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
