 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 799 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
heio金蟲 (小有名氣)
|
[求助]
matlab求助
|
||
|
本人打算用matlab處理數(shù)據(jù)中的趨勢項(xiàng)問題 但無法運(yùn)行 跪求高手指點(diǎn) 程序如下 %滑動平均法消除趨勢項(xiàng) 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ù)此樓» 猜你喜歡
0805求調(diào)劑
已經(jīng)有5人回復(fù)
-
342求調(diào)劑
已經(jīng)有7人回復(fù)
-
調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有3人回復(fù)
-
273求調(diào)劑
已經(jīng)有4人回復(fù)
-
339求調(diào)劑
已經(jīng)有3人回復(fù)
-
267求調(diào)劑
已經(jīng)有7人回復(fù)
-
070300化學(xué)279求調(diào)劑
已經(jīng)有11人回復(fù)
-
318求調(diào)劑
已經(jīng)有3人回復(fù)
-
309分085801求調(diào)劑
已經(jīng)有4人回復(fù)
heio
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1363.9
- 散金: 41
- 帖子: 167
- 在線: 66.1小時
- 蟲號: 1334666
- 注冊: 2011-06-30
- 性別: GG
- 專業(yè): 光學(xué)
送鮮花一朵 
3樓2012-08-09 12:32:10
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
- 注冊: 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
2樓2012-08-09 10:51:34
小新愛小白
金蟲 (小有名氣)
- 應(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('解的變化','種群均值的變化') 這個程序有什么問題，怎么總是提示錯誤，也運(yùn)行不了呢？或者請指教怎么在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('解的變化','種群均值的變化') 這個程序有什么問題，怎么總是提示錯誤，也運(yùn)行不了呢？或者請指教怎么在Matlab 7.7.0（R2008b）里邊具體步驟。 O(∩_∩)O謝謝 |
5樓2012-08-11 20:08:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +3 | 雙馬尾痞老板2 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 考研材料工程351分調(diào)劑 +4 | 整個好的 2026-03-31 | 4/200 |
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +12 | 蛋黃咸肉粽 2026-03-26 | 13/650 |
|
|
|
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|