MATLAB图片为空 5
formatshortsymspvhgmusigmaz2z1wcaikymu1sigma1t1t2q1q2qp=35;%a=0.1;w=25;c=15;k=200;mu=...
format short
syms p v h g mu sigma z2 z1 w c a i k y mu1 sigma1 t1 t2 q1 q2 q
p=35;
%a=0.1;
w=25;
c=15;
k=200;
mu=100;
sigma=20;
g=10;
h=2;
v=10;
mu1=0;
sigma1=1;
t1=normpdf(z1,mu1,sigma1);
t2=normpdf(z2,mu1,sigma1);
T1=int(t1);
T2=int(t2);
i=0.4;
for a=0.75:0.01:0.3;
eq1=sym([(p+g)*(1-a)-w]/[(1-a)*(p-v+h+g)]-T1);
eq2=sym([(p+g)*(1-a+a*i)-w]/[(1-a+a*i)*(p-v+h+g)]-T2);
z1=solve(eq1);
z2=solve(eq2);
q1=sym((mu+sigma*z1)./(1-a));
q2=sym((mu+sigma*z2)./(1-a+a*i));
%axis([0,0.3,-100,100]);
plot(a,q1,'-.rv',a,q2,'--ks')
hold on;
end 展开
syms p v h g mu sigma z2 z1 w c a i k y mu1 sigma1 t1 t2 q1 q2 q
p=35;
%a=0.1;
w=25;
c=15;
k=200;
mu=100;
sigma=20;
g=10;
h=2;
v=10;
mu1=0;
sigma1=1;
t1=normpdf(z1,mu1,sigma1);
t2=normpdf(z2,mu1,sigma1);
T1=int(t1);
T2=int(t2);
i=0.4;
for a=0.75:0.01:0.3;
eq1=sym([(p+g)*(1-a)-w]/[(1-a)*(p-v+h+g)]-T1);
eq2=sym([(p+g)*(1-a+a*i)-w]/[(1-a+a*i)*(p-v+h+g)]-T2);
z1=solve(eq1);
z2=solve(eq2);
q1=sym((mu+sigma*z1)./(1-a));
q2=sym((mu+sigma*z2)./(1-a+a*i));
%axis([0,0.3,-100,100]);
plot(a,q1,'-.rv',a,q2,'--ks')
hold on;
end 展开
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询