2012年建模全国优秀作品程序学习文件.docVIP

程序编号 C1—1 文件名称 Ratio.m 说明 每平方光伏电池的性价比 % 每平方的性价比 clear,clc %% 读入数据 data=xlsread('cumcm2012B_附件3_三种类型的光伏电池(A单晶硅B多晶硅C非晶硅薄膜)组件设计参数和市场价格.xls'); pice=[14.9 12.5 4.8]; long=data(:,2);%长 wide=data(:,3);%宽 U=data(:,4);%电压 I=data(:,5);%电流 eta=data(:,6);%转换率 P=U.*I; S=long.*wide/1000; %% 每平米价格 for i=1:6 p1(i)=P(i)*pice(1)/S(i); end %A单晶硅 for i=7:13 p1(i)=P(i)*pice(2)/S(i); end %B多晶硅 for i=14:24 p1(i)=P(i)*pice(3)/S(i); end %C非晶硅薄膜 %% 每平方的性价比 ratio=eta./p1'; 程序编号 C1—2 文件名称 profit.m 说明 c类电池每块安装在四面墙上的35年利润 clc;clear;close all data=xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls'); data1=data(:,end-3:end);%东南西北的辐射数据 data2=data1; data2(find(data2<30))=0; he=sum(data2); mpower=he./1000; %每平米一年的发电量 power=mpower*10+mpower*15*0.9+mpower*10*0.8; %每平米35年的发电量 没有加入逆变效率 price=power*0.5; % 每平米的面积35年的经济效益 没有加入转换率 data3=xlsread('cumcm2012B_附件3_三种类型的光伏电池(A单晶硅B多晶硅C非晶硅薄膜)组件设计参数和市场价格.xls'); pice=4.8; long=data3(:,2);%长 wide=data3(:,3);%宽 U=data3(:,4);%电压 I=data3(:,5);%电流 eta=data3(:,6);%转换率 P=U.*I; S=long.*wide/1000000; price1=pice.*P;%每块电池的价格 lr=zeros(24,4); for i=1:24 lr(i,:)=price*S(i)*eta(i)-price1(i); %每块电池不考虑逆变器时35年的利润 end clr=lr(14:24,:) % c类电池每块安装在四面墙上的35年利润 程序编号 C1—3 文件名称 Account1.m 说明 第一问的相关计算程序 clc;clear;close all %% 数据的读入 data=xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls'); data1=data(:,3);%水平面总辐射强度 data2=data(:,4);%水平面散射辐射强度 data3=data1-data2;%水平面上直射强度 hpi=40.1*pi/180;%大同的纬度 %% 参数符号说明 %phi是当地纬度；beta是光伏阵列的倾角；delta为太阳赤纬角； %omegap为水平面日落时角；romegat为倾斜面日落时角。 %Rb为倾斜面上的直接辐射量与水平面上直接辐射量之比 % Rb=(cos(hpi-beta).*cos(delta).*sin(omegat)+pi/180*sin(hpi-beta)sin(delta))./(cos(hpi)*cos(delta)*sin(omegap)+pi/180*omegap*sin(hpi)*sin(delta)) % delta=23.5*sin((2*pi*(284+n))/365)*pi/180; % omegap=acos(-tan(hpi)*tan(delta)); % omegat=min(omegap,acos(-tan(hpi-beta)*tan(delta))); %% 南面屋顶 %选用36块B3多晶硅电池 用两个SN14逆变器 % B3的参数U=33.6; I=8.33; 价格12.5 尺寸1482*992 转换率15.98% % 逆变器的价格 price2=15300 逆变效率94% n=1:365; delta=23.5*sin((2*pi*(284+n))/365)*pi/180; omegat=zeros