Міністерство освіти і науки України Національний університет „Львівська політехніка” Кафедра ЕОМ Звіт лабораторної роботи № 2 з дисципліни: «Проектування комп’ютерних засобів обробки сигналів та зображень» на тему: «Діагностика роботи цифрових фільтрів шляхом аналізу їх амплітудно-частотної характеристики» Мета роботи Дослідити і проаналізувати параметри амплітудно-частотної характеристики та вплив віконної обробки при спектральному аналізі сигналів. Варіант завдання до лабораторної роботи № Вар N L Sm Sp S A № Вагової функції  59 28 124 -32 44 8 14 15   Номер Функції Назва Тип функції  15  w(n) = 0.25 + 0,75 cos [ ((n-16)/32]   Хід роботи 1. Формуємо масив синусоїдальної і косинусоїдальної складової згідно з виразом і згідно завдання:
for(i = 0; i < N; i++) { sin0[i] = sin(2*pi*i*l/N); cos0[i] = cos(2*pi*i*l/N); } 2. Сформувати вхідний масив (синусоїдальна і косинусоїдальна складові) згідно з формулою 3.
де
for (n=0;n<=Sp-Sm+1;n++) for(i = 0; i <= N; i++) { sin1[n][i] =Asin[2*pi(Sl+Sm+n)i/SN]; cos1[n][i] =Acos[2*pi(Sl+Sm+n)i/SN]; } 3. Сформувати масив вагової функції W(i) згідно варіанту завдання. Для вагової функції №42 обчислення виконується згідно виразу:
4. Скласти процедуру на мові високого рівня для обчислення АЧХ згідно з формулою 2
a) з ваговою функцією for (n=0;n<=Sp-Sm+1;n++) for(i = 0; i <= N; i++) { Y0[n]=Y0[n]+W[і]*(Acos1[n][i]*cos0[i]+Asin1[n][i]*sin0[i]) Y1[n]=Y1[n]+W[і]*(Asin1[n][i]*cos0[i]-Acos1[n][i]*sin0[i]) } b) без вагової функції for (n=0;n<=Sp-Sm+1;n++) for(i = 0; i <= N; i++) { Y0[n]=Y0[n]+(Acos1[n][i]*cos0[i]+Asin1[n][i]*sin0[i]) Y1[n]=Y1[n]+(Asin1[n][i]*cos0[i]-Acos1[n][i]*sin0[i]) } 5. Скласти процедуру графічного виводу значень АЧХ з використанням вагової функції і без її використання та порівняти отримані результати. Значення координати Y для відображення результату обчислюємо за формулою: Y[n] = sqrt(pow(Y0[n][l],2) + pow(Y1[n][l],2)); Результати виконання індивідуального завдання
Рис.1. Результат виконання програми Лістинг Програми using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms; using ZedGraph; using System.Numerics; namespace ll2 { public partial class Form1 : Form { int N = 0; public Form1() { InitializeComponent(); } private void zedGraph_Load(object sender, EventArgs e) { // DrawGraph(); } private void DrawGraph() { ZedGraph.MasterPane masterPane = zedGraph.MasterPane; masterPane.PaneList.Clear(); GraphPane pane = new GraphPane(); GraphPane pane1 = new GraphPane(); // GraphPane pane2 = new GraphPane(); int labelsXfontSize = 25; int labelsYfontSize = 20; int titleXFontSize = 25; int titleYFontSize = 20; int legendFontSize = 25; int mainTitleFontSize = 45; masterPane.Add(pane); masterPane.Add(pane1); // Установим размеры шрифтов для меток вдоль осей pane.XAxis.Scale.FontSpec.Size = labelsXfontSize; pane.YAxis.Scale.FontSpec.Size = labelsYfontSize; // Установим размеры шрифтов для подписей по осям pane.XAxis.Title.FontSpec.Size = titleXFontSize; pane.YAxis.Title.FontSpec.Size = titleYFontSize; pane.Legend.FontSpec.Size = legendFontSize; pane.Title.FontSpec.Size = mainTitleFontSize; pane.Title.FontSpec.IsUnderline = true; pane.XAxis.Title.Text = "X"; pane.YAxis.Title.Text = "Y"; // Установим размеры шрифтов для меток вдоль осей pane1.XAxis.S...