why do i got this error?

Unable to perform assignment because the size of the left side is 153-by-1 and the size of the right side is 51-by-1-by-3.

d =calc(Vp,Vs,Rh,wavelet,ang);

when i want to perfurm for loop i hane got this error knowing that the calc is function give output 3darray

[KSSV]

You are trying to save wrong number of elements then the array is initialized for.

Show us the full code.

[RADWAN A F ZEYADI]

for i=1:71
    [Vp, Vs, Rho]=applico_RPM_II([Swmean(:,i), shmean(:,i), Swmean(:)]);
end
d =calc(Vp,Vs,Rh,wavelet,ang)

knowing that smean phimean shmean has 51*1

[Image Analyst]

Give us the function line (or the whole function) where you define calc(). Are you sure it wants 5 inputs and not 1 to 4?

And what does this say

which -all calc

[RADWAN A F ZEYADI]

function d=calc(Vp,Vs,Rho,wavelet,ang)

[M,N,Q]=size(Vp);
d_pre=single(zeros(M,N,length(ang),Q));

for i=1:Q
    for j=1:N
        pippo=CMP_zoepprtiz([Vp(:,j,i) Vs(:,j,i) Rho(:,j,i)],wavelet,ang);
        pippo=[pippo(1,:);pippo];        
        d_pre(:,j,1:length(ang),i)=single(pippo);
    end
    
end

end

[Image Analyst]

Why are you overwriting Vp, Vs, and Rho. each time? You're not indexing them to save the separate values so at the end of the loop the're all same scalar value, so you might as well have done the loop only once.

Attach Shmean and swmean in a .mat file so we can run your code.

[RADWAN A F ZEYADI]

ok but the file is large size zipfile because i have used 3 functions wuth main script

[Image Analyst]

I find it hard to believe all your scripts zip up to over 5 MB in size (the attachmente limit). Go ahead and attach them if you want more help. By the way, did you try my suggestion in my Answer below? Or did you even see it (scroll down)?

[RADWAN A F ZEYADI]

dear sir ,

i mean my code and the theories behind not only the scripts ,itried but still have the same errors

[Image Analyst]

@RADWAN A F ZEYADI, I assume you have it figure out now since nothing is attached.

[Image Analyst]

Which m-file is the main program? Do any of the m-files require anything from the wavelet or symbolic toolbox? If so I won't be able to run it.

[RADWAN A F ZEYADI]

untitled777 is the main script and the other the functions ,no noneed to toolbox

[Image Analyst]

>> united777

Error using load

Unable to find file or directory 'model_petro'.

Error in united777 (line 7)

load model_petro

Please attach any files needed at runtime, such as "model_petro.mat".

Also I searched for the line of code you first mentioned and "calc(Vp" does not appear in any of your files. Where is it?

[RADWAN A F ZEYADI]

Cannot attach this file because:

• You are limited to 10 daily uploads. If you need to upload additional files, delete one or more files now or wait 24 hours to upload more files.

about clc is the same calclo_dati_presack

[Image Analyst]

OK, I've downloaded them all, so edit your post and delete one attachment like it says, and upload the mat file.

[RADWAN A F ZEYADI]

i have upload it,

Answer 1

[Image Analyst]

You can try

t = calc(Vp,Vs,Rh,wavelet,ang) % t should by 51x1x3.
% Squeeze to get a 51x3
t = squeeze(t); % A 51c3
whos t
d = t(:, 1); % Take column 1 of t and put it into d

Answer 2

[Image Analyst]

I don't know exactly what you're doing but look at this:

function [updatemodel, k_gain] = ESMDA_inversion(w, data, d_pred, alfa, Cd)
[nd, ne] = size(d_pred);
% data perturbation
term1 = repmat(data,1,ne);
term2 = sqrt(alfa*Cd);
term3 = randn(nd, ne);
term4 = term2 * term3;
whos data
whos term1
whos term2
whos term3
whos term4
dpert=repmat(data,1,ne) + sqrt(alfa*Cd)*randn(nd, ne);
% mean models
mm = mean(w, 2);
md = mean(d_pred, 2);
% covariance matrices
cmd = 1/(ne-1) * (w - mm) * (d_pred - md)';
cdd = 1/(ne-1) *(d_pred - md) * (d_pred - md)';
% Kalman Gain
k_gain = cmd  * inv(cdd + alfa*Cd);
% Updated models
updatemodel = w + k_gain*(dpert - d_pred);
end

