реклама на сайте
подробности

 
 
> Перенос спектра в Matlab
deneb13
сообщение Mar 7 2011, 18:03
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 32
Регистрация: 16-02-11
Пользователь №: 63 022



Здравствуйте! Возникла проблема. Реализую перенос спектра sin на 0-ю частоту. Для этого умножаю на экспоненту. Вопрос, почему из за переноса изменяются амплитудные соотношения?Совсем немного,но меняются!ведь по идее должны изменяться только частоты гармоник и все. Код привожу ниже
Fs=10*10^3;
f1=2*10^3;

t=0:1/Fs:1;
signal=1*sin(2*pi*f1*t);
y1=fft(signal);

figure(1)
plot(abs(y1)),grid


s1=signal.*exp(j*2*pi*f1*t);
figure(2)
y2=fft(s1);
plot(abs(y2)),grid
Спасибо!
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
deneb13
сообщение Mar 8 2011, 17:04
Сообщение #2


Участник
*

Группа: Участник
Сообщений: 32
Регистрация: 16-02-11
Пользователь №: 63 022



Спасибо огромное!действительно помогло!

Если можно, ответьте пожалуйста еще на один вопрос. Наверно опять проблема в погрешности fft, но тогда не понимаю, почему она возникает. С этим же сигналом, после фильтрации, выполняю обратный перенос на первоначальную частоту. Амплитуда опять меняется. Длины векторов остаются прежними. Что опять было упущено? Код прилагаю

Fs=10*10^3;
f1=2*10^3;
f2=4*10^3;
t=0:1/Fs:1-1/Fs;

signal=1*sin(2*pi*f1*t)+2*sin(2*pi*f2*t);
y1=fft(signal);

figure(1)
plot(abs(y1)),grid
title('Спектр исходного сигнала');


s1=signal.*exp(-j*2*pi*f1*t);
figure(2)
y2=fft(s1);
plot(abs(y2)),grid
title('Спектр перенесенного сигнала');

