请帮忙修改一个优化程序

li343

这是一个MATLAB的M文件!

在保证Pf最小的情况下,用matlab优化Rf,Hpf,Hcf的值?

%***************************************************************
% This function is used to caculate the measurement of solar
% collector.
%***************************************************************
function [Pf,Rf,Hpf,Hcf]=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=[1:15];
Hp=[0.1:0.02:2.5];
%**************************************************************
%           Grid Paremeters
%**************************************************************
x=R;
y=Hp;
[X,Y]=meshgrid(x,y);
P=roc.*X.^2+hoc.*Y.^(0.80)./X;
%***************************************************************
%          Get mininum value of P
%***************************************************************
a=min(P);
Pf=min(a);
[s,t]=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 编辑 ]

风花雪月

调整初始值看看

94117239

不知道楼主是程序是解决什么问题的，最好把问题完整的贴上来才能帮你，只看到程序大家也不知道那里有问题

eight

原帖由 li343 于 2007-9-17 07:46 发表

                               
登录/注册后可看大图

这是一个MATLAB的M文件!

在保证Pf最小的情况下,用matlab优化Rf,Hpf,Hcf的值?

%***************************************************************
% This function is used to caculate the measurement ...

请具体说明结果如何不合理