www.gusucode.com > 太阳能小屋的设计-数学建模【论文+程序】 > 2012高教社杯全国大学生数学建模竞赛太阳能小屋的设计【论文+程序】/文章+源码/B149源程序/M/zuijiajiaodu.m

    clear;clc
fushe=xlsread('cumcm.xls','sheet','E4:K8763');
dc=xlsread('cumcm.xls','sheet1','B1:J24');
sp_zs=fushe(:,1)-fushe(:,2);
n_zs=fushe(:,5)-0.5*fushe(:,2);
d_zs=fushe(:,4)-0.5*fushe(:,2);
x_zs=fushe(:,6)-0.5*fushe(:,2);
fdl=[];

%for i=1:91
    for j=1:181
    i=91;
        a=(i-1)*pi/180;b=(j-91)*pi/180;
        sa=sin(a);ca=cos(a);
        sb=sin(b);cb=cos(b);
        if sb<0
            fushe_ry=-d_zs*sa*sb+n_zs*sa*cb+sp_zs*ca+fushe(:,2)*(pi-a)/pi;
        else fushe_ry=x_zs*sa*sb+n_zs*sa*cb+sp_zs*ca+fushe(:,2)*(pi-a)/pi;
        end
        %计算月发电量
        S(1,j)=sum(fushe_ry(1:744));
        S(2,j)=sum(fushe_ry(745:1416));
        S(3,j)=sum(fushe_ry(1417:2160));
        S(4,j)=sum(fushe_ry(2161:2880));
        S(5,j)=sum(fushe_ry(2881:3624));
        S(6,j)=sum(fushe_ry(3625:4344));
        S(7,j)=sum(fushe_ry(4345:5088));
        S(8,j)=sum(fushe_ry(5089:5832));
        S(9,j)=sum(fushe_ry(5833:6552));
        S(10,j)=sum(fushe_ry(6553:7296));
        S(11,j)=sum(fushe_ry(7297:8016));
        S(12,j)=sum(fushe_ry(8017:8760));

    %end
end
x=1:181
plot(x,S)

%[x y]=find(S==max(max(S)))