zhifangtu.m

x=2*sin(4*pi*t)+5*cos(8*pi*t) 对其进行快速傅立叶变换 加入高斯噪声 再进行频谱分析

g=imread (00.jpg);
I=rgb2gray(g);
J=I;
%计算灰度图象的直方图数据
L=256; %灰度级 
Ps = zeros(L,1); %出现频率%统计直方图结果数据
nk=zeros(L,1);%出现次数
[row,col]=size(I);
n=row*col; %总像素个数
for i = 1:row
for j = 1:col
num = double(I(i,j))+1; %获取像素点灰度级
nk(num) = nk(num)+1; %统计nk 
end
end
%计算直方图概率估计
for i=1:L
Ps(i)=nk(i)/n; 
end
figure;
subplot(3,1,1);imshow(J),title('原图');
subplot(3,1,2),plot(nk),title('直方图（nk）');
subplot(3,1,3),plot(Ps),title('直方图（Ps）');