laurence06
2009-05-11, 17:03
通过下面的程序，我本意是希望b取不同的值：0，0.05，0.10，0.15一直到0.55，其他参数值不变来计算f1，f2一直到f6；但是不知道是否for循环发生问题，居然不计算。请知道的朋友帮忙指出问题。另外，我还希望看f1随b变化的趋势图，怎么用matlab实现？
--------------------------------------------------------------------------
for b=0:0.05:0.6
syms theta a h w c;
a=100;
theta=2.5;
c=1;
h=0.25;
w=1.6;
f1=((a*b)/(h*theta))^(1/(1-b))*((theta-1)/(theta*w))^((theta-1)/(1-b));
f2=w*theta/(theta-1);
f3=a*f1^b*f2^(-theta);
f4=((1-b)/(theta-1))*a*((a*b)/(h*theta))^(b/(1-b))*((theta-1)^2/(theta*c*(theta-b)))^((theta-b)/(1-b));
f5=(f2-w)*a*f1^b*f2^(-theta)-h*f1;
[f1 f2 f3 f4 f5]
>> f6=f4+f5;
[f6]
end;
--------------------------------------------------------------------------
--------------------------------------------------------------------------
for b=0:0.05:0.6
syms theta a h w c;
a=100;
theta=2.5;
c=1;
h=0.25;
w=1.6;
f1=((a*b)/(h*theta))^(1/(1-b))*((theta-1)/(theta*w))^((theta-1)/(1-b));
f2=w*theta/(theta-1);
f3=a*f1^b*f2^(-theta);
f4=((1-b)/(theta-1))*a*((a*b)/(h*theta))^(b/(1-b))*((theta-1)^2/(theta*c*(theta-b)))^((theta-b)/(1-b));
f5=(f2-w)*a*f1^b*f2^(-theta)-h*f1;
[f1 f2 f3 f4 f5]
>> f6=f4+f5;
[f6]
end;
--------------------------------------------------------------------------