how to make a subplot.i did it but it has little problem.can anyone please guide me.here are codes,i combined them but didn't give the correct subplot.
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
waqas muhammad
am 30 Jan. 2018
Kommentiert: waqas muhammad
am 30 Jan. 2018
here is the code.
m0=0.938;R=1.2*196.97^(1/3);
t0=0.130;beta=0.375;n0=2;
x2=[0.644417475728155;0.753640776699029;0.851941747572815;0.950242718446602;1.04854368932039;1.14684466019417;1.25606796116505;1.35436893203883;1.45266990291262;1.55097087378641;1.66019417475728;1.75849514563107;1.85679611650485;1.95509708737864;2.09708737864078;2.30461165048544;2.50121359223301;2.70873786407767;2.90533980582524;3.25485436893204;3.75728155339806;4.25970873786408]
y2=[2.02167630196276;1.92179861093512;1.73660235151668;1.49172630869016;1.21807538896743;1.10069420009424;0.898776476692455;0.733899710733955;0.630413465400298;0.442180172338893;0.361064078764099;0.280262865396340;0.217543860882062;0.160518229479539;0.0967134573521131;0.0612989735304463;0.0317252106715399;0.0181703510369072;0.0104069177104193;0.00418075489660610;0.00130367268973768;0.000550940754545585]
y4pos=[2.78255940220713;2.78255940220713;2.38658978685858;2.04696827180752;1.66810053720006;1.58489319246111;1.22712523985119;1;0.857695898590894;0.630957344480193;0.488527357151939;0.359381366380463;0.278255940220713;0.204696827180752;0.129154966501488;0.0814912746902075;0.0419007910578667;0.0238658978685858;0.0135935639087853;0.00569581081073769;0.00175567629127500;0.000698947320727348]
y4neg=[2.15443469003188;2.04696827180752;1.84784979742229;1.58489319246111;1.29154966501488;1.16591440117983;1;0.774263682681127;0.664082785063485;0.464158883361278;0.359381366380463;0.292864456462524;0.226754312587080;0.166810053720006;0.105250028527773;0.0630957344480193;0.0324422607917163;0.0194486243893736;0.0110775685050971;0.00441005945417674;0.00129154966501488;0.000541169526546464]
i=0;
for pt=0.01:0.01:4.5
i=i+1;
mt=sqrt(pt^2+m0^2);
j=0;
for r=0.01:0.01:R
j=j+1;
rho=atanh(0.5*beta*(n0+2)*(r/R)^n0);
w=pt*sinh(rho)/t0;
z=mt*cosh(rho)/t0;
Q(j)=pt*mt*r*besseli(0,w)*besselk(1,z)*0.01;
end
y10(i)=sum(Q);
end
pt=[0.01:0.01:4.5];
y11=y10./sum(y10);
y12=335*y11./pt;
m0=0.938;R=1.2*196.97^(1/3);
t0=0.130;beta=0.375;n0=2;
x4=[0.644417475728155;0.753640776699029;0.851941747572815;0.950242718446602;1.04854368932039;1.14684466019417;1.25606796116505;1.35436893203883;1.45266990291262;1.55097087378641;1.66019417475728;1.75849514563107;1.85679611650485;1.95509708737864;2.09708737864078;2.30461165048544;2.50121359223301;2.70873786407767;2.90533980582524;3.25485436893204;3.75728155339806;4.25970873786408]
y4=[2.02167630196276;1.92179861093512;1.73660235151668;1.49172630869016;1.21807538896743;1.10069420009424;0.898776476692455;0.733899710733955;0.630413465400298;0.442180172338893;0.361064078764099;0.280262865396340;0.217543860882062;0.160518229479539;0.0967134573521131;0.0612989735304463;0.0317252106715399;0.0181703510369072;0.0104069177104193;0.00418075489660610;0.00130367268973768;0.000550940754545585]
y4pos=[2.78255940220713;2.78255940220713;2.38658978685858;2.04696827180752;1.66810053720006;1.58489319246111;1.22712523985119;1;0.857695898590894;0.630957344480193;0.488527357151939;0.359381366380463;0.278255940220713;0.204696827180752;0.129154966501488;0.0814912746902075;0.0419007910578667;0.0238658978685858;0.0135935639087853;0.00569581081073769;0.00175567629127500;0.000698947320727348]
y4neg=[2.15443469003188;2.04696827180752;1.84784979742229;1.58489319246111;1.29154966501488;1.16591440117983;1;0.774263682681127;0.664082785063485;0.464158883361278;0.359381366380463;0.292864456462524;0.226754312587080;0.166810053720006;0.105250028527773;0.0630957344480193;0.0324422607917163;0.0194486243893736;0.0110775685050971;0.00441005945417674;0.00129154966501488;0.000541169526546464]
i=0;
for pt=0.01:0.01:4.5
i=i+1;
mt=sqrt(pt^2+m0^2);
j=0;
for r=0.01:0.01:R
j=j+1;
rho=atanh(0.5*beta*(n0+2)*(r/R)^n0);
w=pt*sinh(rho)/t0;
z=mt*cosh(rho)/t0;
Q(j)=pt*mt*r*besseli(0,w)*besselk(1,z)*0.01;
end
y30(i)=sum(Q);
end
pt=[0.01:0.01:4.5];
y11=y30./sum(y30);
y32=335*y11./pt;
x4i=x4
y4i=interp1(pt,y32,x4i)
kafa2=y4./y4i
subplot(2,1,1)
semilogy(pt,y12,'k')
%axis([0 4.5 1e-3 1e3])
xlabel('P_{T}')
ylabel('1/2\pip_{T} d^{2}N/dp_{T}dy [c^{2}/GeV^{2}]')
semilogy(x2,y2,'r*')
subplot(2,1,2)
scatter(x4,kafa2)
xlabel('p_{T}')
ylabel('data/fit')
c2=y4pos-y4neg
d1=c2./y4i% error=
errorbar(x4,kafa2,d1,'o')
%f=c./kafa2
kafang2=((y4-kafa2)/c2).^2
S2=sum(kafang2)
%now for finding chi^2/dof.. dof in no.of exp.points-no.of free parameters.
dof=22-3 %16 is the no. of experimental points,3 is the number of free parameters
K2=S2/dof
%err=0.3*ones(size(c))
%errorbar(f,err);
%err=0.3*ones(size(c))
%errorbar(f,err);
there are two graphs i have to plot in single window which are ...
the cobination gives me a graph like this..
now the first part of theis subplot doesn't show me the fitting line as in the i have in the first individual plot...
Akzeptierte Antwort
Abhishek Ballaney
am 30 Jan. 2018
https://in.mathworks.com/help/matlab/ref/subplot.html
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Subplots finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!