compeleteinsert.m

能方便的矩阵进行行列的线性扩展,并把生的矩阵存储也很方便

function out=yzfinsertandsave(fin,fout,rowinsert,colinsert)
%需要给定的参数值
%finame='d:\in.txt';
%foname='d:\out.txt';
%rowinsert 每两行要插几行
%colinsert 每两列要插几列

%从文件fin加载数组fin
fin=load(fin);
fin=[fin fin(:,1)]; %组成环

%确保内环和外环对齐
[rownum,colnum]=size(fin);
finout=fin*GetInsertColumnMatrix(colnum,1);
finout=finout(:,1:end-1);
finout1=finout(1,1:end);
finout2=finout(2,1:end-1);
finout2=[finout(2,end),finout2];
fin=[finout1;finout2];
fin=[fin fin(:,1)] %组成环

%行列插值
[rownum,colnum]=size(fin);
out=GetInsertLayerMatrix(rownum,rowinsert)*fin*GetInsertColumnMatrix(colnum,colinsert);
out=out(:,1:end-1)
[rownum,colnum]=size(out);

%生成文件
foutc=[fout '.txt'];
foutm=[fout '.dat'];

fpoc=fopen(foutc,'wt');
fpom=fopen(foutm,'wt');

fprintf(fpoc,'%d\t%d\n',rownum,colnum);%记录行数，列数


for k=1:rownum;
   for j=1:colnum-1;
    fprintf(fpoc,'%f\t',out(k,j));
    fprintf(fpom,'%f\t',out(k,j));
   end;
    fprintf(fpoc,'%f\n',out(k,j+1));
    fprintf(fpom,'%f\n',out(k,j+1));
end;

fclose(fpoc);
fclose(fpom);