WHEN I AM RUNNING THIS CODE I AM GETTING FOLLOWING ERROR

Question

sabitri 2022-7-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1764210-when-i-am-running-this-code-i-am-getting-following-error

编辑： sabitri 2022-7-20

采纳的回答： VBBV

function [p,s] = perdecomp(u)

[ny,nx] = size(u);

u = double(u);

X = 1:nx;

Y = 1:ny;

v = zeros(ny,nx);

v(1,X) = u(1,X)-u(ny,X);

v(ny,X) = -v(1,X);

v(Y,1 ) = v(Y,1 )+u(Y,1)-u(Y,nx);

v(Y,nx) = v(Y,nx)-u(Y,1)+u(Y,nx);

fx = repmat(cos(2.*pi*(X -1)/nx),ny,1);

fy = repmat(cos(2.*pi*(Y'-1)/ny),1,nx);

fx(1,1)=0.; % avoid division by 0 in the line below

s = real(ifft2(fft2(v)*0.5./(2.-fx-fy)));

p = u-s;

Arrays have incompatible sizes for this operation.

Error in perdecomp (line 17)

p = u-s;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

VBBV 2022-7-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1764210-when-i-am-running-this-code-i-am-getting-following-error#answer_1011545

在 MATLAB Online 中打开

u = rand(1,5) % assume a matrix values 
u = 1×5
    0.1621    0.3146    0.5520    0.9985    0.1474
[p,s]= perdecomp(u)  % call the function giving input matrix u 
p = 1×5
    0.1562    0.3116    0.5520    1.0015    0.1532
s = 1×5
    0.0059    0.0029    0.0000   -0.0029   -0.0059
function [p,s] = perdecomp(u) 
 [ny,nx] = size(u);  
 u = double(u);  
X = 1:nx;  
Y = 1:ny; 
v = zeros(ny,nx); 
v(1,X)  = u(1,X)-u(ny,X);
 v(ny,X) = -v(1,X); 
v(Y,1 ) = v(Y,1 )+u(Y,1)-u(Y,nx); 
v(Y,nx) = v(Y,nx)-u(Y,1)+u(Y,nx); 
fx = repmat(cos(2.*pi*(X -1)/nx),ny,1); 
fy = repmat(cos(2.*pi*(Y'-1)/ny),1,nx); 
fx(1,1)=0.;   % avoid division by 0 in the line below 
s = real(ifft2(fft2(v)*0.5./(2.-fx-fy))); 
p = u-s; 
end

Give input vector when calling function perdecomp

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

sabitri 2022-7-20

编辑：sabitri 2022-7-20

Sir, i have to perform this decomsition on an jpg image. when i am writing u=imread('u.jpg') i am getting a long matrix..How can i solve this? here u is an jpg file saved in the folder. @VBBV

请先登录，再进行评论。