MATLAB: Plotting 3D grating

3d gratingrectangular pulse

Considering the plots below, I have a problem in adding an extra dimension to a 1D rectangular pulse train. In other words, I am trying to implement a 3D grating in matlab, as shown below. Any help would be appreciated!

Best Answer

I think this is what you're looking for

x = [1 2 2 3 3 4 4 5 5 6];
z = [0 0 1 1 0 0 1 1 0 0];
y = [1 4];
[X,Y] = meshgrid(x,y)
Z = repmat(z,[length(y) 1])
surf(X,Y,Z,'FaceColor',[0.5 0.5 0.5])

Related Solutions

MATLAB: Using fft with rectangularPulse

It is now:

f=@(t) rectangularPulse(-1,1,t);
time=-6:0.01:6;
Fs = 1/0.01;
Fn = Fs/2;
Freq = linspace(-Fn, Fn, numel(time));
pulse=f(time);
fourier=fft(pulse);
figure()
subplot(3,1,1)
plot(time,pulse)                                % Pulse
subplot(3,1,2)
plot(Freq,fftshift(real(fourier)));             % Magnitude of Fourier
xlim([-10 10])                                  % (Optional)

subplot(3,1,3)
plot(Freq,fftshift(unwrap(angle(fourier))));    % Phase of Fourier
xlim([-10 10])                                  % (Optional)

Use fftshift to centre the amplitude and phase plots.

I also calculated the frequency vector for the ‘fourier’ plots.

MATLAB: Can you help me? :(

This demo shows how to find the fourier and inverse fourier transforms

Fourier transform:

syms t w
y(t) = rectangularPulse(t);
F(w) = int(y(t)*exp(-1i*w*t), t, -0.5, 0.5); % pulse only exist between -0.5 to 0.5
figure;
fplot(F, [-100 100]);

Inverse:

F_inv(t) = int(F(w)*exp(-1i*w*t), w, -100, 100); % integrate frequencies between -100 to 100
t_val = linspace(-2, 2, 1000);
F_val = F_inv(t_val);
figure;
plot(t_val, F_val)

Best Answer

Related Solutions

MATLAB: Using fft with rectangularPulse

MATLAB: Can you help me? :(

Related Question