lingmindufenxi.m

#include "stdio.h" #include "math.h" #include "malloc.h" #include "string.h" #define m 6 #defin

x=[0.075,0.375,0.123,0.115,1.273,12,0.771];
dh=0.0001;
l=41;
x1=ones(l,7);
x2=x1;
x3=x1;
x4=x1;
x5=x1;
x6=x1;
x7=x1;
for i=1:l
    x1(i,:)=x;
    x2(i,:)=x;
    x3(i,:)=x;
    x4(i,:)=x;
    x5(i,:)=x;
    x6(i,:)=x;
    x7(i,:)=x;
end
for i=-(l-1)/2:1:(l-1)/2
    x1(i+(l-1)/2+1,1)=x1(i+(l-1)/2+1,1)+dh*i;

    x2(i+(l-1)/2+1,2)=x2(i+(l-1)/2+1,2)+dh*i;

    x3(i+(l-1)/2+1,3)=x3(i+(l-1)/2+1,3)+dh*i;

    x4(i+(l-1)/2+1,4)=x4(i+(l-1)/2+1,4)+dh*i;

    x5(i+(l-1)/2+1,5)=x5(i+(l-1)/2+1,5)+dh*i;

    x6(i+(l-1)/2+1,6)=x6(i+(l-1)/2+1,6)+dh*i;

    x7(i+(l-1)/2+1,7)=x7(i+(l-1)/2+1,7)+dh*i;
end
A=0.01/3;
B=0.05/3;
C=0.1/3;
D=[B,B,B,C,C,B,B];
for i=1:l
    z1(i)=Czfee(x1(i,:),D);
    z2(i)=Czfee(x2(i,:),D);
    z3(i)=Czfee(x3(i,:),D);
    z4(i)=Czfee(x4(i,:),D);
    z5(i)=Czfee(x5(i,:),D);
    z6(i)=Czfee(x6(i,:),D);
    z7(i)=Czfee(x7(i,:),D);
end
x=-(l-1)/2*dh:dh:(l-1)/2*dh;
V=[-0.002,0.002,400000,600000];
axis(V);
plot(x,z1,'r');
hold on;
plot(x,z2,'black');
hold on;
plot(x,z3,'m'); 
hold on;
plot(x,z4,'b');
hold on;
plot(x,z5,'c');
hold on;
plot(x,z6,'g');
hold on;
plot(x,z7,'y');