小妹最近才開始寫程式,碰到好多問題,煩請版上大大幫忙指點一下~^~,以下為程式碼:
clc;
clear all;
fidx = fopen('kk3.txt');                   
test1 = fscanf(fidx,'%g',[15 inf]);
fclose(fidx);

tot=size(test1);
tot=tot(1,2);

cd=0.98;

for i=1:tot
    if test1(9,i) >2 
        a1=1;
    else
        a1=0;
    end
    if test1(10,i) >2 
        a2=1;
    else
        a2=0;
    end
    if test1(11,i) >2 
        a3=1;
    else
        a3=0;
    end
    if test1(12,i) >2 
        a4=1;
    else
        a4=0;
    end
    if test1(13,i) >2 
        a5=1;
    else
        a5=0;
    end
    
    a(i)=3.141596*(a1*(0.038/2)^2+a2*(0.065/2)^2+a3*(0.1/2)^2+a4*(0.1/2)^2+a5*(0.125/2)^2);
    vn(i)=1.2931*273/(273+test1(2,i))*test1(1,i)/760;
    Q(i)=60*cd*a(i)*sqrt(2*9.81*test1(8,i)*1/vn(i));
    if test1(14,i)<2
            Q(i)=0;
    end
    pn(i)=test1(7,i);
end

plot(Q,pn,'+')
axis([5,15,-5,15])
xlabel('Q(m^3/min)')
ylabel('靜壓(pa)')

===========================================================================================
程式主要是將kk3.txt檔的資料做判斷與計算,然在把圖畫出來,程式邏輯前半部是判斷開與關(txt檔中的9~13行)
a(i)是由前半部程式判斷所得(txt檔中的9~13行),
vn(i)是由(txt檔中的第1,2行所算得)
Q(i)是由a(i), vn(i), 與txt檔中的第8行所得
pn(i)是txt檔中的第7行資料
kk3.txt檔裡的資料共15行,其中第14行是判斷資料是否為我們所要的,所以資料是一段一段的,一段要一段不要,
其實程式就是除了要判斷要與不要外,還要把要的這一段的Q(i)與pn(i)做平均,在劃出來,
小妹只完成把資料都都抓起來,不要的資料當成0,
資料沒有平均,都推在一起@@  (X軸是Q, Y軸是pn)
  
希望趕上大大,能指導一下小妹,該怎樣處理這些資料@@...把每一段資料做平均,然後劃出來^^...麻煩大家了
這是kk3.txt檔:
[url=http://0rz.tw/164Fm]http://0rz.tw/164Fm[/url]