
MATLAB
Хайде
愿你对未来满怀期待,也能活在当下,把握现在~
展开
专栏收录文章
- 默认排序
- 最新发布
- 最早发布
- 最多阅读
- 最少阅读
-
正弦信号的频谱分析
代码】正弦信号的频谱分析。原创 2022-08-12 18:36:16 · 5911 阅读 · 0 评论 -
基2fft实现、二次最佳平方逼近多项式、牛顿迭代法
根据信号与系统的相关理论[4],根据本学期的数字信号处理作业实验得出的相关结论,对正弦信号进行采样,采样频率在满足采样定理的前提下,还需要满足信号频率不能是信号频率的2倍,这样会导致采样点信息一致,不能恢复出原始的信息。正弦信号2倍采样结果如图2.3所示,当 的时候,会在信号的每一个周期采两个点,采到的所有点值一致,故不能得到正确的频谱。当 固定, 减小信号频谱特性会变好。由FFT基本原理可以知道,当N为2的幂次方时,FFT运算包含 级运算,并将每一级编号 ,每一级包含 组,每组包含 个鲽形单元,...原创 2022-08-12 18:25:25 · 1911 阅读 · 1 评论 -
正弦信号滤波、交通流信号滤波
信号滤波,周期计算原创 2022-08-12 18:16:13 · 831 阅读 · 0 评论 -
K均值据类C语言和matlab实现
C语言或matlab实现K均值聚类原创 2022-08-12 18:03:49 · 461 阅读 · 0 评论 -
以zigzag方式对矩阵数据遍历
遍历顺序:代码: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原创 2022-01-05 13:13:30 · 376 阅读 · 0 评论 -
将N个学生分为G组
代码:% %%考试分组%%考试分组%随机性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);%最后一组的人原创 2021-12-02 15:57:22 · 121 阅读 · 0 评论 -
求解二次最佳平方逼近多项式
求解二次最佳平方逼近多项式%以后没事不要编程解决数学练习册上边的题目了%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原创 2021-11-27 11:38:28 · 261 阅读 · 0 评论 -
高斯顺序消元解线性方程组
%高斯消元法解线性方程组%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原创 2021-11-25 22:29:50 · 68 阅读 · 0 评论 -
LU分解求解线性方程组
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原创 2021-11-25 11:23:48 · 80 阅读 · 0 评论 -
编程实现fft
结果代码: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原创 2021-11-24 22:25:11 · 126 阅读 · 0 评论 -
余弦信号的频谱分析
本博客讨论了以下问题:采样点数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原创 2021-11-18 19:51:57 · 334 阅读 · 0 评论 -
复化辛普森公式计算定积分
求积分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(原创 2021-11-10 17:02:22 · 348 阅读 · 0 评论 -
使用龙贝格方法求定积分
使用龙贝格方法求定积分本代码是模拟对函数在区间[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+原创 2021-11-09 23:11:26 · 151 阅读 · 0 评论 -
使用霍夫变换来检测直线
代码: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原创 2021-11-04 21:11:30 · 80 阅读 · 0 评论 -
用均值滤波对含噪信号进行处理
利用matlab对含噪正弦波信号进行滤波处理这是本人的一个小作业,程序中可能会有各种各样的小错误。首先生成正弦信号,分别设置不同的信噪比,设计滑动平均滤波器,对含噪信号进行滤波处理。对固定的信噪比,分析讨论滑动平均滤波器阶数对滤波性能的影响,计算降噪后信号与真实信号的均方根误差,绘制均方根误差随滤波器阶数变化的关系曲线,并绘制出最佳滤波器阶数下的滤波信号与原始信号以及含噪信号的对比曲线。原理:首先生成一个正弦波信号,设置相应的参数,将每一个参数设置成变量,同时给信号叠加噪声,使用awgn函数对信号进原创 2021-10-24 20:16:37 · 237 阅读 · 0 评论