Подскажите как исправить ошибку в коде, чтобы корректно выводило графики

@Zenuk18 · Регистрация: 12.08.2016

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

У меня был пример кода как построить график с 2 отрезков, а мне нужно 3 отрезков как на скриншоте, пробовал переделать код, но постоянно выбивает ошибки

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
50
51
52
clc
clear 
close, close
T=12;m=8;N=2^m;
p=N/32;k1=N/2-p; k2=N/2+p;
dt=T/N;
t1=0:dt:4-dt;
t2=4:dt:8-dt;
t3=8:dt:12-dt;
t=[t1,t2,t3];
x1=4+0*t1;
x2=8+0*t2;
x3=4+0*t3;
x=[x1,x2,x3];
dw=2*%pi/T;
w=1:1:N;
Sx = fft(x)/(N);
figure(1,'BackgroundColor',[1,1,1]);
subplot(3,1,1),   plot(t, x),          title('Вхідний сигнал - x(n)','position',[T/2,max(x)/2]);
subplot(3,1,2),   plot(w, real(Sx)),   title('Дійсна частина ШПФ','position',[N/2,max(real(Sx))/2]);
subplot(3,1,3),   plot(w, imag(Sx)),   title('Уявна частина ШПФ','position',[N/2,max(imag(Sx))/2]);
Sx1 = [Sx(1 : k1) zeros(1, k2 - k1) Sx(k2 + 1 : N)];
Sx2 = [zeros(1, p) Sx(p+1 :N-p) zeros(1, p)];
Sx3 = [Sx(1 : N/4-p/2),zeros(1, p),Sx( (N/4+p/2)+1:3*N/4-p/2), zeros(1, p), Sx(3*N/4+p/2+1:N)];
Sx4 = [zeros(1,  N/4-p/2),Sx((N/4-p/2)+1 :N/4+p/2), zeros(1, N/2-p), Sx(3*N/4-p/2+1:3*N/4+p/2),zeros(1,N/4-p/2)];
figure(2,'BackgroundColor',[1,1,1]);
subplot(5,2,1),plot(w, abs(Sx)),title('Амплітудний спектр X(k)','position',[N/4,max(abs(Sx))/2]);;
subplot(5,2,3),plot(w, abs(Sx1)),title('Модифікований амплітудний спектр Xm1(k)','position',[N/20,max(abs(Sx1))/2]);
subplot(5,2,5),plot(w, abs(Sx2)),title('Модифікований амплітудний спектр Xm2(k)','position',[N/20,max(abs(Sx2))/2]);
subplot(5,2,7),plot(w, abs(Sx3)),title('Модифікований амплітудний спектр Xm3(k)','position',[N/20,max(abs(Sx3))/2]);
subplot(5,2,9),plot(w, abs(Sx4)),title('Модифікований амплітудний спектр Xm4(k)','position',[N/20,max(abs(Sx4))/2]);
x1= N*ifft(Sx);
xm1= N*ifft(Sx1);
xm2= N*ifft(Sx2);
xm3= N*ifft(Sx3);
xm4= N*ifft(Sx4);
subplot(5, 2, 2),   plot(t, x1),         title('Відновлений сигнал x1(n)','position',[T/4,max(abs(x1))/2]);
subplot(5, 2, 4),   plot(t, real(xm1)),  title('Відновлений сигнал xm1(n)','position',[T/4,max(abs(xm1))/2]);
subplot(5, 2, 6),   plot(t, real(xm2)),  title('Відновлений сигнал xm2(n)','position',[T/4,max(abs(xm2))/2]);
subplot(5, 2, 8),   plot(t, real(xm3)),  title('Відновлений сигнал xm3(n)','position',[T/4,max(abs(xm3))/2]);
subplot(5, 2, 10), plot(t, real(xm4)),   title('Відновлений сигнал xm4(n)','position',[T/4,max(abs(xm4))/2]);
onePr=100/(N*max(abs(8),abs(4))); 
bm1=onePr*(sum(abs(x-xm1)))
bm2=onePr*(sum(abs(x-xm2))) 
bm3=onePr*(sum(abs(x-xm3)))
bm4=onePr*(sum(abs(x-xm4))) 
disp("N=",N)
disp("p=",p)
disp("bm1=",bm1)
disp("bm2=",bm2)
disp("bm3=",bm3)
disp("bm4=",bm4)

@АлександрКом · 25.02.2021, 11:24

Zenuk18, попробуйте t1 сделать ровно до 4

@Krasme · 25.02.2021, 11:32

строка 9

Code
1
t3=8:dt:12;

Новые блоги и статьи Все статьи Все блоги /
Мысли в слух kumehtar 07.11.2025 Заметил среди людей, что по-настоящему верная дружба бывает между теми, с кем нечего делить.	Новая зверюга volvo 07.11.2025 Подарок на Хеллоуин, и теперь у нас кроме Tuxedo Cat есть еще и щенок далматинца: Хочу еще Симбу взять, очень нравится. . .	Инференс ML моделей в Java: TensorFlow, DL4J и DJL Javaican 05.11.2025 Python захватил мир машинного обучения - это факт. Но когда дело доходит до продакшена, ситуация не так однозначна. Помню проект в крупном банке три года назад: команда data science натренировала. . .	Mapped types (отображённые типы) в TypeScript Reangularity 03.11.2025 Mapped types работают как конвейер - берут существующую структуру и производят новую по заданным правилам. Меняют модификаторы свойств, трансформируют значения, фильтруют ключи. Один раз описал. . .	Адаптивная случайность в Unity: динамические вероятности для улучшения игрового дизайна GameUnited 02.11.2025 Мой знакомый геймдизайнер потерял двадцать процентов активной аудитории за неделю. А виновником оказался обычный генератор псевдослучайных чисел. Казалось бы - добавил в карточную игру случайное. . .
Протоколы в Python py-thonny 31.10.2025 Традиционная утиная типизация работает просто: попробовал вызвать метод, получилось - отлично, не получилось - упал с ошибкой в рантайме. Протоколы добавляют сюда проверку на этапе статического. . .	C++26: Read-copy-update (RCU) bytestream 30.10.2025 Прошло почти двадцать лет с тех пор, как производители процессоров отказались от гонки мегагерц и перешли на многоядерность. И знаете что? Мы до сих пор спотыкаемся о те же грабли. Каждый раз, когда. . .	Изображения webp на старых x32 ОС Windows XP и Windows 7 Argus19 30.10.2025 Изображения webp на старых x32 ОС Windows XP и Windows 7 Чтобы решить задачу, использовал интернет: поисковики Google и Yandex, а также подсказки Deep Seek. Как оказалось, чтобы создать. . .	Passkey в ASP.NET Core identity stackOverflow 29.10.2025 Пароли мертвы. Нет, серьезно - я повторяю это уже лет пять, но теперь впервые за это время чувствую, что это не просто красивые слова. В . NET 10 команда Microsoft внедрила поддержку Passkey прямо в. . .	Последние результаты исследования от команды MCM (октябрь 2025 г.) Programma_Boinc 29.10.2025 Последние результаты исследования от команды MCM (октябрь 2025 г. ) Поскольку мы продолжаем изучать гены, которые играют ведущую роль в развитии рака, в рамках проекта "Картирование раковых. . .

Подскажите как исправить ошибку в коде, чтобы корректно выводило графики

Решение

Решение

@АлександрКом 1298 / 927 / 449 Регистрация: 21.10.2012 Сообщений: 2,604
	25.02.2021, 11:24
	Сообщение было отмечено Zenuk18 как решение Решение Zenuk18, попробуйте t1 сделать ровно до 4 1

Опции темы