Как получить спектр сигнала?

I=imread('213.tif'); %так я хотел получить значения яркости
I2 = im2double(I);
r1 = smooth(I2,20);
R1=reshape(r1,1024,768); 
 
figure(1);
 
subplot(2,1,1)
plot (I2);
subplot(2,1,2)
plot(R1,'-');
 
signal= R1(1:1024,2); %исследуемый сигнал 
dsignal=diff(signal); %его производная 
 
figure(2);
 
subplot(2,1,1);
plot(signal);
subplot(2,1,2);
plot(dsignal);
 
figure(3);
fsignal=fft(signal);
subplot(2,1,1);
plot(fsignal);
fdsignal=fft(dsignal);
subplot(2,1,2);
plot(fdsignal);

I=imread('213.tif'); %так я хотел получить значения яркости
I2 = im2double(I);
r1 = smooth(I2,20);
R1=reshape(r1,1024,768); 
 
figure(1);
 
subplot(2,1,1)
plot (I2);
subplot(2,1,2)
plot(R1,'-');
 
signal= R1(100:300,2); %исследуемый сигнал 
signal1=R1(1:1024,2);
dsignal=diff(signal); %его производная 
 
figure(2);
 
subplot(2,1,1);
plot(signal);
subplot(2,1,2);
plot(dsignal);
 
figure(3);
t = 0:1024; % время
s = signal; % сигнал
dt = t(2)-t(1); % интервал дискретизации
Fs = 1/dt; % частота дискретизации
N = length(s); % длина сигнала
f = (0:N-1)/N*Fs; % вектор частот для БПФ
plot(f, abs(fft(s)))

I=imread('213.tif'); %так я хотел получить значения яркости
I2 = im2double(I);
r1 = smooth(I2,20);
R1=reshape(r1,1024,768); 
 
figure(1);
 
subplot(2,1,1)
plot (I2);
subplot(2,1,2)
plot(R1);
 
signal= R1(1:1024,2); %исследуемый сигнал 
 
dsignal=diff(signal); %его производная 
 
figure(2);
 
subplot(2,1,1);
plot(signal);
subplot(2,1,2);
plot(dsignal);
 
figure(3);
subplot(2,1,1);
t = 1:1024; % время
s = signal; % сигнал
dt = t(2)-t(1); % интервал дискретизации
Fs = 1024;%1/dt; % частота дискретизации
N = length(s); % длина сигнала
f = (0:N-1)/N*Fs; % вектор частот для БПФ
plot(f, abs(fft(s)))
Y=fft(s);
 
subplot(2,1,2);
t1 = 1:1022; % время
s1 = dsignal; % сигнал
dt1 = 1; % интервал дискретизации
Fs1 = 1024;%1/dt; % частота дискретизации
N1 = length(s1); % длина сигнала
f1 = (0:N1-1)/N1*Fs1; % вектор частот для БПФ
plot(f1, abs(fft(s1)))
Y1=fft(s);
 
figure(4);
subplot(2,1,1);
Fp=Y/Y1;
pX = ifft(Fp);
plot(pX)
 
figure(5);
 
subplot(1,1,1)
R2=reshape(R1,768,1024); 
ss = fft2(R2);
imagesc(abs(fftshift(ss)))
%subplot(2,1,2)
%plot(R1);