毕业回老家-CSDN博客

原创 from Crypto.Cipher import AES 报错

报错原因：默认下载的Crypto是最新版本1.4.1，可以打开其内部关于aes模块的文件，其实什么代码都没有，所以我们安装旧版本就可以了。解决方案：安装低版本的Crypto，比如1.3.3 ，如下图；

2023-05-01 19:32:40 814 1

原创基于模糊C均值聚类（FCM）的图像分割原理+python代码详解

“模糊”：一个元素可以不同程度的属于某几个子集，也就是说元素对于集合的隶属度可以在[0,1]上取连续值。翻译一下： S1：初始化参数：加权指数m，簇心数目C，以及迭代停止阈值ε。S2：随机初始化隶属度矩阵U，注意满足式（2-2）。S3：式（2-3）更新簇心c。S4：式（2-4）更新隶属度矩阵US5：如果隶属度矩阵U满足式（2-5）则返回U并结束算法，否则转到S2各种式如下：流程图比步骤多的一项是计算目标函数J，这一步可有可无，为了观察目标函数的变化趋势，我们在流程图和代码里都加上了计算目标

2022-07-11 18:57:54 4853 5

原创基于K-means的图像分割（python代码详解）

基于K-means的图像分割+代码详解

2022-07-11 15:35:43 13552 3

原创基于Qt平台利用UDP协议传输大文件（百度网盘免费代码）

百度网盘分享代码链接：https://pan.baidu.com/s/1cwcAH4ugobYqQ0c4UYdCiw提取码：1234

2022-06-04 10:16:10 4015 4

原创基于Qt平台利用TCP协议传输文件（百度网盘分享项目）

与UDP协议单次传输数据大小要小于64KB不同，TCP协议利编程使用的是TCP SOCKET套接字。项目文件用百度网盘分享：链接：https://pan.baidu.com/s/1kkDPIsK4jBxM1_m9iGs05Q 提取码：12341. 客户端UI 有几点注意事项： 1.客户端如果要传输文件到服务器端，需要指明服务器端的IP和Port。 ......

2022-06-04 10:02:10 2259

原创基于Qt平台利用UDP协议传输文件

毕业设计用到了Qt，需要用UDP协议组播传输文件。现在介绍一下遇到的一些坑。1、传输文件一般应该用TCP协议，但是我毕设的需求是要组播分发视频，所以必须用UDP协议传输文件。2、利用udpsocket创建套接字对象，利用成员函数writeDatagram能传输的包<64kb。不过TCP协议的套接字对象tcpsocket的成员函数write可以没有大小的限制。3、等最近实现UDP传大文件，会把这些TCP/UDP传输文件的项目放到网盘，大家可以看看，都是比较简单的东西。https://

2022-04-11 17:28:59 4369

原创 matlab 幅频特性（dB）y =20*log10(y)

单位转换为分贝(dB)1、20lg(|H|)，底数为10.2、matlab中的对数函数：log10()：即lglog()：以自然指数e为底log2()：以2为底3、matlab画图y = fft(x)y = abs(y)y = log10(y)

2022-01-11 15:41:43 9080

原创 matlab filter freqz

1.滤波函数filter利用滤波器H(Z)=b/a，对X滤波ba为向量滤波器传输函数H（Z）的分子分母：b,a利用b，a表征H即滤波器2.freqzH=b/a，由b和a提取H和数字角频率w2.1[h,w] =freqz(b,a,n)n默认512，w=[0,pi]2.2[h,w] = freqz(b,a,'whole')w=[-pi,pi];传输函数是对称函数，幅频响应|H|关于w=pi偶对称，相频特性关于w=pi...

2022-01-11 15:35:49 1885

原创数字信号频域分析：DFT与FFT

一、前言连续信号的频域分析：将信号分为周期与非周期利用FS和FT分析，连续信号的频谱是关于纵轴对称的，且是非周期的。但是计算机不能处理连续信号，只可处理数字信号。所以数字信号频谱特性更重要：二、DFT(离散傅里叶变换) 数字信号的频谱是周期信号，周期是 2π rad/s（数字角频率），N点DFT是对数字信号的DTFT（傅里叶变换）在[0，2π]上的N点等间隔采样。FFT是在DFT的基础上简化了计算，提高了速度。...

