matlab 智能卡pr3,难道是matlab的bug

这部分的所有代码，在matlab中运行没有问题：

20090519224437390.jpg (20.44 KB, 下载次数: 2)

2009-5-20 22:45 上传

但是pr3和pr4的计算结果一样，不对呀：

20090519224346390.jpg (38.04 KB, 下载次数: 0)

2009-5-20 22:45 上传

syms theta a b h w c;

a=100;

b=0.4;

theta=2.5;

c=1;

h=0.16:0.02:0.34;

w=(theta-b).*c./(theta-1);

s= h-(h.*w ./((theta-1).*(w-c)+w.*b));

s(s<0)=0;

sp=((a*b) ./((h-s).*theta)).^(1 ./(1-b)).*((theta-1) ./(theta.*w)).^((theta-1) ./(1-b));

[sp]

p=w.*theta ./(theta-1);

[p]

d=a*sp.^b.*p.^(-theta);

[d]

pr1=(w-c).*a.*sp.^b.*p.^(-theta);

[pr1]

pr2=(p-w).*a.*sp.^b.*p.^(-theta)-h.*sp;

[pr2]

prt1=pr1+pr2;

[prt1]

wr=(c./(theta-1)+theta.*c)./(theta-1);

[wr]

s2=h.*theta./((1+b).*theta-1);

s2(s2<0)=0;

[s2]

sr=((a.*b)./(s2.*theta)).^(1./(1-b)).*((theta-1)./((c+c./(theta-1)).*theta)).^((theta-

1)./(1-b));

[sr]

p2=theta.^2.*c./((theta-1).^2);

dr=a*sr.^b.*p2.^(-theta);

pr3=(wr-c)*a.*sr.^b.*p2.^(-theta)-s2.*sr;

[pr3]

pr4=c./(theta-1)*a.*sr.^b.*p2.^(-theta)-(h-s2).*sr;

[pr2]

prt2=pr3+pr4;

[prt2]

sp2=((a*b) ./(h.*theta)).^(1 ./(1-b)).*((theta-1) ./(theta*c)).^((theta-1) ./(1-b));

[sp2]

p3=c.*theta./(theta-1);

[p3]

prc=(p3-c)*a.*sp2.^b.*p3.^(-theta)-h.*sp2;

[prc]

[本帖最后由 laurence06 于 2009-5-20 22:45 编辑]