B1=[0.000074867323718623291
0.000029843947016906191
0.000035649251137170525
0.000042091749772785898
0.000049207651685282051
0.000057025721074006166
0.000065579982606295609
0.000074894688600103538
0.000084998464592263249
0.000095906944599545826
0.00010764057377583412
0.00012020685903585766
0.00013362042900918558
0.00014788529717745042
0.00016301296552518509
0.00017900206017028612
0.00019585384356087681
0.00021355033345440522
0.00023207543321342132
0.00025139354695834167
0.00027147928036987378
0.00029229420683612908
0.00031381075827335641
0.00033597892536717954
0.00035874915207741365
0.00038204131310657659
0.00040578976612778761
0.00042991522772108071
0.00045436483596751089
0.00047904333427598982
0.00050384929579228282
0.00052863752153019803
0.00055333850511602887
0.00057786619647086498
0.00060210611510683437
0.0006257905071241995
0.00064900905656298583
0.00067146426488381777
0.00069307362108383518
0.00071371214273047146
0.00073320062701199738
0.00075142572778662492
0.00076821265544567776
0.00078343209641189076
0.00079691819245689447
0.0008085381039427609
0.00081813426609220782
0.00082557461176558104
0.00083070708850672937
0.00083340327735858914
0.00083352260465997428
0.00083095026276507012
0.00082556391924235831
0.00081726029275528485
0.00080593050712979731
0.00079148599522964802
0.0007738417585242208
0.00075293820561771814
0.00072871712417941971
0.00070113590200281811
0.00067015698634955963
0.00063576966644221988
0.00059798041723357963
0.00055680986912850182
0.00051228299703510137
0.00046444844677513198
0.00041338731501109667
0.00035918843129745774
0.00030194252247403786
0.00024177341975123654
0.00017885598182614156
0.000113306600823241
0.000045333361450898598
-0.000024869351949700824
-0.000097103599499594588
-0.00017111187096671693
-0.00024666699809691337
-0.00032348451646964198
-0.00040129972048327968
-0.00047980547332992493
-0.00055870331403961237
-0.00063766033643191202
-0.00071635231963890027
-0.00079442819345953986
-0.00087154701127434156
-0.00094734204318419027
-0.0010214560517119006
-0.0010935128252451962
-0.0011631496035010649
-0.0012299924905660802
-0.0012936801049696345
-0.0013538395168376132
-0.0014101167323676916
-0.0014621561022228591
-0.001509623166065856
-0.0015521807461206609
-0.0015895143416746335
-0.0016213205887912342
-0.0016473241868954092
-0.0016672543865457272
-0.0016808686184831172
-0.0016879494710142569
-0.0016883117045997894
-0.0016817789958335427
-0.0016682176295261186
-0.0016475294459246312
-0.0016196272022528976
-0.0015844757988635676
-0.0015420758109769049
-0.0014924412224791917
-0.0014356555425266994
-0.0013718032498286298
-0.0013010394487680812
-0.0012235320329592899
-0.0011395078131137791
-0.0010492133697008541
-0.00095294911593929671
-0.00085103874162868699
-0.00074385900503890368
-0.00063180747056448285
-0.00051533159990874123
-0.00039489798082718256
-0.00027102046219501839
-0.00014423261171928799
-0.000015107192143127688
0.00011576835014704229
0.00024776794565646359
0.00038025412485131652
0.00051255881681363938
0.00064401087734885208
0.00077391268795518411
0.00090156324510963168
0.0010262451175025195
0.0011472502297884227
0.0012638570940360445
0.0013753543177702373
0.0014810315797351504
0.0015801993917129518
0.0016721654941180143
0.0017562700070420914
0.0018318714495729322
0.0018983499480612295
0.001955115604617092
0.0020016210874932226
0.0020373357023431358
0.0020617968844452677
0.0020745577194886465
0.0020752427668374626
0.0020635077831940296
0.0020390777645383951
0.0020017215556572092
0.0019512785593668654
0.0018876379184231378
0.00181076522969225
0.0017206794094045876
0.0016174801503727888
0.0015013219439372187
0.0013724410035904861
0.0012311356991687385
0.0010777854567875196
0.00091282935024940249
0.00073678829856948342
0.00055024358978284697
0.00035385856813305495
0.00014835278897694663
-0.000065477663074763916
-0.00028678000255334582
-0.00051462874765795794
-0.00074805087988135809
-0.00098600485267714174
-0.0012274029075307119
-0.0014711000738689599
-0.001715917082149265
-0.0019606173227907416
-0.0022039393041064555
-0.0024445823154837826
-0.0026812209683040001
-0.0029125026476055651
-0.0031370686755146132
-0.0033535313249636353
-0.0035605194051896976
-0.0037566406312513341
-0.0039405274575157303
-0.004110810694194075
-0.004266150890035114
-0.0044052253371144896
-0.00452675145053082
-0.0046294747186861962
-0.0047121940643470245
-0.0047737470119906299
-0.004813040330222243
-0.0048290301594689439
-0.0048207490678348949
-0.0047872940617394782
-0.0047278489212654352
-0.0046416710564801374
-0.0045281149425436956
-0.0043866162362661155
-0.0042167179759658125
-0.0040180521133887058
-0.0037903617578293621
-0.0035334877881856596
-0.0032473894309709458
-0.0029321244684466911
-0.002587872236078986
-0.0022149177580068837
-0.0018136676421232414
-0.0013846306049196367
-0.00092844222474954496
-0.00044583971879492545
0.000062320900346972839
0.00059507946179006342
0.0011513620223204378
0.0017300003802191962
0.002329709842908615
0.0029491251376265581
0.0035867707448209128
0.0042410962523754371
0.0049104570303590027
0.0055931368835668918
0.0062873399938202052
0.0069912108562696761
0.0077028263936973936
0.0084202192531677243
0.0091413635116792211
0.0098642072937224021
0.010586655725801431
0.011306599068558625
0.012021904004335312
0.012730436769847217
0.013430055818563241
0.014118637311338231
0.014794062900710676
0.015454250025499331
0.016097139668424538
0.016720721187340245
0.017323024712336078
0.017902147059561155
0.018456238197103358
0.018983528040541858
0.019482317990833237
0.01995100336042771
0.020388059889961278
0.02079207270759072
0.021161720312999361
0.0214957976242044
0.021793207009755969
0.02205297610813545
0.022274244895276123
0.022456291164941175
0.022598506895452168
0.0227004288620027
0.022761715857729192
0.022782168419279283
0.022761715857729192
0.0227004288620027
0.022598506895452168
0.022456291164941175
0.022274244895276123
0.02205297610813545
0.021793207009755969
0.0214957976242044
0.021161720312999361
0.02079207270759072
0.020388059889961278
0.01995100336042771
0.019482317990833237
0.018983528040541858
0.018456238197103358
0.017902147059561155
0.017323024712336078
0.016720721187340245
0.016097139668424538
0.015454250025499331
0.014794062900710676
0.014118637311338231
0.013430055818563241
0.012730436769847217
0.012021904004335312
0.011306599068558625
0.010586655725801431
0.0098642072937224021
0.0091413635116792211
0.0084202192531677243
0.0077028263936973936
0.0069912108562696761
0.0062873399938202052
0.0055931368835668918
0.0049104570303590027
0.0042410962523754371
0.0035867707448209128
0.0029491251376265581
0.002329709842908615
0.0017300003802191962
0.0011513620223204378
0.00059507946179006342
0.000062320900346972839
-0.00044583971879492545
-0.00092844222474954496
-0.0013846306049196367
-0.0018136676421232414
-0.0022149177580068837
-0.002587872236078986
-0.0029321244684466911
-0.0032473894309709458
-0.0035334877881856596
-0.0037903617578293621
-0.0040180521133887058
-0.0042167179759658125
-0.0043866162362661155
-0.0045281149425436956
-0.0046416710564801374
-0.0047278489212654352
-0.0047872940617394782
-0.0048207490678348949
-0.0048290301594689439
-0.004813040330222243
-0.0047737470119906299
-0.0047121940643470245
-0.0046294747186861962
-0.00452675145053082
-0.0044052253371144896
-0.004266150890035114
-0.004110810694194075
-0.0039405274575157303
-0.0037566406312513341
-0.0035605194051896976
-0.0033535313249636353
-0.0031370686755146132
-0.0029125026476055651
-0.0026812209683040001
-0.0024445823154837826
-0.0022039393041064555
-0.0019606173227907416
-0.001715917082149265
-0.0014711000738689599
-0.0012274029075307119
-0.00098600485267714174
-0.00074805087988135809
-0.00051462874765795794
-0.00028678000255334582
-0.000065477663074763916
0.00014835278897694663
0.00035385856813305495
0.00055024358978284697
0.00073678829856948342
0.00091282935024940249
0.0010777854567875196
0.0012311356991687385
0.0013724410035904861
0.0015013219439372187
0.0016174801503727888
0.0017206794094045876
0.00181076522969225
0.0018876379184231378
0.0019512785593668654
0.0020017215556572092
0.0020390777645383951
0.0020635077831940296
0.0020752427668374626
0.0020745577194886465
0.0020617968844452677
0.0020373357023431358
0.0020016210874932226
0.001955115604617092
0.0018983499480612295
0.0018318714495729322
0.0017562700070420914
0.0016721654941180143
0.0015801993917129518
0.0014810315797351504
0.0013753543177702373
0.0012638570940360445
0.0011472502297884227
0.0010262451175025195
0.00090156324510963168
0.00077391268795518411
0.00064401087734885208
0.00051255881681363938
0.00038025412485131652
0.00024776794565646359
0.00011576835014704229
-0.000015107192143127688
-0.00014423261171928799
-0.00027102046219501839
-0.00039489798082718256
-0.00051533159990874123
-0.00063180747056448285
-0.00074385900503890368
-0.00085103874162868699
-0.00095294911593929671
-0.0010492133697008541
-0.0011395078131137791
-0.0012235320329592899
-0.0013010394487680812
-0.0013718032498286298
-0.0014356555425266994
-0.0014924412224791917
-0.0015420758109769049
-0.0015844757988635676
-0.0016196272022528976
-0.0016475294459246312
-0.0016682176295261186
-0.0016817789958335427
-0.0016883117045997894
-0.0016879494710142569
-0.0016808686184831172
-0.0016672543865457272
-0.0016473241868954092
-0.0016213205887912342
-0.0015895143416746335
-0.0015521807461206609
-0.001509623166065856
-0.0014621561022228591
-0.0014101167323676916
-0.0013538395168376132
-0.0012936801049696345
-0.0012299924905660802
-0.0011631496035010649
-0.0010935128252451962
-0.0010214560517119006
-0.00094734204318419027
-0.00087154701127434156
-0.00079442819345953986
-0.00071635231963890027
-0.00063766033643191202
-0.00055870331403961237
-0.00047980547332992493
-0.00040129972048327968
-0.00032348451646964198
-0.00024666699809691337
-0.00017111187096671693
-0.000097103599499594588
-0.000024869351949700824
0.000045333361450898598
0.000113306600823241
0.00017885598182614156
0.00024177341975123654
0.00030194252247403786
0.00035918843129745774
0.00041338731501109667
0.00046444844677513198
0.00051228299703510137
0.00055680986912850182
0.00059798041723357963
0.00063576966644221988
0.00067015698634955963
0.00070113590200281811
0.00072871712417941971
0.00075293820561771814
0.0007738417585242208
0.00079148599522964802
0.00080593050712979731
0.00081726029275528485
0.00082556391924235831
0.00083095026276507012
0.00083352260465997428
0.00083340327735858914
0.00083070708850672937
0.00082557461176558104
0.00081813426609220782
0.0008085381039427609
0.00079691819245689447
0.00078343209641189076
0.00076821265544567776
0.00075142572778662492
0.00073320062701199738
0.00071371214273047146
0.00069307362108383518
0.00067146426488381777
0.00064900905656298583
0.0006257905071241995
0.00060210611510683437
0.00057786619647086498
0.00055333850511602887
0.00052863752153019803
0.00050384929579228282
0.00047904333427598982
0.00045436483596751089
0.00042991522772108071
0.00040578976612778761
0.00038204131310657659
0.00035874915207741365
0.00033597892536717954
0.00031381075827335641
0.00029229420683612908
0.00027147928036987378
0.00025139354695834167
0.00023207543321342132
0.00021355033345440522
0.00019585384356087681
0.00017900206017028612
0.00016301296552518509
0.00014788529717745042
0.00013362042900918558
0.00012020685903585766
0.00010764057377583412
0.000095906944599545826
0.000084998464592263249
0.000074894688600103538
0.000065579982606295609
0.000057025721074006166
0.000049207651685282051
0.000042091749772785898
0.000035649251137170525
0.000029843947016906191
0.000074867323718623291];
B=B1';
z=filtfilt(B,1,s1);


