Ci siamo quasi...
ora ho 2 caratteristiche della potenza...verde e con asse verde che si intromette...
non so come si postano le immagini ma se fai partire tutto il codice lo puoi notare tu stesso:
function corpot
axis([0 2.5 0 1.4]);hold on
Isc=1.27; %corrente di corto circuito, [mA]
Io=6.938299e-20; %corrente di saturazione del diodo, [mA]
q=1.602e-19; %carica dell'elettrone, [C]
k=1.381e-23; %costante di Boltzmann, [J/K]
T=298; %temperatura, [K]
V=0:.001:2.28; %tensione circuito aperto [V]
n=2; %fattore di idealità
for i=1:length(V)
I(i)=Isc-Io*(exp(V(i)*q/(n*k*T))-1);
P(i)=V(i)*I(i);
end
[AX,H1,H2] = plotyy(V,I,V,P,'plot');
set(get(AX(1),'Ylabel'),'String','Current Density (mA/cm^2)','Color','k')
set(get(AX(2),'Ylabel'),'String','Power Density (mW/c^2)','Color','k')
xlabel('Voltage (V)')
title('Caratteristica J-V e P-V in STC')
%set(H1,'LineStyle','-','Color','k')
%set(H2,'LineStyle','-','Color','r')
%legend({'Current Density','Potenz Density'})
ax1=gca;
gr1=line(V,I, 'Color', 'k');
ax2 = axes('Position',get(ax1,'Position'), 'YAxisLocation','right', 'Color','none', 'XColor','k','YColor','k');
gr2 = line(V,P,'Color','r','Parent',ax2);
set(ax2, 'YLim',[0 2.8], 'YTick',0:0.4:2.8);
non so come ringraziarti per l'aiuto