求助MATLAB编程问题

对q(j)*e(k-j)求和（j=0,1,2,...k）其中q(j)=[1-(1+gama)/j]*q(j-1)

初值是什么呢？

不好意思忘写了，q(0)=1

e和gamma分别指的是什么呢

e是一个误差表示 初值应该是0

gama自己来选定范围吗

gama就是一个任意的常数

这样可以吗？
add_sum.m

1. 
   
2. function y=add_sum()
   
3. gama=input('Please input gama:');
   
4. k=input('Please input k:');
   
5. q=zeros(1,k+1);
   
6. e=zeros(1,k+1);
   
7. q(1)=1;
   
8. e(1)=0;
   
9. y=0;
   
10. for i=2:k+1
    
11.     q(i)=(1-(1+gama)/(i-1))*q(i-1);
    
12.     e(i)=q(i)-q(i-1);
    
13. end
    
14. q
    
15. e=e(end:-1:1)
    
16. y=sum(q.*e);
    
17. end
    

复制代码

result:

1. 
   
2. y=add_sum
   
3. Please input gama:0.5
   
4. Please input k:6
   
5. q =
   
6.     1.0000   -0.5000   -0.1250   -0.0625   -0.0391   -0.0273   -0.0205
   
7. 
   
8. e =
   
9.     0.0068    0.0117    0.0234    0.0625    0.3750   -1.5000         0
   
10. 
    
11. y =
    
12.     0.0205
    

复制代码

题目是(1+gama)/i)，为什么却:( 写成(1+gama)/(i-1))？达人指点哈！

因为i是从2开始的，在matlab中数组索引最小是1，不是0。所以得往后递推一个

HOHO~~太谢谢了！！！:lol :victory: