一元函数的图像显示

tongxin28

此函数是y关与n（其实是x）的一元函数，n（既x）取值0:1000,就是运算量有点大，小弟刚学matlab，不太清楚要是想显示此函数的图像，该如何输入此函数，只知道最后用plot(x,y),请高手不吝指教一下。
为难各位了，表达式如何输入啊？

[ 本帖最后由 ChaChing 于 2010-7-27 23:59 编辑 ]

ericlin

我觉得应该分解开吧，比如说先循环算出n×（n－1）（n－2）……这样？

wei343

<P>我想要是把分子和分母都单独运算应该行的，尝试用符号运算然后再将符号转化为数值估计也可以，作的好会比直接用数值好。</P>

branch

function doing()
y=[];
for n=13:1000
midst=ones(1,13);
for i=1:13
midst(1,i)=(n-i)*10/(i*n);
end
middle=cumprod(midst);
a=middle(13)/sum(middle);
y=[y,a];
end
n=13:1000;
plot(n,y)

tongxin28

感谢楼上的 ，我验证一下！
感谢branch，经验证很ok。帮小弟解释其中的几条主要语句吧。

[ 本帖最后由 ChaChing 于 2010-7-28 00:15 编辑 ]

ericlin

^_^，由见branch使用cumprod这个函数了。
顶一下。

branch

呵呵，cumprod这比较好用啊

midst=ones(1,13);可省去
for i=1:13
midst(1,i)=(n-i)*10/(i*n);循环赋值，产生一个矩阵元素为每相乘数
end
middle=cumprod(midst);一维行累乘，如对a=[1 2 3 4]行累乘得[1 2 6 24]
a=middle(13)/sum(middle);得到总公式

[ 本帖最后由 ChaChing 于 2010-7-28 00:18 编辑 ]