clear
clf
close all
nn=[21,50,150];
dt=0.1;
for i=1:length(nn)
n=nn(i);
t=0:dt:n;
for j=1:length(t)
f(j)=1/(n+1); %如果有阶乘改这一句
for k=0:n
f(j)=f(j)*(t(j)-k);
end
end
figure
plot(t,f);
title(n);
xlabel('t');
ylabel('f(t)');
end
我做出来的值比较大,我值计你的表达式里面有点小问题。
如果有阶乘,将上面有注明的那句改成
f(j)=1/factorial(n+1);
n=21; % it can be changed as 50 150
t=1:0.01:n;
m=length(t);
A=zeros(m,n+1);
for i = 1:n+1
A(:,i)= (t-i+1)';
end
B=prod(A');
B=B';
plot(t,B./(n+1));
f(t)好像一直等于0 啊