2Â¥: Originally posted by hppdyx at 2014-03-14 19:49:34
Á½²ãÑ­»·£º
     1. µÚÒ»²ãÑ­»·alpha
     2. µÚ¶þ²ãÑ­»·ÀÛ¼Ó·ûºÅ£¨È¡Ç°nÏî¼´¿É£¬±ÈÈçn=50£©

3Â¥: Originally posted by ÎÒÒª·É at 2014-03-14 20:03:23
Çë½Ì...

function hello
%%
beta = 1;
alpha = [0.1 0.3 0.5 0.8];
n = length(alpha);
x = linspace(0, 100, 500);   %xÈ¡Öµ·¶Î§ÉèΪ0µ½100
m = length(x);
k = 30;
z = zeros(m, 1);
%%
for nn = 1 : n
    for mm = 1 : m
        for kk = 0 : k
            zz = - (x(mm))^alpha(n) / gamma(alpha(n)*kk+beta);
            z(mm) = z(mm) + zz;
        end
    end
    figure; plot(x, z);
end

4Â¥: Originally posted by hppdyx at 2014-03-14 20:34:04
function hello
%%
beta = 1;
alpha = ;
n = length(alpha);
x = linspace(0, 100, 500);   %xÈ¡Öµ·¶Î§ÉèΪ0µ½100
m = length(x);
k = 30;
z = zeros(m, 1);
%%
for nn = 1 : n
    for mm = 1 : m
   ...

5Â¥: Originally posted by ÐìÒê at 2014-03-14 21:30:32
ÄãµÄKÖ»ÊÇ´Ó0È¡µ½ÁË30£¬Õâ¸öµÄ³öÀ´µÄ½á¹ûÓ¦¸Ã²»×¼È·µÄ°É£¿º¯ÊýÀïÃæÊÇ0µ½ÎÞÇî...

4Â¥: Originally posted by hppdyx at 2014-03-14 20:34:04
function hello
%%
beta = 1;
alpha = ;
n = length(alpha);
x = linspace(0, 100, 500);   %xÈ¡Öµ·¶Î§ÉèΪ0µ½100
m = length(x);
k = 30;
z = zeros(m, 1);
%%
for nn = 1 : n
    for mm = 1 : m
   ...

7Â¥: Originally posted by ÎÒÒª·É at 2014-03-14 22:22:44
ÇëÎÊÈçºÎ½«ÕâЩ·ÅÔÚÒ»¸öͼƬÀïÃ棬Ҳ¾ÍÊÇÖ»»­Ò»·ùͼ...

function hello
%%
beta = 1;
alpha = [0.1 0.3 0.5 0.8];
n = length(alpha);
x = linspace(0, 100, 500);
m = length(x);
k = 30;
z = zeros(m, n);  % °´Áд洢ËùµÃ½á¹û
%%
for nn = 1 : n
    for mm = 1 : m
        for kk = 0 : k
            zz = - (x(mm))^alpha(nn) / gamma(alpha(nn)*kk+beta);   % Éϸö¸øÄãµÄ³ÌÐòÀïÃæµÄalpha(n)д´íÁË£¬¸ÄΪalpha(nn)
            z(mm, nn) = z(mm, nn) + zz;
        end
    end
end
%% »æͼ
plot(x, z);