 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學(xué)術(shù)科研互動平臺 » 計算模擬區(qū) » 程序語言 » 其它 » matlab求助
| 5 | 1/1 | 返回列表 |
| 查看: 790 | 回復(fù): 4 | |||
heio金蟲 (小有名氣)
|
[求助]
matlab求助
|
|
本人打算用matlab處理數(shù)據(jù)中的趨勢項問題 但無法運行 跪求高手指點 程序如下 %滑動平均法消除趨勢項 clear clc close all hidden fni=input('輸入數(shù)據(jù)文件名稱','s'); fid=fopen(fni,'r') sf=fscanf(fid,'%f',1); %讀入采樣頻率值 l=fscanf(fid,'%d',1); %讀入滑動階次 m=fscanf(fid,'%d',1); %讀入平滑次數(shù) fno=fscanf(fid,'%d',1); %讀入輸出數(shù)據(jù)文件名 x=fscanf(fid,'%f',inf); %讀入時程數(shù)據(jù)存成列向量 status=fclose(fid); n=length(x); t=(0:1/sf : (n-1)/sf); b=ones(1,l); a=[b*x(1),x,b*x(n)]; b=a; for k=1:m for j=l+1:n-1 b(j)=mean(a(j-1:j+1)); end a=b; end y=x(1:n)-a(l+1:n+l); plot(t,x,':',t,y,t,a(l+1:n+l),'-.'); xlabel('時間(s)'); ylabel('位移 mm'); legend('輸入','輸出','趨勢'); grid on; fid=fopen(fno,'w'); for k=1:n fprintf(fid,'%f%f\n',t(k),y(k)); end status=fclose(fid); [ Last edited by heio on 2012-8-9 at 10:29 ] |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有7人回復(fù)
-
272求調(diào)劑
已經(jīng)有7人回復(fù)
-
352分 化工與材料
已經(jīng)有3人回復(fù)
-
0703化學(xué)
已經(jīng)有6人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿南師大0703化學(xué) 275求調(diào)劑
已經(jīng)有4人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有10人回復(fù)
csgt0
榮譽版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗: +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
- 注冊: 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計算模擬
heio
2樓2012-08-09 10:51:34
heio
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1363.9
- 散金: 41
- 帖子: 167
- 在線: 66.1小時
- 蟲號: 1334666
- 注冊: 2011-06-30
- 性別: GG
- 專業(yè): 光學(xué)
3樓2012-08-09 12:32:10
小新愛小白
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1365.5
- 紅花: 1
- 帖子: 80
- 在線: 51.8小時
- 蟲號: 1919956
- 注冊: 2012-08-01
- 性別: GG
- 專業(yè): 工程管理
|
Function [eval]=targetalloc (chrom) [m,n]=size(chrom); p=[87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;62 87 70 22 80 42 43 90 13 95 18 19 12 61 35;48 20 42 16 43 58 69 03 34 72 15 24 29 30 75]; w=[47 97 76 62 48 77 33 74 54 65 43 35 63 66 57]; for i=1:m for j=1:15 chrom(i,j)=p(chrom(i,j),j); end; end eval=chrom*w'; NIND=40; MAXGEN=50; GGAP=0.9; trace=zeros(MAXGEN,2); BaseV=crtbase(15,8); Chrom=crtbp(NIND,BaseV)+ones(NIND,15); gen=0; ObjV=taretalloc(Chrom); while gen SelCh=select('sus',Chrom,FitnV,GGAP); SelCh=recombin('xovsp',SelCh,0.7); f=rep([1;8],[1,15]); SelCh=mutbga(SelCh,f);SelCh=fix(SelCh); ObjVSel=taretalloc(SelCh); [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); gen=gen+1; trace(gen,1)=max(ObjV); trace(gen,2)=sum(ObjV)/length(ObjV); end [Y,I]=max(ObjV);Chrom(I,  ,Yplot(trace(:,1),'-.');hold on; plot(trace(:,2));grid Legend('解的變化','種群均值的變化') 這個程序有什么問題，怎么總是提示錯誤，也運行不了呢？或者請指教怎么在Matlab 7.7.0（R2008b）里邊具體步驟。 O(∩_∩)O謝謝 |
4樓2012-08-11 20:08:25
小新愛小白
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1365.5
- 紅花: 1
- 帖子: 80
- 在線: 51.8小時
- 蟲號: 1919956
- 注冊: 2012-08-01
- 性別: GG
- 專業(yè): 工程管理
|
Function [eval]=targetalloc (chrom) [m,n]=size(chrom); p=[87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;87 52 11 78 72 69 94 72 36 28 27 74 24 78 45;62 87 70 22 80 42 43 90 13 95 18 19 12 61 35;48 20 42 16 43 58 69 03 34 72 15 24 29 30 75]; w=[47 97 76 62 48 77 33 74 54 65 43 35 63 66 57]; for i=1:m for j=1:15 chrom(i,j)=p(chrom(i,j),j); end; end eval=chrom*w'; NIND=40; MAXGEN=50; GGAP=0.9; trace=zeros(MAXGEN,2); BaseV=crtbase(15,8); Chrom=crtbp(NIND,BaseV)+ones(NIND,15); gen=0; ObjV=taretalloc(Chrom); while gen SelCh=select('sus',Chrom,FitnV,GGAP); SelCh=recombin('xovsp',SelCh,0.7); f=rep([1;8],[1,15]); SelCh=mutbga(SelCh,f);SelCh=fix(SelCh); ObjVSel=taretalloc(SelCh); [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); gen=gen+1; trace(gen,1)=max(ObjV); trace(gen,2)=sum(ObjV)/length(ObjV); end [Y,I]=max(ObjV);Chrom(I, ,Yplot(trace(:,1),'-.');hold on; plot(trace(:,2));grid Legend('解的變化','種群均值的變化') 這個程序有什么問題，怎么總是提示錯誤，也運行不了呢？或者請指教怎么在Matlab 7.7.0（R2008b）里邊具體步驟。 O(∩_∩)O謝謝 |
5樓2012-08-11 20:08:42
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 heio 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 零八# 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 292求調(diào)劑 +13 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 14/700 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|