function [r,R,S]=imnoise3(M,N,C,A,B) %產生一個大小為M*N的正弦噪聲模型r,R為傅里葉變換，S為正弦噪聲模型的傅里葉的頻譜 %C為沖擊位置的坐標，A是1*k維向量包含振幅的沖擊對，B是由k*2矩陣組成的沖擊對%處理輸入參數 [K,n]=size(C);%矩陣C的行數返回給K，矩陣C的列數返回給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);%構造R for j=1:K%j從1到K取遍，步長為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);%復共軛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)));%產生空間正弦圖形

總結

以上是生活随笔為你收集整理的imnoise3.m的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。