weixin_40933496的博客

代码】正弦信号的频谱分析。

根据信号与系统的相关理论[4]，根据本学期的数字信号处理作业实验得出的相关结论，对正弦信号进行采样，采样频率在满足采样定理的前提下，还需要满足信号频率不能是信号频率的2倍，这样会导致采样点信息一致，不能恢复出原始的信息。正弦信号2倍采样结果如图2.3所示，当 的时候，会在信号的每一个周期采两个点，采到的所有点值一致，故不能得到正确的频谱。当 固定， 减小信号频谱特性会变好。由FFT基本原理可以知道，当N为2的幂次方时，FFT运算包含 级运算，并将每一级编号 ，每一级包含 组，每组包含 个鲽形单元，...

信号滤波，周期计算

C语言或matlab实现K均值聚类

遍历顺序：代码：clear,clc;%实现zigzag方式提取矩阵的的元素len1=4;len2=4;m1=cell(len1,len2);m2=cell(1,len1*len2);m3=ones(len1,len2);for m4=1:len1 for lp2=1:len2 m1{m4,lp2}=[char('A'+(m4-1)*len2+lp2-1)]; endendm5=1;m6=1;m3(1,1)=0;m2(1,1)=m1(1,1);c

代码：% %%考试分组%%考试分组%随机性N=69;%69个人G=4;%分组数目data=1:N;%每个人的初始索引%%给每个人随机分配索引data_index=zeros(1,N);for loop=1:N index=randi(N-loop+1); data_index(loop)=data(index); data=data(find(data~=data(index)));enddata_indexlast_num=fix(N/G);%最后一组的人

求解二次最佳平方逼近多项式%以后没事不要编程解决数学练习册上边的题目了%49-7 38-5都可以用clear,clc%求解二次逼近问题N=2;low=0;high=1;%积分的区间syms x;syms y;y=1/(1*x^2+1);p_x=x;%积分的权值a=zeros(1,N);b=zeros(1,N-1);g0=1;%要用元胞数组来定义吗？% g=1;% g=[cell(g0),cell(g0),cell(g0)];%暂时只能先用这个办法了for loop=1:N+1

%高斯消元法解线性方程组%A=[1 2 3;4 5 6;7 8 9];A=[3 -1 4;-1 2 -2;2 -3 -2];b=[7 -1 0]';A(:,end+1)=b; B = gauss_matrix( A ) %求解线性方程组 b=B(:,end); A=B(:,1:end-1); len1=length(A); x=[zeros(1,len1-1),b(end)/A(end)]; for loop=len1-1:-1:1 x(loop)= ( b(loop)-A(l

LU分解求解线性方程组：clear,clc,close all;% A=[1,0,2,0;0,1,0,1;1,2,4,3;0,1,0,3];% b=[5 3 17 7]';A=[-1 0 -1 3;1 1 4 -1;2 4 1 0;6 2 1 -1];b=[-5 5 -1 6]';[ L,U ] = LU( A );%然后通过LU分解求矩阵方程%首先求解Ly=blen2=length(L);y=[b(1,1)/L(1,1),zeros(1,len2-1)]';for loop3=2:l

结果代码：clear,clc,close all;%自己编程实现时间抽取基2fft算法f=50;%频率fs=10000;%采样频率Ts=1/fs;%采样间隔N=128;n=1:N;data1=sin(2*pi*f*n*Ts);% figure;% plot(n*Ts,data1);title('原始信号');data1;%原始信号%%将数据长度变为2的幂次方if(mod(log2(N),1)~=0) temp=ceil(log2(N)); data1=[data

本博客讨论了以下问题：采样点数N以及采样频率Fs对信号频谱的影响？最大频谱所对应频率？产生频率为20Hz和20.02的两个余弦信号的和信号，并分析信号加窗对频率分辨率的影响？结果1：结果2：结果3：结果4：结果5：代码：clear,clc,close all;t_start=tic;%数字信号处理第二次作业tl=0.5;f1=20;%频率f2=20.02;fs=60;N=60;O0=0;cnt1=1;fs=50;for N=10:30:1000 d

求积分clear,clc,close all;%f_function输入x的取值（支持输入向量）输出函数的值%T_function输入[a,b,n]分别为积分的起始点、结束点、份数E4.1a=0b=1;n=8;h=(b-a)/n;data=a:h:b;%使用复化梯形公式来计算data_T=[data;f_function(data)];T_sum=h/2*(data_T(2,1)+sum(2*data_T(2,2:end-1))+data_T(2,end));disp(strcat(

使用龙贝格方法求定积分本代码是模拟对函数在区间[0,1]上求定积分clear,clc,close all;%用龙贝格算法计算定积分%f=4/(1+x.*x);error_F=1/2*10^-5;%精确度a=0;b=1;%改变函数需要更改后面的T_function里边的函数flag=1;loop=0;while(flag||abs(data(loop,5)-data(loop-1,5))>error_F) n=2^loop;%n=2^loop; data(loop+

代码：clear,clc,close all;%霍夫变化的仿真 %data=[1 2 3;-1 -2 -3]; data=[-1 0 1 2 0 2 3 8 9 6 -1 -3 -5 -7 6 9 8 7 1.1 2.2 3.1;3 2 1 0 0 2 3 8 9 6 -1 -3 -5 -7 0 2 5 9 1 2 3]; %data=[0:10;1:-1:-9];% data=xlsread('data.xls');% data=data';%data=[1 2 3 5 6;2 4

利用matlab对含噪正弦波信号进行滤波处理这是本人的一个小作业，程序中可能会有各种各样的小错误。首先生成正弦信号，分别设置不同的信噪比，设计滑动平均滤波器，对含噪信号进行滤波处理。对固定的信噪比，分析讨论滑动平均滤波器阶数对滤波性能的影响，计算降噪后信号与真实信号的均方根误差，绘制均方根误差随滤波器阶数变化的关系曲线，并绘制出最佳滤波器阶数下的滤波信号与原始信号以及含噪信号的对比曲线。原理：首先生成一个正弦波信号，设置相应的参数，将每一个参数设置成变量，同时给信号叠加噪声，使用awgn函数对信号进