h = 6.626e-34; c = 3e8; k = 1.380e-23; T = 2800; L = 1; D = 0.02e-3; l = linspace(300e-9,900e-9,1000); P = 8*pi*h*c^2*L*D*pi ./ l.^5 ./ ( exp( h*c/k./l./T ) - 1 ); y1 = 0.5; dist = 0.3; figure('Position',[50,50,400,350],'Color','w'); plot(l*1e9,P*1e-9,'LineWidth',2); hold on; plot(l*1e9,0.3*P*1e-9,'LineWidth',2,'Color','r'); hold on; xlabel('Wellenlänge \lambda (nm)'); ylabel({'Strahlungsleistung \Phi_\lambda(\lambda) (W/nm)'}) set(gca,'Position',[0.15 0.15 0.71 0.81]); ylim([0,y1]) h_leg = legend({'\Phi_\lambda(\lambda), I_\lambda(\lambda) Planckscher Körper'},'\Phi_\lambda(\lambda), I_\lambda(\lambda) Emissionskoeffizient 0.3'); set(h_leg,'Position',[0.2 0.5 0.58 0.15]); legend boxoff; ax2 = axes('Position',get(gca,'Position')... ,'XAxisLocation','bottom'... ,'YAxisLocation','right'... ,'Color','none'... % Hintergrund des Graphen ist farblos, vorheriger Plot scheint durch ); set(ax2,'XTick',[]) set(ax2,'Ylim',[0,y1*100/pi/dist^2]) ylabel('Bestrahlungsstärke I_\lambda(\lambda) 30 cm Abstand (µW/cm^2/nm)')