版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

上一頁(yè)

下一頁(yè)

shengjunjie

木蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3613
散金: 208
紅花: 2
帖子: 320
在線: 128.7小時(shí)
蟲(chóng)號(hào): 988929
注冊(cè): 2010-04-04
性別: GG
專業(yè): 橡膠及彈性體

引用回帖:

8樓: Originally posted by dbb627 at 2012-04-17 10:16:54:
把你畫這個(gè)圖的代碼，貼一下，由于你的條件沒(méi)有給f，我畫不出這樣的圖
得到數(shù)據(jù)后，我可以給你做下fourier分析

頻率那一塊就是含在sin(ωt)里面，其中ω=2*pi*f，頻率f的范圍是0-1000Hz

謝謝

贊一下

回復(fù)此樓

好好學(xué)習(xí)，多運(yùn)動(dòng)，多發(fā)Paper...

11樓2012-04-17 10:27:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲(chóng)號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
shengjunjie: 金幣+5, ★★★很有幫助 2012-04-17 15:05:31
xiegangmai: 金幣+2, 謝謝應(yīng)助！ 2012-05-13 23:04:50

odefun=@(t,x)[x(2);3*x(2).^2/(2*x(1)*(1+0.08*x(1).^3))-1.342*10^7*(x(1).^4-x(1)+(61.98/6.71)*(0.185*sin(2*pi*335*t)).^2)./(1+0.08*x(1).^3)];
[t,x]=ode45(odefun,[0 0.06],[1 0]);% 自變量的范圍和初始值
figure(1);
plot(t,x(:,1))
Y=x(:,1)
pn=floor(log2(length(Y)));
N=2^pn;
S=Y(1:N);
Y1 = fft(S);%進(jìn)行N個(gè)點(diǎn)的fft變換
figure;
f1 = 1/(t(2)-t(1))*(0:N/2)/N;
Yabs=abs(Y1(1:N/2+1));
Ym=find(diff(sign(diff(Yabs)))<0)+1;%求幅值極大值的索引
Ymax=max(Yabs(Ym));
Yabs(Yabs Ym1=find(diff(sign(diff(Yabs)))<0)+1%求幅值極大值的索引
;plot(f1,Yabs(1:N/2+1))
hold on;plot(f1(Ym1),Yabs(Ym1),'r*');
axis([0 f1(Ym1(end)) 0 1.1*Ymax]);
xlabel('頻率 (Hz)')
ylabel('振幅|A| (m)')

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

12樓2012-04-17 10:43:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shengjunjie

木蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3613
散金: 208
紅花: 2
帖子: 320
在線: 128.7小時(shí)
蟲(chóng)號(hào): 988929
注冊(cè): 2010-04-04
性別: GG
專業(yè): 橡膠及彈性體

引用回帖:

12樓: Originally posted by dbb627 at 2012-04-17 10:43:35:
odefun=@(t,x);
=ode45(odefun,,);% 自變量的范圍和初始值
figure(1);
plot(t,x(:,1))
Y=x(:,1)
pn=floor(log2(length(Y)));
N=2^pn;
S=Y(1:N);
Y1 = fft(S);%進(jìn)行N個(gè)點(diǎn)的fft變換
figure;
f1 = 1/(t(2) ...

我在公式里面加了阻尼那個(gè)選項(xiàng)，和別人的公式一樣，但是感覺(jué)結(jié)果和他的差距很大�。ㄎ灰魄€一樣，但是振幅曲線差距很大），不知道什么原因
這是原圖

然后我根據(jù)你給我程序的結(jié)果是
odefun=@(t,x)[x(2);3*x(2).^2/(2*x(1)*(1+0.08*x(1).^3))-1.342*10^7*(x(1).^4-x(1)+(61.98/6.71)*(0.185*sin(2*pi*410*t)).^2)./(1+0.08*x(1).^3)-(5*10^3)*x(1).^3*x(2)./(25+2*x(1).^3)];
[t,x]=ode45(odefun,[0 0.06],[1 0]);% 自變量的范圍和初始值
figure(1);
plot(t,x(:,1))
Y=x(:,1)
pn=floor(log2(length(Y)));
N=2^pn;
S=Y(1:N);
Y1 = fft(S);%進(jìn)行N個(gè)點(diǎn)的fft變換
figure;
f1 = 1/(t(2)-t(1))*(0:N/2)/N;
Yabs=abs(Y1(1:N/2+1));
Ym=find(diff(sign(diff(Yabs)))<0)+1;%求幅值極大值的索引
Ymax=max(Yabs(Ym));
Yabs(Yabs Ym1=find(diff(sign(diff(Yabs)))<0)+1;%求幅值極大值的索引
plot(f1,Yabs(1:N/2+1))
hold on;
plot(f1(Ym1),Yabs(Ym1),'r*');
axis([0 f1(Ym1(end)) 0 1.1*Ymax]);
xlabel('頻率 (Hz)')
ylabel('振幅|A| (m)')
結(jié)果

贊一下

回復(fù)此樓

好好學(xué)習(xí)，多運(yùn)動(dòng)，多發(fā)Paper...

13樓2012-04-17 13:53:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲(chóng)號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

引用回帖:

13樓: Originally posted by shengjunjie at 2012-04-17 13:53:19:
我在公式里面加了阻尼那個(gè)選項(xiàng)，和別人的公式一樣，但是感覺(jué)結(jié)果和他的差距很大�。ㄎ灰魄€一樣，但是振幅曲線差距很大），不知道什么原因
這是原圖
71/2c/988929_1334641795_185.jpg
然后我根據(jù)你給我程序的 ...

我不知道，你對(duì)這些一點(diǎn)都不了解嗎？我只是了解一點(diǎn)這方面的理論，他后面那個(gè)圖應(yīng)該是有濾波處理吧，怎么那么平啊。具體的還是靠你自己查資料吧

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

14樓2012-04-17 14:24:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲(chóng)號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★ ★
shengjunjie: 金幣+2, ★有幫助 2012-04-17 15:05:45
xiegangmai: 金幣+1, 謝謝參與！ 2012-05-13 23:05:05

引用回帖:

13樓: Originally posted by shengjunjie at 2012-04-17 13:53:19:
我在公式里面加了阻尼那個(gè)選項(xiàng)，和別人的公式一樣，但是感覺(jué)結(jié)果和他的差距很大啊（位移曲線一樣，但是振幅曲線差距很大），不知道什么原因
這是原圖
71/2c/988929_1334641795_185.jpg
然后我根據(jù)你給我程序的 ...

[t,x]=ode45(odefun,[0 0.06],[1 0]);% 自變量的范圍和初始值的
t這里面的時(shí)間步長(zhǎng)可能不相等，對(duì)于fourier變換有影響，你重新調(diào)整為等步長(zhǎng)，再進(jìn)行fourier分析試試。

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

15樓2012-04-17 14:37:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shengjunjie

木蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3613
散金: 208
紅花: 2
帖子: 320
在線: 128.7小時(shí)
蟲(chóng)號(hào): 988929
注冊(cè): 2010-04-04
性別: GG
專業(yè): 橡膠及彈性體

引用回帖:

14樓: Originally posted by dbb627 at 2012-04-17 14:24:49:
我不知道，你對(duì)這些一點(diǎn)都不了解嗎？我只是了解一點(diǎn)這方面的理論，他后面那個(gè)圖應(yīng)該是有濾波處理吧，怎么那么平啊。具體的還是靠你自己查資料吧

我也是剛接觸這些非線性振動(dòng)方面的問(wèn)題，還是非常感謝你。我再查查資料，你辛苦了

贊一下

回復(fù)此樓

好好學(xué)習(xí)，多運(yùn)動(dòng)，多發(fā)Paper...

16樓2012-04-17 15:03:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shengjunjie

木蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3613
散金: 208
紅花: 2
帖子: 320
在線: 128.7小時(shí)
蟲(chóng)號(hào): 988929
注冊(cè): 2010-04-04
性別: GG
專業(yè): 橡膠及彈性體

引用回帖:

15樓: Originally posted by dbb627 at 2012-04-17 14:37:21:
=ode45(odefun,,);% 自變量的范圍和初始值的
t這里面的時(shí)間步長(zhǎng)可能不相等，對(duì)于fourier變換有影響，你重新調(diào)整為等步長(zhǎng)，再進(jìn)行fourier分析試試。

謝謝，我再看看，麻煩你拉

贊一下

回復(fù)此樓

好好學(xué)習(xí)，多運(yùn)動(dòng)，多發(fā)Paper...

17樓2012-04-17 15:04:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sugerbee

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.9
紅花: 2
帖子: 22
在線: 8.4小時(shí)
蟲(chóng)號(hào): 1772811
注冊(cè): 2012-04-23
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

【答案】應(yīng)助回帖

★ ★
xiegangmai: 金幣+1, 謝謝參與！ 2012-05-13 23:05:15
shengjunjie: 金幣+1, ★有幫助 2012-05-14 11:08:51

你的頻譜圖應(yīng)該去掉前面部位定的部分才對(duì)。因?yàn)槟愕姆l圖里有一部分是不穩(wěn)定的部分。所以在時(shí)域圖一致的情況下才出現(xiàn)不正確的幅頻圖。希望對(duì)你有用

贊一下

回復(fù)此樓

18樓2012-05-13 11:25:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sugerbee

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.9
紅花: 2
帖子: 22
在線: 8.4小時(shí)
蟲(chóng)號(hào): 1772811
注冊(cè): 2012-04-23
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

【答案】應(yīng)助回帖

★ ★ ★
xiegangmai: 金幣+1, 謝謝參與！ 2012-05-14 23:16:02
shengjunjie: 金幣+2, ★有幫助 2012-05-15 08:52:33

昨天忘了說(shuō)，論文里的那個(gè)幅頻圖不是傅里葉變換得到的。而是通過(guò)筆算得到幅值和頻率的函數(shù)之后作圖的。建議你先看看機(jī)械振動(dòng)，先給線性的大致了解一下在做。這樣方向比較明確。

贊一下

回復(fù)此樓

19樓2012-05-14 17:02:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shengjunjie

木蟲(chóng) (正式寫手)

引用回帖:

19樓: Originally posted by sugerbee at 2012-05-14 17:02:10:
昨天忘了說(shuō)，論文里的那個(gè)幅頻圖不是傅里葉變換得到的。而是通過(guò)筆算得到幅值和頻率的函數(shù)之后作圖的。建議你先看看機(jī)械振動(dòng)，先給線性的大致了解一下在做。這樣方向比較明確。

20樓2012-05-15 08:52:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shengjunjie 的主題更新

返回列表

上一頁(yè)

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +5	簡(jiǎn)單點(diǎn)0 2026-03-17	5/250	2026-03-18 08:13 by zhukairuo
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	5/250	2026-03-17 22:01 by 睿08
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂(lè) 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky