Matlab 版 (精华区)

发信人: hahn (有奇☆闭关造文中), 信区: Matlab
标  题: [合集] 如何画出这样的一个函数
发信站: 哈工大紫丁香 (Wed Oct 11 19:20:27 2006), 站内

────────────────────────────────────────
  logocde (吉人自有天相)           于  (Thu Jun  1 10:38:08 2006)  说道:

附件中是一个概率密度函数：
n1＝11；n2＝33；n3＝33；
c12＝1；c13＝2；L＝3；
下面是我用符号运算作的 ：
syms x t c m n L
g=c^(2*n*L)*(1-x)^(2*m*L-1)/((1-x)^2*m+n*c^2)^((m+n)*L)+c^(-m*L)*(1-x)^(2*n*L-
1)/((1-x)^2*n+m/c^2)^((m+n)*L);
h1=subs(g,[m,n,c,L],[11,33,1,3]);
h2=subs(g,[m,n,c,L],[11,33,2,3]);
k1=int(h1);
k2=int(h2);
f=h1*(subs(k1,x,1)-k1)+h2*(subs(k2,x,1)-k2);
Const=4*(gamma((m+n)*L))^2*m^(2*L*m)*n^(2*L*n)/(gamma(L*m)^2*gamma(L*n)^2);
P=subs(Const,[m,n,L],[11,33,3]);
fh=P*f;
ezplot(fh,[0,1])

P是前边的一大陀系数，我尝试了几个n1，n2，n3的取值，发现对于给定的n1＝11；n2＝3
3；n3＝33；p巨大无比，而x在0，1之间时f值又极小，matlab显示为0，都超出了正常的数
据范围。
我把p中的指数部分移到了f中计算，仍然不能解决问题，因为gamma函数那部分的值太大了
。
请大家都帮我看一看吧，想不出啥办法了