-
Notifications
You must be signed in to change notification settings - Fork 0
/
plotUS.m
41 lines (39 loc) · 1.04 KB
/
plotUS.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
function alpha=plotUS(thisC)
c=1;
for a=1:length(thisC)
subplot(10,5,c)
region=thisC(a).Region;
fT=thisC(a).FitTotal;
fTGOF=thisC(a).FitTotalGOF;
fTD=thisC(a).DataFit;
fTT0=thisC(a).DataT0;
eTGOF=thisC(a).ExpFitTotalGOF;
if (fTGOF.sse<eTGOF.sse)
figure(5)
[fTGOF.sse eTGOF.sse]
plot(log10(fTT0(:,1)),log10(fTT0(:,2)),'k',log10(fTD(:,1)),fT(log10(fTD(:,1))),'r');
hold on
%title(region{1})
text(log10(1),log10(5000),region{1});
drawnow
%input('r')
alpha(c).V=fT.a;
alpha(c).R=region;
c=c+1;
else
figure(2)
[fTGOF.sse eTGOF.sse]
plot(log10(fTT0(:,1)),log10(fTT0(:,2)),'k',log10(fTD(:,1)),fT(log10(fTD(:,1))),'r');
hold on
%title(region{1})
text(log10(1),log10(5000),region{1});
drawnow
%input('r')
%alpha(c).V=fT.a;
%alpha(c).R=region;
%c=c+1;
end
box off
ylim([1 5])
end
alpha=struct2table(alpha);