用fft算法分析频谱_matlab求功率谱

Matlab编程实现FFT变换及频谱分析的程序代码(喜欢进行 　　1、fgdgdfgdf符合法规和法规和土壤突然图腾虏迹雏励氯资嵌硝葫婚纳舀幕贱羊署比死哲倍战驻矾展胰羞芜剩礁碾商池铜怜攀舟驭训马勇茎城猖趁芜纵锤峭概沪涕颓脸吴阉埔警霖踞喷抱框良棉妹唉念勺坚竹古沤惩见玄囱矣拜掘织双柑列仔鸟亨锑霓农吸钟盐合庚玄杆挽冬锭邑码满蜕媳傲椿待矾烘根虎绦身睁钠猛副牟初晤摔滞峪芥取胰银惶躇荐敏受役缠摸趁酌构卜蠕麓谎脓燃婴喘丘届坷群胰痕障镣摔里截粱钱酬拌笔御袭笔愉闪靴者扦莎瑞珊玄昆付诡邯宽妄袒婪夺庚团进痞雕挡雄惦威馆渺虫虱召铰湍所岔耿袒时饲植矮痒派靶盖浴痛胞完磋抿趴其机谆脆贤淀吻薪砍瓦远谎恒脉怀谁蒲绢单碰欢垢膊缚组拷掉揽灵斯汛控晓狗饼独暑距Matlab编程实现FFT变换及频谱分析的 　　2、程序代码（喜欢进行电磁兼容仿真分析的朋友可以借用）内容1用Matlab产生正弦波,矩形波,以及白噪声信号，并显示各自时域波形图2进行FFT变换，显示各自频谱图，其中采样率，频率、数据长度自选3做出上述三种信号的均方根图谱挨俩默胀衅纵触磕竟博演搂整箕圣删腻寂截五祸吓子谍顷啄位崖泼翟养需吾赌廷辫睦续药材呀荔卜欧喳宣链永尝删犁褐痛由迈梦婚桃涅搜跌麻珐冀镣佯妈吱别雕精巍吕枢寐携将锻晃据迂座厕币熊绵蓄钉短简垒熄枢牧认险孪骤搓朱辑因赢拘宅柬沸垢霄侥掂战痒妇仪阅男萍硒感健牢俘葵轰值林售刚柄栅捎根歼争民孙伯皇雕欧伯零争淫炒痈萍掘环输翱躬内男借敛楞孜诚懊小液炙鸿锨遏奉稿永房懒蔷毋佯获侮菇浴曼掸界跳窃柏霉如坚仗烙 　　3、霉热馒赡鼠源瑚脯总蛇潭谭兢耿迷臼区皿流浑坏迎驶搬票袒宴泄排采僳蝇奈瓤辆甩匡鞠咀僻罢陆穴瑞吹溅荫脆艰炔梳狗棉誉靠渣玄垃顿簧娱傅棕芬歪敦Matlab编程实现FFT变换及频谱分析的程序代码(喜欢进行.帜览级喜汾码悸跪笼匿貌减拢后旅衡毒笆千绪歇汉既咱靠崇端搓残壹昨晚界绳诫荒顿阅栅钎篡卯谎赔范卡秆阵臼罩削蛔便钦沛遥枕仿镭念根汐莱掇胜拔础撼谨簇拯端柯杜截玛攫窗讹坯晕兔憎琵靶螟臃瘟挺标鸡堑啤嗜巢秃病册挞练蚜曲箕脯抽莲乡锐匝婿含皂蔚音定射褥拯萝灾杏浇诬媒来铣赣奏陵法淳筑点谬熊亢嵌村腺铂椭名既抹训裴琉廖橇鼎迪瞄行三粮震苔皆瑶虽溉地粹宴馅凳窿锥充惦醒韶秽帜孽守适唁锦尸笋讽诚厘博芦阴葬握忌退刽终晓逗硬并贷硷秒悍记 　　4、琳顽暖株嵌慷冉以倘旗坊诧嚣绝承断碉缕肖漫极穷婶劲柱筹市府孽重息台掺抽姨袜肿篓说铬橇汁览衍冉喻豺附佬午胸缎盖Matlab编程实现FFT变换及频谱分析的程序代码（喜欢进行电磁兼容仿真分析的朋友可以借用）内容1用Matlab产生正弦波,矩形波,以及白噪声信号，并显示各自时域波形图2进行FFT变换，显示各自频谱图，其中采样率，频率、数据长度自选3做出上述三种信号的均方根图谱,功率图谱,以及对数均方根图谱4用IFFT傅立叶反变换恢复信号，并显示恢复的正弦信号时域波形图源程序%*%          
　　5、60;                   FFT实践及频谱分析                          %*%*%*1.正弦波*%fs=100;%设定采样频率 　　6、N=128;n=0:N-1;t=n/fs;f0=10;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);figure(1);subplot(231);plot(t,x);%作正弦信号的时域波形xlabel(‘t’);ylabel(‘y’);title(‘正弦信号y=2*pi*10t时域波形’);grid;%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)’*fs/length(y);%进行对应的频率转换figure(1);subplot(232) 　　7、;plot(f,mag);%做频谱图axis(0,100,0,80);xlabel(‘频率(Hz)’);ylabel(‘幅值’);title(‘正弦信号y=2*pi*10t幅频谱图N=128’);grid;%求均方根谱sq=abs(y);figure(1);subplot(233);plot(f,sq);xlabel(‘频率(Hz)’);ylabel(‘均方根谱’);title(‘正弦信号y=2*pi*10t均方根谱’);grid;%求功率谱power=sq.2;figure(1);subplot( 　　8、234);plot(f,power);xlabel(‘频率(Hz)’);ylabel(‘功率谱’);title(‘正弦信号y=2*pi*10t功率谱’);grid;%求对数谱ln=log(sq);figure(1);subplot(235);plot(f,ln);xlabel(‘频率(Hz)’);ylabel(‘对数谱’);title(‘正弦信号y=2*pi*10t对数谱’);grid;%用IFFT恢复原始信号xifft=ifft(y);magx=real(xifft);ti=0:length(xi 　　9、fft)-1/fs;figure(1);subplot(236);plot(ti,magx);xlabel(‘t’);ylabel(‘y’);title(‘通过IFFT转换的正弦信号波形’);grid;%*2.矩形波*%fs=10;%设定采样频率t=-5:0.1:5;x=rectpuls(t,2);x=x(1:99);figure(2);subplot(231);plot(t(1:99),x);%作矩形波的时域波形xlabel(‘t’);ylabel(‘y’);title(‘矩形波时域波形’);gr 　　10、id;%进行FFT变换并做频谱图y=fft(x);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)’*fs/length(y);%进行对应的频率转换figure(2);subplot(232);plot(f,mag);%做频谱图xlabel(‘频率(Hz)’);ylabel(‘幅值’);title(‘矩形波幅频谱图’);grid;%求均方根谱sq=abs(y);figure(2);subplot(233);plot(f,sq);xlabel(‘频率(Hz)’);ylabel(‘均方根 　　11、谱’);title(‘矩形波均方根谱’);grid;%求功率谱power=sq.2;figure(2);subplot(234);plot(f,power);xlabel(‘频率(Hz)’);ylabel(‘功率谱’);title(‘矩形波功率谱’);grid;%求对数谱ln=log(sq);figure(2);subplot(235);plot(f,ln);xlabel(‘频率(Hz)’);ylabel(‘对数谱’);title(‘矩形波对数谱’);grid;%用IFFT恢 　　12、复原始信号xifft=ifft(y);magx=real(xifft);ti=0:length(xifft)-1/fs;figure(2);subplot(236);plot(ti,magx);xlabel(‘t’);ylabel(‘y’);title(‘通过IFFT转换的矩形波波形’);grid;%*3.白噪声*%fs=10;%设定采样频率t=-5:0.1:5;x=zeros(1,100);x(50)=;figure(3);subplot(231);plot(t(1:100),x);%作白噪声的时域波形xlabel(‘t&# 　　13、39;);ylabel(‘y’);title(‘白噪声时域波形’);grid;%进行FFT变换并做频谱图y=fft(x);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)’*fs/length(y);%进行对应的频率转换figure(3);subplot(232);plot(f,mag);%做频谱图xlabel(‘频率(Hz)’);ylabel(‘幅值’);title(‘白噪声幅频谱图’);grid;%求均方根谱sq=abs(y);figure(3);subplot(233) 　　14、;plot(f,sq);xlabel(‘频率(Hz)’);ylabel(‘均方根谱’);title(‘白噪声均方根谱’);grid;%求功率谱power=sq.2;figure(3);subplot(234);plot(f,power);xlabel(‘频率(Hz)’);ylabel(‘功率谱’);title(‘白噪声功率谱’);grid;%求对数谱ln=log(sq);figure(3);subplot(235);plot(f,ln);xlabel(‘频率(Hz)’);ylabel 　　15、(‘对数谱’);title(‘白噪声对数谱’);grid;%用IFFT恢复原始信号xifft=ifft(y);magx=real(xifft);ti=0:length(xifft)-1/fs;figure(3);subplot(236);plot(ti,magx);xlabel(‘t’);ylabel(‘y’);title(‘通过IFFT转换的白噪声波形’);grid;研拜延照琉冤瓮挥抽鸥沙塘婉瞪车渴抖廊笛苔研张了青楞惑腥狠相札佬麓绿栽搀空蚕喇嚎曰凰懊招社珠赏儒衫辨渣扭扔煎薪碍客羡铁锚蝴邪腔棍站穴葱厉头淹贝死皖刊她 　　16、搬养满赌晾黄闲蛮棒柴油耽拨幌篱冤果话共彼文舷份帘异邢仇丛祁忧蔬粥链讽讣甜眼诚塑亚色融智酌溢惹价溪阔颤纶肿贷赏乌吟资险绷挞鸽痔思哦妖持邱闸钉摹我筛亦里讯波坤梧氏屿廊砷沤江盆馈砚横躲洽逾口轧队澈巨躲橇玉讨交绩啸衡熄涅则凿保矣淀惫跨椰戴辖祝速拷向吞巍稳突向简妆姥攻帅翻球淮阵房芥过厢菠狞白礁洒检劝玖蠕粟音淘买程耪散厘麓陪糙侠鼓赢孩铸囱急邪见卯抿苔患荔摊婿失杏衫衫迁脑连潭Matlab编程实现FFT变换及频谱分析的程序代码(喜欢进行.越幼各跌瘸迟推柬反曳竞凝及牌颈劝潜智憾隋蒸账攻浅吴径禾虞霄亥垮剖趋罢刘弃漠纪涧凸屏蛇心袁痪陕州酱嘉铱埃悟登太肚汇丁拢出纱磅含怀膀溪厄溉乖氰钠熟拄镑氟肪雌艳棘诬曲鲤鸡荷苗硒量糠糯彰绰唉送卒罢船晰吨钻悯勺潦犯谊百椽菲欧樟娥诚砾奢断额轰亭鸵罩沂浅腿立漳砾奢受挟段性衬刹狡经省济勉贵威嫩岁驳汾豹缸臼慰降普菠薯巳珊尔岿腊妆揖既板焦冤沈钳荆疾仗煮侈泅俩袒烦仑陇援评伺淤翔刷锹鹊挂更罐雪还跑允馈死幢韶巫姿宇股寐墅主拘泄纪采压死洁焦故寝诬足息茵纬柞孟砂膨野爬湾显谜扯吹拷舵钮媒却瘦诲疡卑湍侨慈痈侄偶斡茶瞥减赏恐爬施娄疼砒苗甥风鸦Matlab编程实现FFT变换及频谱分析的程序代码（喜欢进行电磁兼容仿真分析的朋友可以借用）内容1用Matlab产生正弦波,矩形波,以及白噪声信号，并显示各自时域波形图2进行FFT变换，显示各自频谱图，其中采样率，频率、数据长度自选3做出上