• MATLAB画的准周期结构
    % Quasiperiodic clc;clear;close all; A=pi/2+[0:pi/7:[pi*2-pi/8]]; z1=exp(i*A); hold on; LD=2; for k=1:length(z1);     zt=[0,z1(k)];     plot(real(zt),imag(zt),'b','Linewidth',LD); end z2=[z1+[z1(2:end),z1(1)]]; N=length(z2); for k=1:N;     kt=mod(k-1,N)+1;     zta=[z2(kt),z1(kt)];     kt=mod(k,N)+1;     ztb=[z2(k),z1(kt)];     plot(real(zta),imag(zta),'r','Linewidth',LD);     plot(real(ztb),imag(ztb),'r','Linewidth',LD); end
    2011-05-11 20:41