* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
Моделирование в системе Maple + MATLAB
1175
Рис. 15.51. Временная зависимость сигнала с шумом
> ft := fft(noisy_data): > VectorOptions(ft, datatype);
complex8 Эта операция переводит задачу из временного представления сигнала в час тотное, что позволяет использовать частотные методы анализа сигнала. Выделим, к примеру, действительную и мнимую части элементов вектора ft и проверим его размер:
> real_part := map(Re, ft): imag_part := map(Im, ft): > dimensions(ft);
[1500] Теперь выполним обычные операции вычисления спектра и зададим построе ние графика частотного спектра мощности сигнала:
> > > > setvar("FT", ft);setvar("n", num); evalM("result = FT.*conj(FT)/n"); pwr := getvar("result"): VectorOptions(pwr, datatype);
float8
> pwr_list := convert(pwr, list): > pwr_points := [seq([(t-1)/Time[num], pwr_list[t]], t=1..num/2)]: > plots[pointplot](pwr_points, style=line);
Спектр сигнала представлен на рис. 15.52. Из него отчетливо видно, что сигнал представлен двумя частотными состав ляющими с разной амплитудой. Они четко выделяются.
15.5.2. Моделирование линейного осциллятора
Теперь рассмотрим организацию совместной работы систем Maple 9 и MATLAB 6.5 SP1 (это широко распространенная, хотя и не новейшая версия данной системы) на примере моделирования механического осциллятора (маятника).