imnoise3.m

来自「Matlab數位影音處理 影像方面的處理」· M 代码 · 共 18 行

function [r,R,S]=imnoise3(M,N,C,A,B)
[k,n]=size(C);
if nargin==3
    A(1:k)=1.0;
    B(1:k,1:2)=0;
elseif nargin ==4
    B(1:k,1:2)=0
end
R=zeros (M,N);
for j=1:k
    u1=M/2+1+C(j,1);v1=N/2+1+C(j,2);
    R(u1,v1)=i*(A(j)/2)*exp(i*2*pi*C(j,1)*B(j,1)/M);
    u2=M/2+1-C(j,1);v2=N/2+1-C(j,2);
    R(u2,v2)=-i*(A(j)/2)*exp(i*2*pi*C(j,2)*B(j,2)/N);
end
S=abs(R);
r=real(ifft2(ifftshift(R)));