Design og modellering af metanolanlæg til VEnzin-visionen Bilag

18-03-07 23:10 D:\DTU\Eksamensprojekt\bilag\elektrolysedrift.m 4 of 12
k))*DogV*levetid;
151 DogV_samlet2(k)=DogV_samlet(iterationer,k);
152 sparet(iterationer,k)=sparet_el(iterationer,k)-(ekstraanlagsudgift(k)
+Lagerprisen(iterationer,k))*(1+DogV*levetid);
153 sparet2(k)=sparet(iterationer,k);
154 andel_sparet(k)=sparet(iterationer,k)/ref_omkostninger;
155 lagerpris_andel(k)=Lagerprisen2(k)/ref_anlagsudgift;
156 tilbagebetalingstid_aar(k)=(ekstraanlagsudgift(k)+Lagerprisen2(k))/
(sparet_el2(k)/levetid-(ekstraanlagsudgift(k)+Lagerprisen2(k))*DogV);
157 forrentning(k)=1/tilbagebetalingstid_aar(k);
158
159 sparet_sum_nutidsvardi(k)=0;
160 stop_nutidsvardi=0;
161 for i=1:50
162 sparet_aarligt=(sparet_el2(k)/levetid-(ekstraanlagsudgift(k)
+Lagerprisen2(k))*DogV)/(1+kalkulationsrente)^i;
163 sparet_sum_nutidsvardi(k)=sparet_sum_nutidsvardi(k)+sparet_aarligt;
164 if sparet_sum_nutidsvardi(k) >= ekstraanlagsudgift(k)+Lagerprisen2(k);
165 andel_af_aar=(sparet_sum_nutidsvardi(k)-(ekstraanlagsudgift(k)
+Lagerprisen2(k)))/sparet_aarligt;
166 tilbagebetalingstid2_aar(k)=i-andel_af_aar;
167 stop_nutidsvardi=1;
168 break
169 end
170 i=i+1;
171 end
172 if stop_nutidsvardi == 0
173 tilbagebetalingstid2_aar(k)=99;
174 end
175 sparet_sum_nutidsvardi(k)=0;
176 for i=1:levetid
177 sparet_sum_nutidsvardi(k)=sparet_sum_nutidsvardi(k)+(sparet_el2(k)
/levetid-(ekstraanlagsudgift(k)+Lagerprisen2(k))*DogV)/(1+kalkulationsrente)^i;
178 end
179 sparet_nutidsvardi(k)=sparet_sum_nutidsvardi(k)-ekstraanlagsudgift(k)-
Lagerprisen2(k);
180 andel_sparet_nutidsvardi(k)=sparet_nutidsvardi(k)/ref_omkostninger;
181
182 % Lokaliserer "peak"- og "off-peak"-områder i lagerbeholdningen henover
alle
183 % timerne. Et "peak"-område er et område hvor der forekommer en stigning i
lagerbeholdningen fra minimumbeholdningen til maksimumbeholdningen.
184 % Omvendt er et "off-peak"-område et område hvor der forekommer et fald i
lagerbeholdningen fra maksimumbeholdningen til minimumbeholdningen.
185 [minlager,minnr]=min(lager(:,k));
186 [maxlager,maxnr]=max(lager(:,k));
187 i=1;
188 s=1;
189 f=1;
190 lagerhistorie=1;
191 while i maxlager-lagerreduktion2*3
193 if lagerhistorie == 0
194 stigning(s,2)=i;
195 s=s+1;
196 elseif lagerhistorie == 1