 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學(xué)術(shù)科研互動平臺 » 計算模擬區(qū) » 程序語言 » 其它 » Matlab
| 6 | 1/1 | 返回列表 |
| 查看: 669 | 回復(fù): 5 | |||
[求助]
Matlab
|
|
如下表，希望通過循環(huán)數(shù)和步數(shù)的控制，通過設(shè)置電壓(第二列)的誤差，實現(xiàn)數(shù)據(jù)的挑選。如果不加循環(huán)和步數(shù)項，下面code可以實現(xiàn)單步工作，但是，增加了循環(huán)(n)和步數(shù)(m)后，就不工作了。 另外，還有一些小問題（見code標注），如：如果小數(shù)點后的位數(shù)很多時，似乎zx(find(zx==0))=[]; zy(find(zy==0))=[]; NewArray=[zx,zy]存在問題。急求高手指導(dǎo)！謝謝！ y = data(:,2); %對應(yīng)y值 u = data(:,3); %對應(yīng)步數(shù) v = data(:,4); %對應(yīng)循環(huán)數(shù) zx =zeros(length(x),1); zy =zeros(length(y),1); h = 1; zx(1) = x(1);%保留第一步的數(shù)據(jù)（如何實現(xiàn)保留其他步的第一數(shù)據(jù)？） zy(1) = y(1); a=[10,1,0.1]; B=[0,(-0.1)*a(1),(-1)*a(1),0.1*a(1),1*a(1),(-0.1)*a(2),(-1)*a(2),0.1*a(2),1*a(2),(-0.1)*a(3),(-1)*a(3),0.1*a(3),1*a(3)];%B內(nèi)部數(shù)的個數(shù)與步數(shù)對應(yīng) A=B*0.02; for k = 2:length(x) for n=1:max(v); for m=2:max(u); if abs(y(k)-y(k-1))>=abs(A(m)); %現(xiàn)在好像A(2)等同時影響其他步的誤差？） h = h+1; zx(h) = x(k-1); zy(h) = y(k-1); end end end end zx(find(zx==0))=[]; zy(find(zy==0))=[]; NewArray=[zx,zy]; disp(NewArray); plot(zx,zy,zx,zy,'o'); 時間 電壓 步數(shù) 循環(huán)數(shù) 5.008124936 -6.0395E-05 1 1 10.00032039 0.499913335 2 1 30.01483456 0.500076175 2 1 … … … … … … … … 310.0146066 0.500076175 2 1 310.0304182 5.000491142 3 1 … … … … … … … … 610.0291657 5.000491142 3 1 610.0294107 -0.49970841 4 1 630.0439261 -0.49987125 4 1 … … … … … … … … 890.2321767 -0.49987125 4 1 910.0436822 -0.49987125 4 1 910.059488 0.499913335 2 2 930.0738865 0.500076175 2 2 … … … … … … … … 1170.247626 0.499913335 2 2 1190.262159 0.500076175 2 2 1210.073952 0.499913335 2 2 1210.074271 5.000654221 3 2 … … … … … … … … 1490.276902 5.000491142 3 2 1510.088659 5.000654221 3 2 1510.088878 -0.49954554 4 2 1530.103367 -0.49987125 4 2 … … … … … … … … 1770.277052 -0.49987125 4 2 1790.291507 -0.49987125 4 2 1810.103251 -0.49987125 4 2 1810.119112 0.499750465 2 3 1830.133624 0.500076175 2 3 … … … … 2090.321718 0.500076175 2 3 2110.133399 0.500076175 2 3 … … … … |
回復(fù)此樓» 猜你喜歡
316求調(diào)劑
已經(jīng)有6人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
-
一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！
已經(jīng)有3人回復(fù)
-
310求調(diào)劑
已經(jīng)有3人回復(fù)
-
機械專碩299求調(diào)劑至材料
已經(jīng)有4人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有12人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
風泊
鐵桿木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 8548.1
- 帖子: 62
- 在線: 394.5小時
- 蟲號: 1016311
- 注冊: 2010-05-11
- 專業(yè): 金屬結(jié)構(gòu)材料
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
不確定你的誤差數(shù)組B如何和步數(shù)對應(yīng)，所以下面代碼只考慮了單一誤差，請自行測試 load data; data123=data(:,1:3); %取原始數(shù)據(jù)中的1、2、3列 Data4col_uni=unique(data(:,4)); NewArray=[]; for i=1:length(Data4col_uni) datafind=find(data(:,4)==i); Temp=data123(datafind,  ; %把循環(huán)i的數(shù)據(jù)取出來data12=Temp(:,1:2); %取Temp中的1、2列數(shù)據(jù) Data3col_uni=unique(Temp(:,3)); for j=1:length(Data3col_uni) datafind_1=find(Temp(:,3)==j); dataSingleStep=data12(datafind_1, ;%把第i循環(huán)，第j步時的數(shù)據(jù)取出來if(length(dataSingleStep(:,1))==1) NewArray=[NewArray;dataSingleStep]; else x = dataSingleStep(:,1); y = dataSingleStep(:,2); zx =zeros(length(x),1); zy =zeros(length(y),1); h = 1; zx(1) = x(1); zy(1) = y(1); for k = 2:length(x) if abs(y(k)-y(h)) >= 0.0002 zx(h) = x(k); zy(h) = y(k); h = h+1; end end zx(zx==0)=[]; zy(zy==0)=[]; NewArrayTemp=[zx,zy]; NewArray=[NewArray;NewArrayTemp]; end end end subplot(2,2,1); plot(data(:,1),data(:,2),data(:,1),data(:,2),'o'); subplot(2,2,2); plot(NewArray(:,1),NewArray(:,2),NewArray(:,1),NewArray(:,2),'o'); PS：請給出足夠的數(shù)據(jù)以便進行測試，謝謝。 |
2樓2012-12-10 19:43:46
風泊
鐵桿木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 8548.1
- 帖子: 62
- 在線: 394.5小時
- 蟲號: 1016311
- 注冊: 2010-05-11
- 專業(yè): 金屬結(jié)構(gòu)材料
3樓2012-12-10 19:52:58
4樓2012-12-10 20:40:50
realghost828
金蟲 (小有名氣)
- 應(yīng)助: 18 (小學(xué)生)
- 金幣: 1191.2
- 帖子: 76
- 在線: 57.5小時
- 蟲號: 532963
- 注冊: 2008-03-25
- 性別: GG
- 專業(yè): 運籌學(xué)
5樓2012-12-10 21:46:56
風泊
鐵桿木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 8548.1
- 帖子: 62
- 在線: 394.5小時
- 蟲號: 1016311
- 注冊: 2010-05-11
- 專業(yè): 金屬結(jié)構(gòu)材料
6樓2012-12-10 23:01:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 btliu 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 279分求調(diào)劑 一志愿211 +11 | chaojifeixia 2026-03-19 | 12/600 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|