有關數(shù)據(jù)讀入matlab

來源：edatop 更新時間：2024-07-12 閱讀：

hfss方法CST前些日子在論壇見不少人都對讀入數(shù)據(jù)有些迷茫
cst、hfss的數(shù)據(jù)基本格式是描述性語言＋數(shù)據(jù)的格式
下面是一個cst的數(shù)據(jù)輸出
       Frequency / GHz        SZmin(1),Zmin(1)/abs,linear
-----------------------------------------------------------------
                         5                           0.4557996
                 5.029                           0.4553181
                 5.058                           0.4547905
                 5.087                           0.4542164
                 5.116                           0.4535954
                 5.145                           0.452927
                 5.174                           0.4522107
                 5.203                           0.4514461
但是matlab沒有這么高的領悟了，不會自己讀數(shù)據(jù)
方案不少說一說常用的。
1自己手動刪除
       Frequency / GHz        SZmin(1),Zmin(1)/abs,linear
-----------------------------------------------------------------
然后使用
load (文件名).txt
不錯的方法，如果你文件名是純數(shù)字的，比如是24，在matlab的workspace中會有一個X24的數(shù)組。
下面對x24處理就是了。
2一些人比較懶，不愿意刪東西
可以這樣：
[f,tr]=textread('24.txt','%f%f',2);
先到上面在看一下cst的數(shù)據(jù)結構，這里的2 是說從第3行開始讀數(shù)。
如果你的數(shù)據(jù)是大于兩列的，有幾列，在[a1,a2,a3……]就是幾個向量，同時，"%f%f……"中的%f的個數(shù)和前面的列向量的個數(shù)相同。
3有些時候儀器的輸出比較牛，我們的古董級別的失網(wǎng)的輸出是csv格式的帶編碼的數(shù)據(jù)。
# Version 1.00
#
freq[Hz];re:Trc1_b2d1sam[V];im:Trc1_b2d1sam[V];
1.700000000000000E+011;7.317389423714706E-003;-3.586198237922638E-003;
1.701500000000000E+011;-8.718488268764849E-003;5.495278251522440E-003;
1.703000000000000E+011;-1.050072140670953E-002;3.472822925688487E-003;
可以這么處理
fid=fopen('C:Documents and SettingsAdministratorMy Documents58.csv','r');% open file
i=0;
while 1                                                       %get the length of the data
    tline = fgetl(fid);
    if ischar(tline),   break,   end
    i=i+1;
end
fclose(fid);                                              %the array can't go back, when it get the bottom of the file. so open the file again
fid=fopen('C:Documents and SettingsAdministratorMy Documents58.csv','r');
temp=fgetl(fid);                                        % due to the structure of the data, first three line is useless.
temp=fgetl(fid);
temp=fgetl(fid);
for j=1:1:i-3                                               %read the data
    temp=fgetl(fid);
    a=strrep(temp,';',' ');                               % use the 'table' to replace ';'
    [b,c]=strtok(a);
    [e,f]=strtok(c);
    data(j,1)=str2double(b);
    data(j,2)=str2double(e);
    data(j,3)=str2double(f);
end
4事實上，我們的數(shù)據(jù)量比較大，可以一次打開多個文件
clc
close all
clear all
t1=27;
t2=30;
for t=t1:1:t2
f=strcat(num2str(t),'.cvs');
g='C:Documents and SettingsAdministratorMy Documents4.13';
fid=fopen(strcat(g,f),'r+');
t1,t2可以自己指定。因為實驗數(shù)據(jù)是要測60次取均值的，導師一直是這么要求的。如果讓你刪60個文件的說明的話，會發(fā)瘋的。
5上面的這些都是對不會使用dos的人寫的
如果你會使dos，十分簡單，寫一個批處理文件，在matlab的m文件的一開始執(zhí)行一下
       Frequency / GHz        SZmin(1),Zmin(1)/abs,linear
-----------------------------------------------------------------
                     5                           0.4557996
                 5.029                           0.4553181
                 5.058                           0.4547905
馬上全部變成
                     5                           0.4557996
                 5.029                           0.4553181
                 5.058                           0.4547905
現(xiàn)在我們的實驗室都是在用批處理的方式來寫的，因為數(shù)據(jù)太多了。。。。
先寫到這里了，希望對大家數(shù)據(jù)處理有幫助。