Changeset 761 for applications/bdmtoolbox/tutorial
- Timestamp:
- 01/06/10 22:30:00 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/bdmtoolbox/tutorial/mpdm/dist_estim_example.m
r759 r761 16 16 fz.rvc = RVtimes([u,u], [0, -1]); 17 17 fz.A = [1.8, -0.9]; 18 fz.const = -0.2;18 fz.const = 0.2; 19 19 fz.R = 1e-2; 20 20 … … 35 35 A1.rgr = RVtimes([z,z],[0,-1]) ; % correct structure is {y,y} 36 36 A1.options ='logfull,logll'; 37 A1.rv_param = RV({'a','b','c','r'}); 37 38 A1.frg = 0.95; 38 39 A1.constant = 1; … … 40 41 A2=A1; 41 42 A2.rv = z; 43 A2.rv_param = RV({'a','b','c','r'}); 42 44 A2.rgr = RVtimes([u,u],[0,-1]) ; % correct structure is {y,y} 43 45 … … 52 54 Merger.stop_niter=10; 53 55 Merger.effss_coef=0.9; 56 %Merger.dbg_file = 'm.it'; 57 58 M2.class = 'merger_base'; 59 M2.method='lognormal'; 60 M2.beta=1.2; 54 61 55 62 support.class='discrete_support'; 56 support.epdf= struct('class','enorm<ldmat>','mu',[ 1,1,1,1,0.1],'R',eye(5));57 support.npoints=[ 100];63 support.epdf= struct('class','enorm<ldmat>','mu',[2,-1,0,0.1],'R',eye(4)); 64 support.npoints=[300]; 58 65 66 support2.class='rectangular_support'; 67 support2.ranges= {[0,5],[-2,2],[-2,2],[0.001,3]}; 68 support2.gridsizes = [10,10,10,10]; 59 69 60 70 … … 69 79 A1_mean = zeros(4,ndat); 70 80 A2_mean = zeros(4,ndat); 71 MG_mean = zeros( 5,ndat);81 MG_mean = zeros(4,ndat); 72 82 73 83 for t=1:ndat 74 84 f1=Set.Est0_apost{t}; 75 f1.rv = RV({'a','b','c','r'});76 85 f2=Set.Est1_apost{t}; 77 f2.rv = RV({'a','b','d','r'}); 78 86 if t>5 87 % support.epdf = f1; 88 end 89 79 90 [res,ftilde] = merger({f1,f2},support,Merger); 80 91 81 92 A1_mean(:,t) = epdf_mean(f1); 82 93 A2_mean(:,t) = epdf_mean(f2); … … 100 111 101 112 figure(2) 102 subplot(2,2,1); 103 plot(A1_mean'); 104 subplot(2,2,2); 105 plot(A2_mean'); 106 subplot(2,2,3); 113 plot(A1_mean',':'); 114 hold on; 115 plot(A2_mean','--'); 107 116 plot(MG_mean');
