Matlab 版 (精华区)

发信人: hahn (有奇☆闭关造文中), 信区: Matlab
标  题: [合集] 再问关于符号运算的问题
发信站: 哈工大紫丁香 (Wed Oct 11 19:19:28 2006), 站内

────────────────────────────────────────
  logocde (吉人自有天相)           于  (Wed May 31 21:21:56 2006)  说道:

刚刚接触matlab的符号运算，自己找了半天也没觉得哪里有问题，大家帮帮找找错，谢谢

！
要画出附件中的函数的图形，因为还要对它作积分，所以用符号运算。
其中c，n，L等等都可以用常数代换，变量只有一个。这是一个概率密度函数
我写的程序如下：
syms x t c m n L
g=c^(2*n)*(1-x)^(2*m-1)/((1-x)^2*m+n*c^2)^(m+n)+c^(-m)*(1-x)^(2*n-1)/((1-x)^2*

n+m/c)^(m+n);
h1=subs(g,[m,n,c],[33,22,1]);
h2=subs(g,[m,n,c],[33,22,2]);
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,100])

可是画不出来不说，P居然是inf，我取m，n和L的值没问题，不知道哪里有问题。
请大家帮忙看一下，谢谢！！！


────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed May 31 22:10:12 2006)  说道:

给int函数一个积分上下限
【 在 logocde (吉人自有天相) 的大作中提到: 】
: 刚刚接触matlab的符号运算，自己找了半天也没觉得哪里有问题，大家帮帮找找错，谢谢
: ！
: 要画出附件中的函数的图形，因为还要对它作积分，所以用符号运算。
: ...................



────────────────────────────────────────
  logocde (吉人自有天相)           于  (Wed May 31 22:17:44 2006)  说道:

求不定积分，然后在代入积分上下限求差不也一样吗？
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 给int函数一个积分上下限





────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed May 31 22:18:48 2006)  说道:

不定积分没意思吧,得出来的数值也是不确定的
变上限积分还有点意思
你考虑把积分限确定吧
【 在 logocde (吉人自有天相) 的大作中提到: 】
: 求不定积分，然后在代入积分上下限求差不也一样吗？




────────────────────────────────────────
  logocde (吉人自有天相)           于  (Wed May 31 22:29:26 2006)  说道:

int好像不能作变上限积分，所以我作不定积分，得出的k1，k2是符号表达式，然后我再把
积分限代入的。
我试了用int作变上限积分，int（f，t，1）是不行的，请问有什么问题吗？
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 不定积分没意思吧,得出来的数值也是不确定的
: 变上限积分还有点意思
: 你考虑把积分限确定吧