[求助]新手求助@@
小妹最近才開始寫程式,碰到好多問題,煩請版上大大幫忙指點一下~^~,以下為程式碼: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)
[url=http://www.]http://www.[/url]
這是放大圖
[url=http://www.]http://www.[/url]
希望趕上大大,能指導一下小妹,該怎樣處理這些資料@@...把每一段資料做平均,然後劃出來^^...麻煩大家了
這是kk3.txt檔:[url=http://www.]http://www.[/url]