【正文】 。[t,x]=ode45(39。q5*x(1)]輸入：ts=::10。q3*x(1)。f=[u*x(1)。q4=+u*。q2=+u/。subplot(2,2,4)fly(x,n+15)。subplot(2,2,2)fly(x,n+5)。)axis([26,26,4,1])。,x,y2,39。plot(x,y,x,z,x,y1,39。R=max(r)。y1=ones(1,length(x))*pi。y=fly1(x,n)。for i=1:n y=y+((1cos(i*pi))/(i^2*pi))*cos(i*x)+(((1)^(i+1))/i)*sin(i*x)。w=zeros(1,l)。3．function y=fly1(x,n)。for i=1:n y=y+((1cos(i*pi))/(i^2*pi))*cos(i*x)+(((1)^(i+1))/i)*sin(i*x)。w=zeros(1,l)。2．function y=fly(x,n)。 endendplot(x,y)。 if x(i)(2*m+1)*pi0 y(i)=0。 end end end end endend二題：1．function y=fd(x)n=length(x)。39。 disp(D(n,1:2))。D(n,m)~=0 if gjwtf(G(i,j),D(n,m)) disp(G(i,1:2))。k=1。 end end endendx=length(G(:,1))。 D(k,:)=LD(i,:)。 k=k+1。 end end endendk=1。 G(k,:)=LG(i,:)。 k=k+1。k=1。G=[]。l=length(L(1,:))。LG=L。39。for i=1:h for j=3:l if L(i,j)==a for k=j:l if L(i,k)==b f=1。l=length(L(1,:))。39。end 5.function[f]=gjwtf(a,b)f=0。39。163。206。180。200。202。 r=r+1。)。186。206。181。201。187。 disp(39。39。163。213。179。179。181。 disp(G(1,1:2))。39。163。180。179。179。203。 disp(r)。39。163。G(i,j)~=0 f=0。r=1。 end end endendx=length(G(:,1))。 D(k,:)=LD(i,:)。 k=k+1。 end end endendk=1。 G(k,:)=LG(i,:)。 k=k+1。f=1。D=[]