|
Revision 225, 0.6 kB
(checked in by smidl, 15 years ago)
|
|
prestavba mpf_test do inv-gamma
|
-
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 | for i=1:4 |
|---|
| 11 | subplot(4,1,i) |
|---|
| 12 | plot(xth(:,i)') |
|---|
| 13 | hold on; |
|---|
| 14 | plot(xthE(:,i)','r') |
|---|
| 15 | plot(xthM(:,i+4)','g') |
|---|
| 16 | % plot(xthV(i,:)','m'); |
|---|
| 17 | end |
|---|
| 18 | |
|---|
| 19 | figure(3) |
|---|
| 20 | for i=1:4 |
|---|
| 21 | subplot(4,1,i) |
|---|
| 22 | plot(xthM(:,i)') |
|---|
| 23 | if i<4 |
|---|
| 24 | %set(gca,'XTick',[]); |
|---|
| 25 | else |
|---|
| 26 | xlabel('sample [t]'); |
|---|
| 27 | end |
|---|
| 28 | ylabel(['Q(' num2str(i) ',' num2str(i) ')']); |
|---|
| 29 | hold on |
|---|
| 30 | plot(Qtr(:,i)','--') |
|---|
| 31 | end |
|---|
| 32 | |
|---|
| 33 | |
|---|
| 34 | % figure(3) |
|---|
| 35 | % ndat = size(ll,2); |
|---|
| 36 | % hold off |
|---|
| 37 | % plot(llgrid,ll(:,floor(ndat/3)),':'); |
|---|
| 38 | % hold on |
|---|
| 39 | % plot(llgrid,ll(:,floor(2*ndat/3)),'--'); |
|---|
| 40 | % plot(llgrid,ll(:,ndat)); |
|---|
