function [r,R,S]=imnoise3(M,N,C,A,B) %產(chǎn)生一個(gè)大小為M*N的正弦噪聲模型r,R為傅里葉變換，S為正弦噪聲模型的傅里葉的頻譜 %C為沖擊位置的坐標(biāo)，A是1*k維向量包含振幅的沖擊對(duì)，B是由k*2矩陣組成的沖擊對(duì)%處理輸入?yún)?shù) [K,n]=size(C);%矩陣C的行數(shù)返回給K，矩陣C的列數(shù)返回給n if nargin==3A(1:K)=1.0;B(1:K,1:2)=0; elseif nargin==4B(1:K,1:2)=0; end R=zeros(M,N);%構(gòu)造R for j=1:K%j從1到K取遍，步長(zhǎng)為1u1=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);%復(fù)共軛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);%計(jì)算光譜 r=real(ifft2(ifftshift(R)));%產(chǎn)生空間正弦圖形

總結(jié)

以上是生活随笔為你收集整理的imnoise3.m的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。