readpch.m

来自「读pch文件的小程序」· M 代码 · 共 41 行

clear
clc
fid = fopen('maax.pch');

% begin to read the first line %
tline = fgetl(fid);
name = tline(9:12);
ifo = str2num(tline(24));
Ndof = str2num(tline((length(tline)-4):length(tline)));
% end of read the first line %

lie = 0;
matrix_lie = zeros(2+Ndof/6,6);
hang = 0;
KK = zeros(Ndof,Ndof);

while ~feof(fid)
    tline= fgetl(fid);
    if tline(1)=='*' 
        Grid = str2num(tline(18:24));
        Component = str2num(tline(25:40));
        KKKK = str2num(tline(41:length(tline)));
        hang = matrix_lie(Grid,Component);
        KK(hang,lie) = KKKK;       
%     elseif tline(1)=='*' & lie>Ndof
%         Grid = str2num(tline(18:24));
%         Component = str2num(tline(25:40));
%         KKKK = str2num(tline(41:length(tline)));
%         hang = matrix_lie(Grid,Component);
%         MM(hang,lie) = KKKK;       
    else
        lie = lie + 1;
        Grid = str2num(tline(24:40));
        Component = str2num(tline(45:length(tline)));
        matrix_lie(Grid,Component) = lie;
    end
   end
fclose(fid);
eval([name,'= KK(1:hang,1:hang) + KK(1:hang,1:hang)''- diag(diag(KK(1:hang,1:hang)));']);
eval(['save',' ','data',' ',name])