imnoise3.m

来自「matlab中产生常用的噪声」· M 代码 · 共 20 行

function [r,R,s]=imnoise3(c,a,b)
c=double(c);
[M,N]=size(c);
if nargin==3
    a(1:M)=1.0;
    b(1:M,1:2)=0;
elseif nargin==4
    b(1:M,1:2)=0;
end
R=zeros(M,N);
for j=1:M
    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,1)*b(j,1)/N);
end
s=abs(R);
r=real(ifft2(ifftshift(R)));