如何在matlab中繪製這個吉布斯傅里葉圖

Question

請有人可以幫助我，我需要用gibbs振盪來繪製這個波閉合的牙齒，我有一天能做到這一點，我從來沒有用Matlab進行圖形繪製。

這是我嘗試編程吉布斯振盪波圖：http://imageshack.us/a/img585/9724/kakai12.jpg

這是我做了什麼，但不是真的，如果它的工作這種方式。

clear all; 
syms x n; 
f=(-4)*(x-pi); 
T=2; w0=((2*pi)/T); 
fb=f*sin(n*w0*x); 
bn=(2/x)*int(fb,x,0,2*pi); 
x=0:1/(1000):15; 
bn_12=bn*sin(n*w0*x); 
c=1; j=100; 
fs=symsum(bn_12,n,c,j) 

請有人幫忙。

Answer 1

現在不能測試此權利，但這樣的事情應該工作：

t = -4:0.01:4 %Create numerical vector of time 
T = 2; %Period, adjust as suitable 
V = 4; %Adjust to appropriate amplitude 
w0 = 2*pi/T; %Angular frequency 

line = (V/pi) * w0 * t - V; %Start with straight line with correct gradient and y intercept 
%Then subtract 2*V for every 2*pi 
step = floor((w0 * t)/(2*pi))*2*V; %So step should be a step function that steps up 2V every 2pi 
sawtooth = line - step; 

plot(w0*t, sawtooth)