Не выводится вид окна Чебышева, в чем ошибка?

@DFF11F · Регистрация: 24.03.2012

Студворк — интернет-сервис помощи студентам

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
%ФУНКЦИЯ РАСЧЕТА СПМ МЕТОДОМ УЭЛЧА
function zadanie4(x)
nxt=1;
ur=input('ВВЕДИТЕ УРОВЕНЬ БОКОВЫХ ЛЕПЕСТКОВ =');
dlit=input('ВВЕДИТЕ ДЛИТЕЛЬНОСТЬ СИГНАЛА =');
while nxt>0
long=input('ВВЕДИТЕ ДЛИНУ ОКНА n=');
ur=ur; %#ok<ASGSL>
beta = cosh(acosh(10^(ur/20))/(long-1));
k = 0:long-1;
W = cos((long-1)*acos(beta*cos(pi*k/long))) ./ ...
    cosh((long-1)*acosh(beta));   %ОКНО ЧЕБЫШЕВА
 
W = W.*exp(1i*pi/long * k);
w = real(fft(W)); 
 
midl = long/2+1;
w = w/w(2);
w = [fliplr(w(2:midl)) w(2:midl)]';
 
plot(w)
xlabel('Число отсчётов');
ylabel('Амплитуда');
grid on
 
for ii=1:3
        subplot(3,2,(2*ii-1))
        pwelch(x(1:dlit,ii),boxcar(dlit),[],[],100);     %СПМ БЕЗ ВЫДЕЛЯЮЩЕЙ ФУНКЦИИ
        title(['СПМ ',num2str(ii),'-ГО КАНАЛА БЕЗ ВЫДЕЛЯЮЩЕЙ ФУНКЦИИ'],'FontName','Times New Roman');
        subplot(3,2,2*ii)
        pwelch(x(1:dlit,ii),w,[],[],100);                %СГЛАЖЕНО С ПОМОЩЬЮ ОКНА ЧЕБЫШЕВА
        title(['СПМ ',num2str(ii),'-ГО КАНАЛА С ИСПОЛЬЗОВАНИЕМ ОКНА ЧЕБЫШЕВА'],'FontName','Times New Roman');
        xlabel('НОМЕР ОТСЧЕТА');
        ylabel('ЗНАЧЕНИЕ СПМ');
     grid on
    end;
    nxt=input('ИЗМЕНИТЬ ДЛИНУ ОКНА И ПРОДОЛЖИТЬ РАСЧЕТ? (1-ДА, 0-НЕТ)  ');
end;
MENU

@vital792 · 30.04.2012, 00:35

DFF11F, все работает, если вводить разумные значения(в коде отсутствуют проверки на корректность вводимых данных). Вот таким кодом я тестировал твою функцию:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
function main
clear; clc;
s = zeros(1024, 3);
s(:, 1) = sin(2*pi*100*(1:1024)/1000);
s(:, 2) = sin(2*pi*200*(1:1024)/1000);
s(:, 3) = sin(2*pi*300*(1:1024)/1000);
zadanie4(s);
end
 
%ФУНКЦИЯ РАСЧЕТА СПМ МЕТОДОМ УЭЛЧА
function zadanie4(x)
nxt=1;
ur=input('ВВЕДИТЕ УРОВЕНЬ БОКОВЫХ ЛЕПЕСТКОВ =');
dlit=length(x);%input('ВВЕДИТЕ ДЛИТЕЛЬНОСТЬ СИГНАЛА =');
while nxt>0
long=input('ВВЕДИТЕ ДЛИНУ ОКНА n=');
ur=ur; %#ok<ASGSL>
beta = cosh(acosh(10^(ur/20))/(long-1));
k = 0:long-1;
W = cos((long-1)*acos(beta*cos(pi*k/long))) ./ ...
    cosh((long-1)*acosh(beta));   %ОКНО ЧЕБЫШЕВА
 
W = W.*exp(1i*pi/long * k);
w = real(fft(W)); 
 
midl = long/2+1;
w = w/w(2);
w = [fliplr(w(2:midl)) w(2:midl)]';
 
plot(w)
xlabel('Число отсчётов');
ylabel('Амплитуда');
grid on
 
