| 5 | 1/1 | 返回列表 |
| 查看: 5094 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
用matlab進(jìn)行批量處理數(shù)據(jù)時(shí)遇到的問題 已有2人參與
|
|||
|
fluent計(jì)算結(jié)果輸出170個(gè)壓力分布數(shù)據(jù),現(xiàn)在需要把數(shù)據(jù)進(jìn)行批量處理求平均值。小女子新蟲一枚,金幣不多,望大神幫助。 下面是編的程序: %-------------------------------------------------------------------------% clear %導(dǎo)入原始數(shù)據(jù) cd('C:\Users\yangliu\Desktop\LES 未定義來流的湍流度\cp');%待處理的數(shù)據(jù)文件所在文件夾 Allname=struct2cell(dir); %得到上述文件夾下的所有文件名 %dir命令,可以得到路徑內(nèi)包括文件名在內(nèi)的文件信息,為struc數(shù)據(jù)結(jié)構(gòu)。 %-------------------------------------------------------------------------% %逐個(gè)讀入文件 [m,n]=size(Allname); y=zeros(300,1);%對(duì)矩陣y進(jìn)行賦值 for i=3:n %前兩個(gè)不是文件名 name=Allname{1,i};%第一行作為文件名用來提取數(shù)據(jù) hang=[5 304];%5行到304行 lie=[1 2];%1到2列 fid1=fopen('C:\Users\yangliu\Desktop\LES 未定義來流的湍流度\cp\name','r'); temp=textscan(fid1,'%f %f ',5,304);%讀取5-304行數(shù)據(jù) fclose(fid1); for k=1 hang(2)-hang(1)+1)for l=1 lie(2)-lie(1)+1)b=temp{l+lie(1)-1}{k+hang(1)-1}; A(k,l)=str2num(b); end % for l=1 lie(2)-lie(1)+1)end % for k=1 hang(2)-hang(1)+1)%將數(shù)據(jù)賦值到矩陣 B(i)=A(:,2); if mod(i,n)==1 fid2=fopen('C:\Users\yangliu\Desktop\LES 未定義來流的湍流度\cp\cpaverage.txt','wt'); D(:,1)=A(:,1); D(:,2)=y(:,1)/(n-2); ![]() ![]() ![]() ![]() ![]() ![]() for j=1:1:300 for m=1:1:2 if m==2 fprintf(fid,'%g\n',D(j,m)); else fprintf(fid,'%g\t',D(j,m)); end end end fclose(fid2);%將數(shù)據(jù)逐個(gè)寫入文件中 else i=i+1; y(:,1)=B(i)+y(:,1); end %for if mod(i,n)==1 end %for i=3:n 這是自己編的程序,讀取5-304行數(shù)據(jù),并對(duì)第二列進(jìn)行平均計(jì)算,提示錯(cuò)誤使用 textscan。文件標(biāo)識(shí)符無效。使用 fopen 生成有效的文件標(biāo)識(shí)符。后面附有一個(gè)待處理的文件。小女子matlab剛開始學(xué)習(xí),實(shí)在是不知道怎么改程序了,求各位大神幫助,不吝賜教啊。 ![]() ![]() |
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +3 |
新蟲 (小有名氣)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +3 |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑,不限專業(yè),有學(xué)上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
[考研] 311求調(diào)劑 +7 | 藍(lán)月亮亮 2026-03-30 | 7/350 |
|
|
[考研] 287求調(diào)劑 +14 | land xuxu 2026-03-26 | 14/700 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|