yzfmakegrid.m

来自「能方便的矩阵进行行列的线性扩展,并把生的矩阵存储也很方便」· M 代码 · 共 33 行

function kk=yzfMakeGrid(outfn)
cenn=30;
thetan=30;

r=-1:1/cenn:1;
r=r';
theta=0:2*pi/thetan:2*pi;
y=sin(theta);
x=cos(theta);
y=y';
x=x';
[xr,xc]=size(r);
[tr,tc]=size(y);
    z1=y.*r(1);
    z2=x.*r(1);
    kk=[z2 z1];
    for m=2:xr;
        z1=y.*r(m);
	    z2=x.*r(m);
        tt=[z2 z1];
        tt=tt';
        kk=kk';
        kk=[kk tt];
        kk=kk';
    end;
    [r,c]=size(kk);
    rt=rand(r,1);
    kk=[kk rt];
    kk=kk';
    fo=fopen(outfn,'wt');
    fprintf(fo,'%f\t%f\t%f\n',kk);
    fclose(fo);