 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 789 | 回復(fù): 4 | |||
heio金蟲 (小有名氣)
|
[求助]
matlab求助
|
|
本人打算用matlab處理數(shù)據(jù)中的趨勢項(xiàng)問題 但無法運(yùn)行 跪求高手指點(diǎn) 程序如下 %滑動(dòng)平均法消除趨勢項(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); %讀入滑動(dòng)階次 m=fscanf(fid,'%d',1); %讀入平滑次數(shù) fno=fscanf(fid,'%d',1); %讀入輸出數(shù)據(jù)文件名 x=fscanf(fid,'%f',inf); %讀入時(shí)程數(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('時(shí)間(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ù)此樓» 猜你喜歡
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有13人回復(fù)
-
289求調(diào)劑
已經(jīng)有11人回復(fù)
-
網(wǎng)絡(luò)空間安全0839招調(diào)劑
已經(jīng)有5人回復(fù)
-
材料277求調(diào)劑
已經(jīng)有4人回復(fù)
-
274求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑
已經(jīng)有3人回復(fù)
-
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有8人回復(fù)
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號(hào): 449046
- 注冊(cè): 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
heio
2樓2012-08-09 10:51:34
heio
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1363.9
- 散金: 41
- 帖子: 167
- 在線: 66.1小時(shí)
- 蟲號(hào): 1334666
- 注冊(cè): 2011-06-30
- 性別: GG
- 專業(yè): 光學(xué)
3樓2012-08-09 12:32:10
小新愛小白
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1365.5
- 紅花: 1
- 帖子: 80
- 在線: 51.8小時(shí)
- 蟲號(hào): 1919956
- 注冊(cè): 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('解的變化','種群均值的變化') 這個(gè)程序有什么問題，怎么總是提示錯(cuò)誤，也運(yùn)行不了呢？或者請(qǐng)指教怎么在Matlab 7.7.0（R2008b）里邊具體步驟。 O(∩_∩)O謝謝 |
4樓2012-08-11 20:08:25
小新愛小白
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1365.5
- 紅花: 1
- 帖子: 80
- 在線: 51.8小時(shí)
- 蟲號(hào): 1919956
- 注冊(cè): 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('解的變化','種群均值的變化') 這個(gè)程序有什么問題，怎么總是提示錯(cuò)誤，也運(yùn)行不了呢？或者請(qǐng)指教怎么在Matlab 7.7.0（R2008b）里邊具體步驟。 O(∩_∩)O謝謝 |
5樓2012-08-11 20:08:42
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 heio 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +11 | 碩星赴 2026-03-23 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +3 | Mr. Z 2026-03-25 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +3 | 鈕恩雪 2026-03-25 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|