请帮忙修改一个优化程序
这是一个MATLAB的M文件!在保证Pf最小的情况下,用matlab优化Rf,Hpf,Hcf的值?
%***************************************************************
% This function is used to caculate the measurement of solar
% collector.
%***************************************************************
function =SoclCal(ta,ra,pr,pc,r,yt,pout)
%***************************************************************
% Tempretures
%***************************************************************
ts=348;
tm=324;
ti=(ta+tm)/2;
tt=(ts-ti)/ta;
%**************************************************************
% Physical Coefficient
%**************************************************************
rm=1.093;
k=2.76e-2;
g=9.8;
cp=1.005e3;
rf=24.3e-6;
v=16.96e-6;
pm=(rm^8*k^5*g^15)/(ra^15*cp^15*rf^(5/3)*v^(7/3));
%**************************************************************
% Coefficient
%**************************************************************
roc=pi*pr/cos(2*pi/45);
hoc=(pc/tt)*(pout^11/((0.018*yt)^11*pm*r^16))^(1/15);
%**************************************************************
% Independent Variable
%**************************************************************
R=;
Hp=;
%**************************************************************
% Grid Paremeters
%**************************************************************
x=R;
y=Hp;
=meshgrid(x,y);
P=roc.*X.^2+hoc.*Y.^(0.80)./X;
%***************************************************************
% Get mininum value of P
%***************************************************************
a=min(P);
Pf=min(a);
=size(P);
zb=0;
for i=1:s
for j=1:t
if P(i,j)==Pf
zb=1;
break;
end
end
if zb==1
break;
end
end
Rf=X(i,j);
Hpf=Y(i,j);
Hcf=(Pf-pi*pr*Rf^2/cos(2*pi/45))/pc;
以上是我自己做的程序,所得结果非常的不合情理,请高人帮我修改一下!
先谢谢了!
[ 本帖最后由 eight 于 2007-10-12 16:28 编辑 ] 调整初始值看看 不知道楼主是程序是解决什么问题的，最好把问题完整的贴上来才能帮你，只看到程序大家也不知道那里有问题 原帖由 li343 于 2007-9-17 07:46 发表 http://www.chinavib.com/forum/images/common/back.gif
这是一个MATLAB的M文件!
在保证Pf最小的情况下,用matlab优化Rf,Hpf,Hcf的值?
%***************************************************************
% This function is used to caculate the measurement ...
请具体说明结果如何不合理
页:
[1]