|
Revision 229, 0.8 kB
(checked in by smidl, 17 years ago)
|
|
epdf has a new function: variance()
|
-
Property svn:eol-style set to
native
-
Property svn:keywords set to
Rev Author Date
|
| Line | |
|---|
| 1 | close all |
|---|
| 2 | itload('../../mpf_test.it') |
|---|
| 3 | figure(1); |
|---|
| 4 | for i =1:4 |
|---|
| 5 | subplot(4,1,i) |
|---|
| 6 | plot(xth(:,i)') |
|---|
| 7 | end |
|---|
| 8 | |
|---|
| 9 | figure(2) |
|---|
| 10 | StdX = 2*sqrt(VarE); |
|---|
| 11 | for i=1:4 |
|---|
| 12 | subplot(4,1,i) |
|---|
| 13 | plot(xth(:,i)') |
|---|
| 14 | hold on; |
|---|
| 15 | plot(xthE(:,i)','r') |
|---|
| 16 | plot(xthE(:,i) - StdX(:,i),'r:'); |
|---|
| 17 | plot(xthE(:,i) + StdX(:,i),'r:'); |
|---|
| 18 | plot(xthM(:,i+4)','g') |
|---|
| 19 | % plot(xthV(i,:)','m'); |
|---|
| 20 | end |
|---|
| 21 | |
|---|
| 22 | figure(3) |
|---|
| 23 | StdM = 2*sqrt(VarM); |
|---|
| 24 | for i=1:4 |
|---|
| 25 | subplot(4,1,i) |
|---|
| 26 | plot(xthM(:,i)') |
|---|
| 27 | hold on |
|---|
| 28 | plot(xthM(:,i)-StdM(:,i),':') |
|---|
| 29 | plot(xthM(:,i)+StdM(:,i),':') |
|---|
| 30 | if i<4 |
|---|
| 31 | %set(gca,'XTick',[]); |
|---|
| 32 | else |
|---|
| 33 | xlabel('sample [t]'); |
|---|
| 34 | end |
|---|
| 35 | ylabel(['Q(' num2str(i) ',' num2str(i) ')']); |
|---|
| 36 | plot(Qtr(:,i)','--') |
|---|
| 37 | end |
|---|
| 38 | |
|---|
| 39 | |
|---|
| 40 | % figure(3) |
|---|
| 41 | % ndat = size(ll,2); |
|---|
| 42 | % hold off |
|---|
| 43 | % plot(llgrid,ll(:,floor(ndat/3)),':'); |
|---|
| 44 | % hold on |
|---|
| 45 | % plot(llgrid,ll(:,floor(2*ndat/3)),'--'); |
|---|
| 46 | % plot(llgrid,ll(:,ndat)); |
|---|