2022-01-10 15:32:50 5922

原创 matlab：卷积function

function [n,y] = convm(x,nx,h,nh)y = conv(x,h);n = [nx(1)+nh(1):nx(1)+nh(1)+length(y)-1];end卷积后序列的起始序号 = 两个卷积序列起始序号和

2022-01-10 14:51:39 456

原创 matlab：常见周期信号

1.锯齿波：x =sawtooth(2*pi*f*t,xmax)产生频率为f，最大值在xmax*T，幅度在[0，1]的锯齿波。2.方波信号x =square(t)x =square(t,duty)T=2*pi，duty为占空比[0,100]单位是%，即一个周期内方波存在的时间与一个周期的比值，幅度[-1,1]占空比默认是50%...

2022-01-10 14:17:26 2648

原创连续信号频域分析：功率信号的频谱与能量信号频谱密度

一、周期信号的频谱1.1周期信号的虚指数/三角傅里叶级数：其中，1.2 意义：从频域上看，周期信号可以由其所有频率分量叠加而成，并且这些频域分量都是周期信号频率的整数倍，傅里叶级数的作用就是算出这些不同的频率分量的幅度和相位。1.3所以，周期信号的频域分析就是分析频谱：幅度谱和相位谱。两种完备集对应的频谱关系如下： (1)三角函数：频谱都是单边谱(2)虚指数：频谱都是双边.........

2022-01-09 21:01:40 8253

原创基于matlab的多路BPSK调制—为什么采样率=码元速率x一个码元周期内对码元的采样点数

一、为什么“fs = Rb x N” 其中， Rb（码元速率）、fs（采样率）、N（一个码元周期内对码元的采样点数）先说结论，因为调制时，信源序列与调制信号（正弦信号）的点数必须一样，尤其是多路信源序列一起调制时，因为最终的已调信号储存在一个矩阵的不同行里，而矩阵不同行的列数必须一样，也就是说不同码元速率的码元系列，最终的采样点数都是一样的。二、对码元速率分别为1Msymbol/s、2Msymbol/s 4Msymbol/s 和5Msymbol/s的信源经过汉明编码后......

2021-12-21 20:01:39 5972 1

原创基于matlab的多路信道编码——汉明码