y2(10:(length(y2)-10))=0;
figure(3)
plot(abs(y2)),grid
s_f = ifft(y2);
title('Спектр фильтрованного сигнала ОБНУЛЕНИЕМ');


s_obr=z.*exp(j*2*pi*f1*t);
spektr=fft(s_obr);
figure(4)
plot(abs(spektr)),grid
title('Спектр обратного переноса сигнала');

rez=signal-s_obr;
spektr_rez=fft(rez);


z0=fft(z);
figure(5)
plot(abs(z0)),grid
title('Спектр фильтрованного сигнала');

figure(6)
plot(abs(spektr_rez)),grid
title('КОМПЕНСАЦИЯ');

Спасибо!
P.S. Извините за возможно глупые вопросы-Matlab в стадии начального освоения blush.gif
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- deneb13   Перенос спектра в Matlab   Mar 7 2011, 18:03
- - Taradov Alexander   Разница зависит от частоты дискретизации и длинны ...   Mar 7 2011, 19:43
- - Самурай   Цитата(deneb13 @ Mar 7 2011, 21:03) Здрав...   Mar 7 2011, 20:34
- - Taradov Alexander   1. для создания фильтров в момент выполнения есть ...   Mar 8 2011, 17:20
- - deneb13   Исходная задача: сигнал есть сумма полезной состав...   Mar 8 2011, 17:24
|- - Taradov Alexander   Не думаю что так получится в жизни, не хватит точн...   Mar 8 2011, 17:33
- - deneb13   Это задание на курсовой. Нужно сделать именно комп...   Mar 8 2011, 17:38
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 8 2011, 20:38) Это з...   Mar 8 2011, 17:44
- - deneb13   Спасибо! Taradov Alexander, по вашему совету ...   Mar 8 2011, 18:16
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 8 2011, 21:16) Tarad...   Mar 8 2011, 18:25
- - deneb13   Помеха в дальнейшем так и будет sin-узкополосная.С...   Mar 8 2011, 18:47
- - Taradov Alexander   Хорошо, этот метод почти сработал для небольшого к...   Mar 8 2011, 19:03
|- - deneb13   Цитата(Taradov Alexander @ Mar 8 2011, 22...   Mar 9 2011, 11:11
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 14:11) Вопро...   Mar 9 2011, 11:17
- - deneb13   Например,в Тузове "Стат.теория приема сложных...   Mar 8 2011, 19:15
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 8 2011, 22:10) Напри...   Mar 8 2011, 19:31
- - deneb13   Спасибо!Вы мне очень помогли!А насчет книг...   Mar 8 2011, 19:32
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 8 2011, 22:32) Спаси...   Mar 8 2011, 19:37
- - deneb13   Спасибо!   Mar 8 2011, 19:47
- - 314   Извините, если не в тему, но такой метод фильтраци...   Mar 9 2011, 06:38
- - Taradov Alexander   Что-то не выходит у меня из головы этот метод, мож...   Mar 9 2011, 06:39
- - deneb13   Мои знания о фильтре Калмана достаточно поверхност...   Mar 9 2011, 08:33
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 11:33) Перен...   Mar 9 2011, 09:46
|- - deneb13   Цитата(Taradov Alexander @ Mar 9 2011, 12...   Mar 9 2011, 09:55
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 12:55) а как...   Mar 9 2011, 09:59
- - petrov   Узкополосную помеху можно адаптивным фильтром комп...   Mar 9 2011, 09:30
- - deneb13   Да,я знаю об адаптивной фильтрации. Просто по зада...   Mar 9 2011, 09:43
- - deneb13   я описала общий принцип метода компенсации.в своем...   Mar 9 2011, 09:49
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 12:49) Т.е. ...   Mar 9 2011, 09:55
- - deneb13   Спасибо,интересная для меня мысль!   Mar 9 2011, 10:01
- - Herz   Цитата(Taradov Alexander @ Mar 9 2011, 08...   Mar 9 2011, 10:15
- - deneb13   Да,блок всегда фиксированной длины.а если использо...   Mar 9 2011, 11:41
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 14:41) Да,бл...   Mar 9 2011, 11:53
- - petrov   deneb13, делайте с переносом в 0 и обратно, в нуле...   Mar 9 2011, 12:05
|- - Taradov Alexander   QUOTE (petrov @ Mar 9 2011, 15:05) deneb1...   Mar 9 2011, 12:13
|- - petrov   Цитата(Taradov Alexander @ Mar 9 2011, 15...   Mar 9 2011, 12:17
||- - Taradov Alexander   QUOTE (petrov @ Mar 9 2011, 15:17) Это бу...   Mar 9 2011, 12:22
|- - deneb13   Цитата(Taradov Alexander @ Mar 9 2011, 15...   Mar 9 2011, 12:45
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 15:45) именн...   Mar 9 2011, 12:52
- - deneb13   Буду благодарна,если расскажете о своих результата...   Mar 9 2011, 12:58
- - Taradov Alexander   Я решил попробовать для начала просто одну помеху...   Mar 9 2011, 18:43
|- - Самурай   Цитата(Taradov Alexander @ Mar 9 2011, 21...   Mar 9 2011, 19:30
|- - Taradov Alexander   QUOTE (Самурай @ Mar 9 2011, 22:30) Полос...   Mar 9 2011, 19:42
- - deneb13   Скажите, а как можно посчитать комплексный фильтр?...   Mar 9 2011, 20:48
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 23:48) Скажи...   Mar 9 2011, 20:58
- - deneb13   Для того,чтобы узнать насколько следует увеличиват...   Mar 9 2011, 20:59
|- - Taradov Alexander   QUOTE (deneb13 @ Mar 9 2011, 23:59) Для т...   Mar 9 2011, 21:03
|- - petrov   Цитата(Taradov Alexander @ Mar 10 2011, 00...   Mar 9 2011, 21:46
- - deneb13   Спасибо,попробую   Mar 9 2011, 21:05


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 30th July 2025 - 07:47
Рейтинг@Mail.ru


Страница сгенерированна за 0.0144 секунд с 7
ELECTRONIX ©2004-2016