If you run that you'll see

Name Size Bytes Class Attributes

data 51x3 612 single

Name Size Bytes Class Attributes

term1 51x1500 306000 single

Name Size Bytes Class Attributes

term2 153x153 93636 single

Name Size Bytes Class Attributes

term3 153x500 612000 double

Name Size Bytes Class Attributes

term4 153x500 306000 single

and you're trying to add term4 to term1. But you cannot add a 153 row matrix to a 51 row matrix. To add matrices they must have the same number of rows and columns. Examine the logic of your program and try to figure out what you really want to do.

[RADWAN A F ZEYADI]

thank you but i have convert the data to 153*1 and i tried to apply for loop then i have to obtain 51*71 for3 d

this is the main idea in eventually i have to apply dcosine transform to compress the final images

[RADWAN A F ZEYADI]

whos

Name Size Bytes Class Attributes

Cd 153x153 93636 single

Cm 153x153 187272 double

DX 1x1 8 double

DY 1x1 8 double

M 1x1 8 double

N 1x1 8 double

NI_TAU_Y 51x51 20808 double

Phi 51x500 102000 single

Phi_sim 51x500 204000 double

Phi_sim_logit 51x500 204000 double

Phi_true 51x71 28968 double

Phi_true_logit 51x71 28968 double

Phimean 51x1 204 single

Phiup 51x500 102000 single

Rho_true 51x71 28968 double

Rhosim 51x500 204000 double

Rhoup 51x500 204000 double

Sh 51x500 102000 single

Sh_sim 51x500 204000 double

Sh_sim_logit 51x500 204000 double

Sh_true 51x71 28968 double

Sh_true_logit 51x71 28968 double

Shup 51x500 102000 single

Sw 51x500 102000 single

Sw_sim 51x500 204000 double

Sw_sim_logit 51x500 204000 double

Sw_true 51x71 28968 double

Sw_true_logit 51x71 28968 double

Swmean 51x1 204 single

Swup 51x500 102000 single

Vp_true 51x71 28968 double

Vpsim 51x500 408000 double complex

Vpup 51x500 408000 double complex

Vs_true 51x71 28968 double

Vssim 51x500 408000 double complex

Vsup 51x500 408000 double complex

a 1x1 8 double

alfa 1x1 8 double

ang 1x3 24 double

b 1x1 8 double

c 1x1 8 double

cdd 153x153 93636 single

cmd 153x153 93636 single

cov_logit 3x3 72 double

d 1x1 8 double

d_obs 51x71x3 43452 single

d_obs_noise 51x71x3 43452 single

d_pred 153x500 306000 single

data 153x1 612 single

dpert 153x500 306000 single

dt 1x1 8 double

far 51x1 204 single

farp 51x500 102000 single

iter 1x1 8 double

j 1x1 8 double

k_gain 153x153 93636 single

m_sim 153x500 612000 double

md 153x1 612 single

mean_logit 153x1 1224 double

mid 51x1 204 single

midp 51x500 102000 single

mm 153x1 1224 double

mpost 153x1 612 single

nang 1x1 8 double

nd 1x1 8 double

ne 1x1 8 double

near 51x1 204 single

nearp 51x500 102000 single

ni_tau_y 1x101 808 double

shmean 51x1 204 single

std_noise 1x1 4 single

tau_y 1x101 808 double

updatemodel 153x500 306000 single

w 153x500 612000 double

wavelet 1x61 488 double

well_log_phi 51x3 1224 double

well_log_sh 51x3 1224 double

well_log_sw 51x3 1224 double

wup 153x500 306000 single

xax 1x201 1608 double

yax 1x51 408 double

[Image Analyst]

OK but my answer stands. I put a breakpoint on the line where it tried to add two differently shaped matrices. You need to fix that.