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

18-03-07 23:10 D:\DTU\Eksamensprojekt\bilag\elektrolysedrift.m 9 of 12
401 time_doegn=time_doegn+1;
402 if time_doegn == 25
403 time_doegn = 1;
404 end
405 i=i+1;
406 end
407 lagerfejl=min(lager(:,k));
408 iii=iii+1;
409 end
410 timer_i_drift_=timer_i_drift(end,k)
411 Driftstimer
412 iterationer
413 sparet_=sparet2(k)
414 lagerfejl_start_stop=lager(end,k)-lager(1,k)
415
416 iterationer=iterationer-1;
417 iterationer_list(k)=iterationer;
418 andel_list(k)=andel;
419 Driftstimer_list(k)=Driftstimer;
420 Driftstimer=Driftstimer-timestep;
421 k=k+1;
422 end
423 [maxxx,kk]=max(sparet2)
424
425 figure(3);
426 bar(uge_drift(:,kk))
427 axis([0 53*antal_aar 0 24*7])
428 %title('Driftsfordeling på ugebasis')
429 ylabel('Driftstimer per uge [h]' )
430 xlabel('Uge-nr. [-]' )
431
432 figure(4);
433 plot(lager(:,kk))
434 axis tight
435 %title('Variationen i lagerbeholdningen henover året')
436 ylabel('Brintlagerbeholdning [MWh]' )
437 xlabel('Time på året [-]' )
438
439 figure(19);
440 plot(driftstatus_doegn_foer(:,kk),'k' )
441 hold on
442 plot(driftstatus_doegn_efter(:,kk),'b')
443 ylabel('Antal timer [-]' )
444 hold off
445 legend('før' ,'efter' ,1);
446
447 figure(7);
448 plot((ref_eludgift-sparet_el(1:iterationer_list(kk),kk))
/ref_omkostninger_aar*100,'-k' )
449 hold on
450 plot(Lagerprisen(1:iterationer_list(kk),kk)/ref_omkostninger_aar*100, '--r')
451 plot(anlagsudgift2(1:iterationer_list(kk),kk)/ref_omkostninger_aar*100,':g')
452 plot(DogV_samlet(1:iterationer_list(kk),kk)/ref_omkostninger_aar*100, ':b')
453 plot(sparet(1:iterationer_list(kk),kk)/ref_omkostninger_aar*100,'-.c' )
454 hold off
455 legend('El' ,'Gaslagre','Elektrolyseanlæg' ,'D&V','Sparede omkostninger' ,1);