for ii=1:3
        subplot(3,2,(2*ii-1))
        pwelch(x(1:dlit,ii),boxcar(dlit),[],[],100);     %СПМ БЕЗ ВЫДЕЛЯЮЩЕЙ ФУНКЦИИ
        title(['СПМ ',num2str(ii),'-ГО КАНАЛА БЕЗ ВЫДЕЛЯЮЩЕЙ ФУНКЦИИ'],'FontName','Times New Roman');
        subplot(3,2,2*ii)
        pwelch(x(1:dlit,ii),w,[],[],100);                %СГЛАЖЕНО С ПОМОЩЬЮ ОКНА ЧЕБЫШЕВА
        title(['СПМ ',num2str(ii),'-ГО КАНАЛА С ИСПОЛЬЗОВАНИЕМ ОКНА ЧЕБЫШЕВА'],'FontName','Times New Roman');
        xlabel('НОМЕР ОТСЧЕТА');
        ylabel('ЗНАЧЕНИЕ СПМ');
     grid on
    end;
    nxt=input('ИЗМЕНИТЬ ДЛИНУ ОКНА И ПРОДОЛЖИТЬ РАСЧЕТ? (1-ДА, 0-НЕТ)  ');
end;
end
% MENU

Добавлено через 4 минуты
зы надеюсь не забыл, что длина окна должна быть четной? Я об этом уже писал, если не четная, оконная функция вычисляется немного по другому.

@DFF11F · 30.04.2012, 14:47 **[ТС]**

vital792, Не забыл, делаю для различных окон (четной длины 32,64,128 и 256) Выводятся графики СПМ, а вот сам вид окна не выводится... не знаю почему.

Добавлено через 1 минуту
vital792, Нужно, чтобы программа выдавала 2 графика figure 1 - с видом окна и figure 2 непосредственно с анализом реализаций.

@vital792 · 30.04.2012, 15:07

Сообщение от DFF11F

Нужно, чтобы программа выдавала 2 графика figure 1 - с видом окна и figure 2 непосредственно с анализом реализаций.

так поставь перед for команду figure

Matlab M
1
2
3
4
5
6
7
8
9
...
ylabel('Амплитуда');
grid on
 
 figure
 
for ii=1:3
        subplot(3,2,(2*ii-1))
...

@DFF11F · 01.05.2012, 17:36 **[ТС]**

vital792, Неа, не кажет... Пустой первый график, а со вторым всё ОК

@vital792 · 01.05.2012, 20:38

Сообщение от DFF11F

Пустой первый график

какие параметры вводишь? Я вводил уровень лепестков 100, длину окна 64(и другие близкие параметры пробовал) все ок

@DFF11F · 02.05.2012, 22:17 **[ТС]**

vital792, Такие и ввожу, 32 и 64... Второй график - все ОК, а первый пуст... Не знаю почему, уровень ввожу так же 100, как бы и было при стандартном.

Добавлено через 39 минут
vital792, Все, заработало! Спасибо за помощь.

@vital792 · 02.05.2012, 22:18

хм, странно. Тогда можно еще так попробовать:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
...
figure(1)
plot(w)
xlabel('Число отсчётов');
ylabel('Амплитуда');
grid on
 
for ii=1:3
        figure(2)
        subplot(3,2,(2*ii-1))
...

У меня все ок

@vital792 · 02.05.2012, 22:20

не успел. Ну ладно, хорошо что получилось. Будет еще вариант)

@DFF11F · 03.05.2012, 08:58 **[ТС]**

vital792, Спасибо! Как знать, может пригодится!

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .

@DFF11F 29 / 24 / 1 Регистрация: 24.03.2012 Сообщений: 127
	30.04.2012, 14:47 [ТС]
	vital792, Не забыл, делаю для различных окон (четной длины 32,64,128 и 256) Выводятся графики СПМ, а вот сам вид окна не выводится... не знаю почему. Добавлено через 1 минуту vital792, Нужно, чтобы программа выдавала 2 графика figure 1 - с видом окна и figure 2 непосредственно с анализом реализаций. 0

@DFF11F 29 / 24 / 1 Регистрация: 24.03.2012 Сообщений: 127
	01.05.2012, 17:36 [ТС]
	vital792, Неа, не кажет... Пустой первый график, а со вторым всё ОК 0

@DFF11F 29 / 24 / 1 Регистрация: 24.03.2012 Сообщений: 127
	02.05.2012, 22:17 [ТС]
	vital792, Такие и ввожу, 32 и 64... Второй график - все ОК, а первый пуст... Не знаю почему, уровень ввожу так же 100, как бы и было при стандартном. Добавлено через 39 минут vital792, Все, заработало! Спасибо за помощь. 0

@vital792 2014 / 1286 / 61 Регистрация: 05.06.2010 Сообщений: 2,213
	02.05.2012, 22:20
	не успел. Ну ладно, хорошо что получилось. Будет еще вариант) 1

@DFF11F 29 / 24 / 1 Регистрация: 24.03.2012 Сообщений: 127
	03.05.2012, 08:58 [ТС]
	vital792, Спасибо! Как знать, может пригодится! 0