how to plot this equation ?

1 次查看（过去 30 天）
Mohamed Ahmed
Mohamed Ahmed 2022-1-25
编辑： Mohamed Ahmed 2022-1-31
I'm trying to plot the equation shown in the image, but I'm not sure how to apply the Galois Seq. "gf()" here.
Can you please help me with this ?
Expected magnitude waveform.
clear all
close all
clc
%%
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sampling period
L = 1024; % Length of signal
t = 2*(0:L-1)*T; % Time vector
signal = 0;
m = 5;
for k=1:((2^m)-1)
g = gf(2^m);
signal = signal + g(k).*exp(1i*2*pi*k*t)/((2^m)-1);
end
figure
subplot(2,2,1)
plot(t,signal)
title('signal')

请先登录，再回答此问题。

采纳的回答

Sulaymon Eshkabilov
Ran in:
Here you'd need to define the purpose of using GF.
If you are aiming to select any randomly selected values of GF while computing the wave form fcn values, this can be done with the followings:
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sampling period
L = 1024; % Length of signal
t = 2*(0:L-1)*T; % Time vector
signal = 0;
m = 5;
x = [0 1 1; 0 1 0; 1 1 1];
for k=1:((2^m)-1)
g =gf(x,m);
G = g.x(:);
signal = signal + double(G(randi(numel(G))))*exp(1i*2*pi*k*t)/((2^m)-1);
end
figure
subplot(2,2,1)
plot(t,real(signal))
title('signal: RE')
subplot(2,2,2)
plot(t,imag(signal))
title('signal: IM')
  3 个评论
显示 1更早的评论
Sulaymon Eshkabilov
Sulaymon Eshkabilov 2022-1-30
If you want to have different values of x, then change the values of x or if you want to have specific values for x, then speficy them.
See this DOC for more info on x.
Mohamed Ahmed
Mohamed Ahmed 2022-1-31
编辑：Mohamed Ahmed 2022-1-31
Should they all be changed to the same values? and how will that affect the signal ?
And what are these spikes in the middle of the signal ?

请先登录，再进行评论。

更多回答（1 个）

Benjamin Thompson
Benjamin Thompson 2022-1-26
This lines needs an extra set of parenthesis:
signal = signal + g(k).*exp((1i*2*pi*k*t)/((2^m)-1));
But the way you are using the gf function does not match how MATLAB describes it in the documentation.

类别

在 Help Center 和 File Exchange 中查找有关 Digital Filter Analysis 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by