请帮忙修改一个优化程序

li343

在保证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-9-16 23:19 编辑 ]

xjzuo

先把你的问题及参数用word上传一下.---否则很难猜测你在做什么,更不用说提出更好的办法了.
另: 你的代码存在明显问题.

li343

首先感谢xjzuo版主的提示，现把已把问题说明及变量做成word上传！

问题说明.doc (83 KB, 下载次数: 12)

2007-9-17 12:40 上传

点击文件名下载附件

%%%%%%%%%%%变量表%%%%%%%%%%%%%%%%%
%ts:地面黑体温度；                          tm:接近烟筒入口处的最大温度；
%ti:集热棚和地面之间气流的平均温度；ta:集热棚周边入口处的环境温度；
%rm:为tm温度下空气的密度；             ra:为ta温度下空气的密度；
%k: 空气导热系数；                          g: 重力加速度；
%cp:定压比热容；                             rf:热扩散率；
%v: 空气的运动粘度；                       r: 烟筒的半径；
%yt:电机转子效率；                         pout:发电机的输出功率；
%pr:集热棚材料价格；                      pc:烟筒材料价格；
%Rf:集热棚半径；                            Hpf:集热棚距地面的高度；
%Hcf:烟筒的高度；                           Pf:太阳能发电机的总成本；
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%***************************************************************
% 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;

请版主帮我再看看如何修改．

[ 本帖最后由 li343 于 2007-9-17 12:40 编辑 ]

xjzuo

看了一下你的说明文档------这应该是一个“无约束非线性多变量优化问题”。
可以直接利用fminsearch或者fminunc函数求解，你可以试试。
（看看help，或者搜索论坛均可找到例子）

另：现在再返回来看你的代码，似乎离求解问题很远。