2 семестр / Курсовая работа / forkvar / forkrvar32

clear all; clc; close all;
graphics_toolkit('fltk');
a=input('a=','s'); % Ввод номера файла
% или a=input('a='); % Ввод номера файла
% тогда нужна следующая строка a=int2str(a);
fd1 = fopen([a +'.raw'], 'r'); % Дескриптор файла
[X, cnt] = fread(fd1, 40001);% cnt ограничение Количества прочитанных байт
fclose(fd1);
clear fd1;
Signal=X';

%% Параметры
Fd=8000;% Частота дискретизации (Гц)
FftL=1024;% Количество линий Фурье спектра
T=1:130;% Массив отсчетов времени
Signal1=Signal(1:130);
%% Спектральное представление сигнала

FftS=abs(fft(Signal,FftL));% Амплитуды преобразования Фурье сигнала

%% Построение графиков
%subplot(2,1,1);% Выбор области окна для построения
plot(T, Signal1);% Построение сигнала
title('Сигнал');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда, В');% Подпись оси у графика
F=0:Fd/FftL:Fd/2;% Массив частот вычисляемого спектра Фурье
figure
%subplot(2,1,2);% Выбор области окна для построения
for k=1:length(F)
Z(k)=FftS(k);
if Z(k)>=1300
Z(k)=1300;
disp(k*Fd/FftL)
else
Z(k)=round(Z(k));
endif
Z(k);
end
k=1:length(F);
%plot(F, Z(k));% Построение спектра Фурье сигнала
stem(F, Z(k));% Построение спектра Фурье сигнала
%plot(F, FftS(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда, В');% Подпись оси у графика