gain.m

在matlab中实现数字信号处理的一些源程序

f=[0.000000 0.001672 0.003344 0.005017 0.006689 0.008361 0.010033 0.011706 0.013378 0.015050 0.016722 0.018395 0.020067 0.021739 0.023411 0.025084 0.026756 0.028428 0.030100 0.031773 0.033445 0.035117 0.036789 0.038462 0.040134 0.041806 0.043478 0.045151 0.046823 0.048495 0.050167 0.051839 0.053512 0.055184 0.056856 0.058528 0.060201 0.061873 0.063545 0.065217 0.066890 0.068562 0.070234 0.071906 0.073579 0.075251 0.076923 0.078595 0.080268 0.081940 0.083612 0.085284 0.086957 0.088629 0.090301 0.091973 0.093645 0.095318 0.096990 0.098662 0.100334 0.102007 0.103679 0.105351 0.107023 0.108696 0.110368 0.112040 0.113712 0.115385 0.117057 0.118729 0.120401 0.122074 0.123746 0.125418 0.127090 0.128763 0.130435 0.132107 0.133779 0.135452 0.137124 0.138796 0.140468 0.142140 0.143813 0.145485 0.147157 0.148829 0.150502 0.152174 0.153846 0.155518 0.157191 0.158863 0.160535 0.162207 0.163880 0.165552 0.167224 0.168896 0.170569 0.172241 0.173913 0.175585 0.177258 0.178930 0.180602 0.182274 0.183946 0.185619 0.187291 0.188963 0.190635 0.192308 0.193980 0.195652 0.197324 0.198997 0.200669 0.202341 0.204013 0.205686 0.207358 0.209030 0.210702 0.212375 0.214047 0.215719 0.217391 0.219064 0.220736 0.222408 0.224080 0.225753 0.227425 0.229097 0.230769 0.232441 0.234114 0.235786 0.237458 0.239130 0.240803 0.242475 0.244147 0.245819 0.247492 0.249164 0.250836 0.252508 0.254181 0.255853 0.257525 0.259197 0.260870 0.262542 0.264214 0.265886 0.267559 0.269231 0.270903 0.272575 0.274247 0.275920 0.277592 0.279264 0.280936 0.282609 0.284281 0.285953 0.287625 0.289298 0.290970 0.292642 0.294314 0.295987 0.297659 0.299331 0.301003 0.302676 0.304348 0.306020 0.307692 0.309365 0.311037 0.312709 0.314381 0.316054 0.317726 0.319398 0.321070 0.322742 0.324415 0.326087 0.327759 0.329431 0.331104 0.332776 0.334448 0.336120 0.337793 0.339465 0.341137 0.342809 0.344482 0.346154 0.347826 0.349498 0.351171 0.352843 0.354515 0.356187 0.357860 0.359532 0.361204 0.362876 0.364548 0.366221 0.367893 0.369565 0.371237 0.372910 0.374582 0.376254 0.377926 0.379599 0.381271 0.382943 0.384615 0.386288 0.387960 0.389632 0.391304 0.392977 0.394649 0.396321 0.397993 0.399666 0.401338 0.403010 0.404682 0.406355 0.408027 0.409699 0.411371 0.413043 0.414716 0.416388 0.418060 0.419732 0.421405 0.423077 0.424749 0.426421 0.428094 0.429766 0.431438 0.433110 0.434783 0.436455 0.438127 0.439799 0.441472 0.443144 0.444816 0.446488 0.448161 0.449833 0.451505 0.453177 0.454849 0.456522 0.458194 0.459866 0.461538 0.463211 0.464883 0.466555 0.468227 0.469900 0.471572 0.473244 0.474916 0.476589 0.478261 0.479933 0.481605 0.483278 0.484950 0.486622 0.488294 0.489967 0.491639 0.493311 0.494983 0.496656 0.498328 0.500000 ];
x=[0.052632 0.052634 0.052643 0.052658 0.052678 0.052704 0.052736 0.052774 0.052818 0.052867 0.052923 0.052984 0.053052 0.053125 0.053205 0.053290 0.053382 0.053480 0.053584 0.053695 0.053812 0.053936 0.054066 0.054202 0.054345 0.054496 0.054653 0.054816 0.054987 0.055166 0.055351 0.055544 0.055744 0.055952 0.056167 0.056391 0.056623 0.056862 0.057111 0.057367 0.057633 0.057907 0.058190 0.058482 0.058784 0.059096 0.059417 0.059749 0.060091 0.060444 0.060807 0.061182 0.061568 0.061966 0.062375 0.062798 0.063233 0.063681 0.064143 0.064618 0.065108 0.065612 0.066132 0.066668 0.067219 0.067787 0.068373 0.068976 0.069598 0.070239 0.070899 0.071580 0.072283 0.073007 0.073754 0.074524 0.075320 0.076140 0.076988 0.077863 0.078767 0.079700 0.080665 0.081663 0.082695 0.083762 0.084867 0.086011 0.087195 0.088423 0.089695 0.091014 0.092382 0.093803 0.095279 0.096812 0.098406 0.100064 0.101790 0.103588 0.105462 0.107416 0.109456 0.111587 0.113815 0.116145 0.118586 0.121144 0.123827 0.126646 0.129609 0.132727 0.136013 0.139479 0.143141 0.147014 0.151117 0.155470 0.160096 0.165019 0.170269 0.175878 0.181883 0.188325 0.195253 0.202721 0.210792 0.219540 0.229050 0.239420 0.250769 0.263235 0.276980 0.292202 0.309135 0.328063 0.349328 0.373347 0.400626 0.431777 0.467530 0.508736 0.556336 0.611243 0.674064 0.744489 0.820160 0.895023 0.958084 0.995069 0.995069 0.958084 0.895023 0.820160 0.744489 0.674064 0.611243 0.556336 0.508736 0.467530 0.431777 0.400626 0.373347 0.349328 0.328063 0.309135 0.292202 0.276980 0.263235 0.250769 0.239420 0.229050 0.219540 0.210792 0.202721 0.195253 0.188325 0.181883 0.175878 0.170269 0.165019 0.160096 0.155470 0.151117 0.147014 0.143141 0.139479 0.136013 0.132727 0.129609 0.126646 0.123827 0.121144 0.118586 0.116145 0.113815 0.111587 0.109456 0.107416 0.105462 0.103588 0.101790 0.100064 0.098406 0.096812 0.095279 0.093803 0.092382 0.091014 0.089695 0.088423 0.087195 0.086011 0.084867 0.083762 0.082695 0.081663 0.080665 0.079700 0.078767 0.077863 0.076988 0.076140 0.075320 0.074524 0.073754 0.073007 0.072283 0.071580 0.070899 0.070239 0.069598 0.068976 0.068373 0.067787 0.067219 0.066668 0.066132 0.065612 0.065108 0.064618 0.064143 0.063681 0.063233 0.062798 0.062375 0.061966 0.061568 0.061182 0.060807 0.060444 0.060091 0.059749 0.059417 0.059096 0.058784 0.058482 0.058190 0.057907 0.057633 0.057367 0.057111 0.056862 0.056623 0.056391 0.056167 0.055952 0.055744 0.055544 0.055351 0.055166 0.054987 0.054816 0.054653 0.054496 0.054345 0.054202 0.054066 0.053936 0.053812 0.053695 0.053584 0.053480 0.053382 0.053290 0.053205 0.053125 0.053052 0.052984 0.052923 0.052867 0.052818 0.052774 0.052736 0.052704 0.052678 0.052658 0.052643 0.052634 0.052632 ];
y=[3.141593 3.141040 3.140486 3.139933 3.139379 3.138825 3.138270 3.137715 3.137158 3.136601 3.136042 3.135482 3.134921 3.134359 3.133794 3.133228 3.132661 3.132091 3.131519 3.130944 3.130367 3.129788 3.129206 3.128621 3.128033 3.127442 3.126847 3.126249 3.125647 3.125042 3.124432 3.123818 3.123200 3.122578 3.121950 3.121318 3.120681 3.120038 3.119390 3.118736 3.118075 3.117409 3.116737 3.116057 3.115371 3.114677 3.113976 3.113268 3.112551 3.111826 3.111092 3.110350 3.109598 3.108836 3.108064 3.107282 3.106490 3.105686 3.104870 3.104043 3.103203 3.102350 3.101484 3.100604 3.099709 3.098800 3.097875 3.096933 3.095975 3.094999 3.094005 3.092992 3.091959 3.090905 3.089830 3.088733 3.087613 3.086468 3.085298 3.084101 3.082877 3.081623 3.080340 3.079025 3.077677 3.076294 3.074875 3.073418 3.071920 3.070381 3.068798 3.067168 3.065490 3.063760 3.061976 3.060134 3.058233 3.056267 3.054233 3.052129 3.049948 3.047686 3.045339 3.042901 3.040365 3.037727 3.034977 3.032109 3.029115 3.025984 3.022707 3.019273 3.015669 3.011882 3.007896 3.003695 2.999260 2.994570 2.989601 2.984327 2.978718 2.972740 2.966354 2.959516 2.952175 2.944273 2.935743 2.926504 2.916466 2.905518 2.893532 2.880354 2.865798 2.849641 2.831607 2.811360 2.788480 2.762442 2.732578 2.698039 2.657729 2.610230 2.553699 2.485762 2.403418 2.303069 2.180871 2.033802 1.861771 1.670284 1.471308 1.279822 1.107791 0.960721 0.838524 0.738175 0.655831 0.587894 0.531363 0.483863 0.443553 0.409015 0.379151 0.353112 0.330232 0.309985 0.291952 0.275794 0.261239 0.248061 0.236075 0.225127 0.215088 0.205850 0.197319 0.189417 0.182077 0.175239 0.168853 0.162874 0.157265 0.151991 0.147023 0.142332 0.137897 0.133696 0.129711 0.125924 0.122320 0.118886 0.115609 0.112478 0.109483 0.106616 0.103866 0.101227 0.098692 0.096254 0.093906 0.091645 0.089464 0.087359 0.085326 0.083360 0.081458 0.079617 0.077833 0.076103 0.074424 0.072795 0.071211 0.069672 0.068175 0.066718 0.065299 0.063916 0.062568 0.061253 0.059969 0.058716 0.057492 0.056295 0.055125 0.053980 0.052860 0.051762 0.050687 0.049634 0.048601 0.047588 0.046594 0.045618 0.044660 0.043718 0.042793 0.041883 0.040989 0.040108 0.039242 0.038389 0.037550 0.036722 0.035907 0.035103 0.034310 0.033528 0.032757 0.031995 0.031243 0.030500 0.029767 0.029042 0.028325 0.027616 0.026915 0.026222 0.025535 0.024856 0.024183 0.023517 0.022857 0.022203 0.021555 0.020912 0.020275 0.019642 0.019015 0.018392 0.017774 0.017161 0.016551 0.015945 0.015344 0.014746 0.014151 0.013560 0.012972 0.012387 0.011805 0.011225 0.010648 0.010074 0.009502 0.008932 0.008364 0.007798 0.007234 0.006671 0.006110 0.005550 0.004992 0.004434 0.003878 0.003322 0.002768 0.002213 0.001660 0.001106 0.000553 0.000000 ];
subplot 211; 
plot(f,x);
xlabel('归一化频率');
ylabel('幅度增益');
title('数字系统的幅频响应');
subplot 212; 
plot(f,y);
xlabel('归一化频率');
ylabel('相移');
title('数字系统的相频响应');