我们分两个方面分析，一是信道编码原理，二是matlab仿真。一、信源产生了什么？信源的作用：将信息转换成电信号，数学上信源的输出是一个随时间变换的随机函数。也就是说信源序列是一组随机序列数，可以用特定函数产生信源序列：1、randi([a,b],m,n)：产生一个m x n的均匀分布随机矩阵，元素的值属于[a,b]且取整数。2、randn(m,n) > 0.5 函数产生m x n的正态分布的随机矩阵，元素值在[0,1......

2021-12-21 19:30:11 4751 1

原创 Python2.Python与matlab自定义函数(封装)

函数能提高应用的模块性，和代码的重复利用率。你已经知道Python提供了许多内建函数，比如print()。但你也可以自己创建函数，这被叫做用户自定义函数。1、Python：参考Python3 函数 | 菜鸟教程1.1函数定义：def 函数名（x1,x2...）：(1).函数第一行一般写注释(2).return来返回输出。如果没有输出不用写return。如果没有输入：def 函数名（），即空括号egdef max(a,b): if a>=b:

2021-12-07 19:08:28 593

原创 Python笔记1.一篇文章初步入门

现在开始学Python吧！就拿runoob.com入门吧Python3 简介 | 菜鸟教程从我认为最重要的几点开始学：1、输入输出2、数据类型3、程序的基本结构：循环and分支

2021-12-07 13:02:06 1233

原创信噪比(dB)换算公式

SNR(dB) = 10*lg(SNR);1、举例：0(dB) = 1 =10^0;10(dB) = 10 =10^1;20(dB) = 100 =10^2;30(dB) = 1000 =10^3;即10*n (dB) = 10^n;2、转换成dB的意义压缩坐标轴。3、图像SNR(dB)与SNR对比可以看，SNR范围从0-1000，SNR(dB)的最大值只要到404、SNR（dB）作为自变量作图，就不能再用plot了semilogx(X,Y) 在 ..

2021-12-05 18:15:55 70264

原创基于matlab的信道编码实验：多维连续信源与信道的信息量、多维连续信道的信道容量

当年matlab的一个很简单的课程实验，现在拿出来分享一下：一、特殊N维连续信源的熵1、均匀分布的连续信源：1.1一维：1.2N维若N维随机变量X=(X1,X2....XN)中各分量彼此统计独立，且分别在[a1,b1]、[a2,b2].....[aN,bN]的区域内均匀分布，则信源熵为：因此只要知道各随机变量的取值区间，分别求出各个变量在各自区间内均匀分布时的相对熵，再将这些相对熵进行累加求和，就能计算出多维连续信源的相对熵。实现代码；%%特殊N维连续信...

2021-12-05 11:27:28 3066 1

原创 matlab8.input输入函数fprintf、disp输出函数

一、inputx = input(prompt)prompt是文本即字符串，用来提示输入。eg.prompt = 'What is the original value? ';x = input(prompt)y = x*10二、disp输出1、disp(X)显示x的值，而不打印变量名称。即不出现‘x’/2、ans= disp(x)显示变量的另一种方法是键入它的名称，这种方法会在值前面显示一个前导，比如'ans ='三、fprintffprintf(.

2021-12-05 10:41:23 2962

原创 matlab通信系统建模3.信道噪声

1、agwny = awgn(x,SNR)给信号x中加入高斯白噪声，信噪比为SNR(dB)clc;clear;t = 0:0.01:10;x = sawtooth(t);%锯齿波y = awgn(x,10);%加入高斯白噪声plot(t,x,t,y);legend('原信号','叠加高斯白噪声信号');2、噪声发生器2.1均匀分布随机噪声产生器noise = a + (b - a)*rand(m,n)产生在（a,b）服从均匀分布的随机噪声2.2高斯随机.

2021-12-04 14:07:53 2103

原创 matlab通信系统建模2.信源编解码

信源编码的三种常见方式：算术编码、哈夫曼编码、LZ编码1、算术二进制编解码算数编码实际的编译码过程比较复杂，但在性能上具有许多优点，特别是所需要的参数很少，不像哈夫曼编码那样需要一个很大的码表。从理论上说，只要已知信源符号集及其符号概率，算数编码的平均码长可以接近符号熵。(1)编码：arithencocode = arithenco(seq,counts)code：算术编码序列seq：信源产生的待编码序列counts：信源产生序列的不同符合出现次数clc;clear;

2021-12-04 11:40:18 4917

原创 matlab通信系统建模1.信源函数

matlab对通信系统仿真很方便，大多数操作都有对应的函数命令，比较难的是知识点的掌握。下图是数字通信系统模型：下面介绍信源函数的产生。信源的作用：将信息转换成电信号，数学上信源的输出是一个随时间变换的随机函数。提到随机函数，就可以使用rand系列的函数1、randirandi([a,b],m,n)：产生一个m x n的均匀分布随机矩阵，元素的值属于[a,b]且取整数。2、randnrandn(m,n)> 0.5 函数产生m x n的正态分布的随机矩阵，..

2021-12-04 10:23:43 2295

原创 matlab 5 .封装：幅频特性

%x:待处理信号%fs:采样率%N:采样点数function plot_fft(x,fs,N)y = fft(x);y = abs(y);f = (0:N-1)*fs/N;figure,plot(f(1:N/2),y(1:N/2)*2/N),grid on,xlabel('f/Hz'),ylabel('|fft(x)|'),title('信号幅频特性'),grid on;end%% 测试代码clc;clear;close all;N=2048;%采样点数fs=2048;%采.

2021-12-02 17:32:16 2102

原创 MATLAB 6.封装相频特性

%y = fft(x);%N:采样点数%fs:采样率function plot_fft_angle(y,fs,N)f = (0:N-1)*fs/N;ph = 2*angle(y(1:N/2));plot(f(1:N/2),ph(1:N/2));end

2021-12-02 17:21:45 421

原创采样率Fs与采样点数N对数字信号的时域波形和频谱的影响

在数字信号处理中，我们所用的数字信号大多是由模拟信号采样得到的，我们在利用matlab plot函数将这些离散点连成曲线时，会发现采样率与采样点数会对波形有很大的影响，我们下面来分析一下具体影响，并分析一下对于不同信号，Fs与N取多少合适：一、结论：Fs,N对数字信号时域波形和频谱的影响时域频域 Fs 分辨率分辨率 N 采样时长分辨率、幅度二、分析2.1数字频率与模拟频率的关系重点在数字频率的物理意义：(1).数字角频率与模拟

2021-12-01 22:43:16 19263

原创 matlab 5. 信号频谱分析封装

%x:待处理信号%fs:采样率%N:采样点数function plot_fft(x,fs,N)y = fft(x);y = abs(y);f = (0:N-1)*fs/N;figure,plot(f(1:N/2),y(1:N/2)*2/N),grid on,xlabel('f/Hz'),ylabel('|fft(x)|'),title('信号幅频特性');end测试代码：clc;clear;close all;N = 256;%采样点数fs = 128;%采样频率%% 时域

2021-12-01 20:53:24 2428

原创 matlab 4.BEP(带阻)封装巴特沃斯带阻滤波器

%x：待处理信号%fs：采样频率%wp;通带边界频率%ws:阻带边界频率%rp:Rp通带最大衰减(dB)%rs:Rs阻带最小衰减(dB)function y=BPF(x,fs,wp,ws,rp,rs)[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率[b,a] = butter(n,w,'stop');%确定H(S)的分子分母系数向量[H,W] = freqz(b,a);%H即频率响应函数，W即数字角频率figure,plo

2021-11-30 21:00:27 437

原创 matlab 3.BPF封装巴特沃斯带通滤波器

%x：待处理信号%fs：采样频率%wp;通带边界频率%ws:阻带边界频率%rp:Rp通带最大衰减(dB)%rs:Rs阻带最小衰减(dB)function y=BPF(x,fs,wp,ws,rp,rs)[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率[b,a] = butter(n,w);%确定H(S)的分子分母系数向量[H,W] = freqz(b,a);%H即频率响应函数，W即数字角频率figure,plot(W*fs/

2021-11-30 20:58:49 930

原创 matlab 2.封装HPF 巴特沃斯高通滤波器

%x：待处理信号%fs：采样频率%wp;通带边界频率%ws:阻带边界频率%rp:Rp通带最大衰减(dB)%rs:Rs阻带最小衰减(dB)function y=HPF(x,fs,wp,ws,rp,rs)[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率[b,a] = butter(n,w,'high');%确定H(S)的分子分母系数向量[H,W] = freqz(b,a);%H即频率响应函数，W即数字角频率figure,plo

2021-11-30 20:24:59 712

原创 matlab 1.封装LPF 巴特沃斯低通滤波器

%x：待处理信号%fs：采样频率%wp;通带边界频率%ws:阻带边界频率%rp:Rp通带最大衰减(dB)%rs:Rs阻带最小衰减(dB)function y=LPF_function(x,fs,wp,ws,rp,rs)[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率[b,a] = butter(n,w);%确定H(S)的分子分母系数向量[H,W] = freqz(b,a);%H即频率响应函数，W即数字角频率figure,p

2021-11-30 19:19:50 404

原创基于MATLAB的LPF低通滤波器设计(巴特沃斯)

clc;clear;close all;N=500;%采样点数fs=100;%采样频率%% 时域波形t = (0:N-1)/fs;%以fs采N个点x = sin(2*pi*10*t) + sin(2*pi*45*t);figure(1),subplot(211),plot(t,x),xlabel('t'),ylabel('x'),title('时域波形'),axis([0 1 -2 2]);%时域波形图%% 频谱n = length(x);f = (-n/2:n/2-1)*fs/n;%

2021-11-30 17:03:44 4825 1

marujie123的博客