1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
// 08.09.15
// 08.10.10
// 08.10.21
// 09.09.10
function FsL=Wirerawparadata(varargin)
timer();
Nargs=length(varargin);
Fbdy=Projpara(varargin(1));
Fd=varargin(2); N=3;
FdL=Fullformfunc(Fd);
Fxy=Mixop(1,FdL);
Tmp=Mixop(5,FdL);
K=mtlb_findstr(Tmp,'=');
Uname=part(Tmp,1:K-1);
Urg=evstr(part(Tmp,K+1:length(Tmp)));
Tmp=Mixop(6,FdL);
K=mtlb_findstr(Tmp,'=');
Vname=part(Tmp,1:K-1);
Vrg=evstr(part(Tmp,K+1:length(Tmp)));
DuL=varargin(N);
DvL=varargin(N+1);
N=N+2;
Np=[50,50];
if Nargs>=N
Np=varargin(N);
if type(Np)==1 & length(Np)==1
Np=[Np,Np];
end;
end
Eps=0.05;
if Nargs>=N+1
Eps=varargin(N+1);
end;
Eps2=0.2;
if Nargs>=N+2
Eps2=varargin(N+2);
end;
Umin=Urg(1); Umax=Urg(2);
Vmin=Vrg(1); Vmax=Vrg(2);
if type(DuL)==10
DuL=evstr(DuL);
elseif length(DuL)==1
Tmp=(Umax-Umin)/(DuL+1);
Tmp1=Umin+Tmp;
Tmp2=Umin+DuL*Tmp
DuL=[Tmp1:Tmp:Tmp2];
end;
if type(DvL)==10
DvL=evstr(DvL);
elseif length(DvL)==1
Tmp=(Vmax-Vmin)/(DvL+1);
Tmp1=Vmin+Tmp;
Tmp2=Vmin+DvL*Tmp
DvL=[Tmp1:Tmp:Tmp2];
end;
Fbdxy=Makexybdy(Fd,Np);
FsL=[];
for I=1:length(DuL)
U=DuL(I);
Tmp=Paramplot('[U,V]','V=[Vmin,Vmax]','N=Np(1)');
FkL=Cuspsplitpara(Tmp,Fd,Eps);
for J=1:Mixlength(FkL)
Fk=Mixop(J,FkL);
FsL=Mixadd(FsL,Fk);
end;
end;
Nx=Mixlength(FsL);
if Nx>0
Mixdisp(Mix(Uname,'-data=',1,':',Nx,' Time=',timer()));
end;
for I=1:length(DvL)
V=DvL(I);
Tmp=Paramplot('[U,V]','U=[Umin,Umax]','N=Np(2)');
FkL=Cuspsplitpara(Tmp,Fd,Eps);
for J=1:Mixlength(FkL)
Fk=Mixop(J,FkL);
FsL=Mixadd(FsL,Fk);
end;
end;
Ny=Mixlength(FsL)-Nx;
if Ny>0
Mixdisp(Mix(Vname,'-data=',Nx+1,':',Nx+Ny,' Time=',timer()));
end;
endfunction;
|