weixin_42845306的博客

原创 LINUX系统SSH连接安装Matlab、添加环境变量、后台运行Matlab脚本

网上找一个Matlab安装包（iso镜像）上传进服务器端，或者直接使用wget指令在服务器端下载镜像，一定注意要下载Linux版本的而不是window版本的。双击打开，在最后面加上一句话，将matlab的路径添加到系统变量，改完保存，重新登录服务器。通过以上步骤，你可以在任意位置输入matlab回车，即可打开matlab运行界面。即可打开可视化安装界面，一步一步操作即可安装完成（我用的是学校邮箱激活）。若存在挂起的进程，需要手动找到并杀掉后exit。好，可以安心的关机去干别的事情了。

原创 ADS 使用调整器动态调参教程（Tuning）

ADS动态调参教程

原创 matlab 精确线搜索代码：进退法求初始区间、黄金分割法求最小值

精确线搜索、进退法、黄金分割法

原创 matlab 二阶导（海森矩阵）的数值计算（附代码和示例）

海森矩阵数值计算

原创 matlab 导数和梯度的数值计算（附代码和示例）

导数和梯度的数值计算

原创 MATLAB 绘制信号单双边频谱代码

【代码】MATLAB 绘制信号单双边频谱代码

原创 Latex 公式图表按章节编号、顺序编号

latex最前面加上三句代码：\documentclass[a4paper]{article}%核心就是下面三段代码\usepackage{amsmath}\numberwithin{equation}{section}%公式按章节编号\numberwithin{figure}{section}%图表按章节编号%在其他文章类型中，section可以换为chapter\begin{document} \section{principle} %这里只演示公式 \begin{equatio

原创 Matlab APP designer 里的grid on 和hold on

众所周知，MATLAB里显示图窗网格的命令是：grid on连续绘图的命令是hold on然而在MATLAB新推出的图形化设计平台 APP designer中，这两个命令并没有用。由于该平台较新，网上教程基本没找到，还好mathwork把手册更新的比较完善，可以自学。grid on方法：点击坐标轴，右下角对象属性里把网格勾上即可：编程控制的话，就是这两行命令：app.UIAxes.XGrid='on';app.UIAxes.YGrid='on';hold on方法：触类旁通，把图窗的属

原创 电脑连接示波器读取数据

示波器型号：北京普源的RIGOL MSO4034首先连接示波器：连接线与示波器相连的端口如上图USB DEVICE所示，是一个比较奇怪的方形接口。连接线的另一端是标准的USB typeA接口，可以直接插在电脑上，当示波器和电脑连接时，电脑会提示安装驱动。该型号示波器端口是业界通用的标准驱动VISA设计，可以使用SCPI命令控制。首先下载NI的Ultra Sigma软件，安装，启动。下载链接：https://pan.baidu.com/s/1K8u5ZoB2D31aix76tGqZcg提取码：ng

原创 K-L变换原理、图像压缩与人脸识别实战（附matlab代码）

目录摘要K-L变换的基本原理K-L变换在图像压缩中的应用研究背景及意义数字图像压缩与发展现状K-L变换与图像压缩K-L变换在人脸识别中的应用人脸识别K-L变换实现人脸识别总结参考文献matlab代码致谢摘要K-L变换的基本原理K-L变换在图像压缩中的应用研究背景及意义数字图像压缩与发展现状K-L变换与图像压缩K-L变换在人脸识别中的应用人脸识别K-L变换实现人脸识别总结参考文献matlab代码1. clear 2. %set the

原创 【近万字】分数傅里叶变换课程学习笔记

学习自“课堂在线”平台，北京理工大学陶然教授的课程视频，讲解的非常详细全面，数学公式推导都有，以下为学习笔记，仅记录要点部分。注：学习此课程，按重要程度排序，需要有信号与系统、数字信号处理、随机信号分析、矩阵论等课程的基础。目录绪论分数傅里叶变换定义分数傅里叶变换的性质分数卷积与滤波卷积定理卷积的快速计算方法分数卷积功率谱分数功率谱分数域采样与重建分数域谱分析分数域采样后重建频域带通采样定理：分数域带通采样定理离散FRFT快速算法chirp信号检测与参数估计参数估计公式量纲归一化后记绪论傅里叶变换处

原创 基于分数阶傅里叶变换的chirp信号检测与参数估计（原理附代码）

线性调频信号（chirp信号）顾名思义，该信号的频率随着时间线性变换，其复数表达形式如下：s(t)=e2jπ(f0t+0.5μt2)s(t)=e^{2j\pi(f_0 t+ 0.5\mu t^2)}s(t)=e2jπ(f0​t+0.5μt2)根据欧拉公式，其相位项为2π(f0t+0.5μt2))2\pi(f_0 t+ 0.5\mu t^2))2π(f0​t+0.5μt2))。信号的角频率是相位对时间的导数。对相位求导后有2π(f0+μt)2\pi(f_0+\mu t)2π(f0​+μt)，显然，f0f

原创 分数阶傅里叶变换（FrFT）详细原理与matlab代码实现

本文主要是基于Haldun M. Ozaktas, Orhan等人的论文Digital Computation of the Fractional Fourier Transform 翻译而成。如有错误之处还请指出。摘要本文给出了一种高效、精确计算分数傅里叶变换的算法。对于时间带宽积为N的信号，该算法计算复杂度为O（N logN）。我们还讨论了离散分数傅里叶变换的定义。引言略……预备知识分数傅里叶变换定义{Ff}(x)\{\mathcal{F}f\}(x){Ff}(x)表示对f(x)f(x)f

原创 压缩感知与稀疏恢复 matlab示例（附代码）

前言最近因为科研需要，又开始重新研究压缩感知（CS）与稀疏恢复（SR）理论。本人系初学，很多东西都没有学明白，姑且先摸着石头过河，仿照网上的例子，用matlab编程实现最基本的例子，写下了这篇笔记。原理关于压缩感知与稀疏恢复的原理就不再赘述，网上有很多博主写的很详细，这里推荐https://zhuanlan.zhihu.com/p/22445302这篇文章原理写的通俗易懂。本文在这篇文章的基础上结合https://blog.csdn.net/xiahouzuoxin/article/details/

原创 电磁仿真软件CST2020 下载链接及安装破·解教程笔记

众所周知，配环境是安装软件最头疼的事情。CST2020我调了一上午，网上现有的安装教程良莠不齐，很难有一篇文章能够全流程的讲述安装过程以及遇到的问题。正好今天我要给另一台电脑装软件，借此机会把安装过程记录下来，以供将来使用和他人借鉴。首先提供下载地址：百度网盘链接：https://pan.baidu.com/s/1cDJZelufj-nliUfekMNpEA提取码：8ihx附：根据二○○二年一月一日《计算机软件保护条例》规定：为了学习和研究软件内含的设计思想和原理，通过安装、显示、传输或者存储软

原创 MATLAB PDE工具箱电磁场仿真平行电容板及电偶极子详解

首先打开matlab，在命令行输入pdetool打开PDE工具箱工具箱的功能我就不赘述了，大概就是以有限元的数值计算方法，计算二阶偏微分方程，因此可以解决电磁场、热力学等问题。首先点击左上角按钮，画三个矩形画完一个后可以双击矩形编辑，参数如下：其中R1表示上极板，R2表示下极板，R3表示整个空间区域。在上方Set formula为R3-R2-R1，表示求解区域是整个空间区域减去两个极板。下拉框选择Electrostatics表示求解的是静电场问题。接下来点击上方菜单栏中的这个按钮，表

原创 MATLAB 四元数旋转函数

%四元数旋转法，输入theta是旋转角度（右手系），v是旋转轴矢量，u是原始点位置矢量，输出U是旋转后的点位置矢量function U=QuaternionRotate(theta,v,u) v=v./sqrt(sum(v.*v));%归一化 q=[cos(theta/2),v(1)*sin(theta/2),v(2)*sin(theta/2),v(3)*sin(theta/2)]; qt=quatconj(q);%取共轭 w=[0,u];%纯四元数 W=quatm

原创 MATLAB 绕轴旋转图像函数 rotate 详解

介绍以指定原点和方向旋转图像对象。语法rotate(h,direction,alpha)rotate 函数在三维空间中旋转图形对象。rotate(h,direction,alpha) 将图形对象 h 旋转 alpha 度。指定 h 为曲面、补片、线条、文本或图像对象。direction 是一个二元素或三元素向量，它与旋转轴原点共同确定旋转轴。旋转轴的默认原点是图框中心。该点未必是坐标区的原点。正 alpha 定义围绕方向向量从旋转原点伸展的右旋角度。如果 h 为句柄数组，则所有对象必须为相同坐

原创 图像处理课程设计大报告 MATLAB GUI APP实现直方图均衡化、几何变换和加噪滤波

单独写了一篇文章：https://blog.csdn.net/weixin_42845306/article/details/117950570

原创 直方图均衡化 原理、流程、公式推导及matlab实现

附上按照上述流程实现的MATLAB代码：clc;close all;clear all; H= imread('测试图\lena512.bmp'); H=H(:,:,1);%判断是否为三通道彩色图片 若是 则将其灰度化if length(size(H))>2 H=rgb2gray(H); end%获取图片的尺寸 便于计算总像素数 即m*n [m,n]=size(H); %生成一个一行256列的矩阵p=zeros(1,256); % 统计各灰度的像素个数...

原创 Matlab 数字滤波器设计大报告（数字信号处理课程设计）附代码

数字滤波器是由数字乘法器、加法器和延时单元组成的一种算法或装置，随着数字技术的发展而得到了广泛应用，其功能是对输入离散信号的数字序列进行运算处理，以达到改变信号频谱的目的。 本文构建了脑电信号检测的简易模型，设计了一种数字巴特沃斯低通滤波器，旨在滤除工频信号干扰，保留脑电信号。 本文首先建立了信号模型并对信号进行了分析。然后提出了滤波器设计指标，设计出模拟巴特沃斯滤波器。接着使用双线性变换，将模拟滤波器转换为数字滤波器，画出了滤波器信号流图。最后使用滤波器对信号进行了滤波，分析了滤波结果。

原创 巴特沃斯归一化数字低通滤波器参数、原理及使用方法

记录一下，免得使用的时候忘记。数字信号处理（第四版），西安电子科技大学出版社，高西全 丁玉美等编著，161页原理和使用方法见前两页

原创 DSP双线性变换详解

原创 tikz画GPS卫星电文通信结构图

TeX代码如下：

原创 MATLAB 二脉冲对消、重频参差滤波器设计仿真

雷达原理与系统作业，仅供参考，欢迎讨论，杜绝完全照搬。一、 问题提出二、问题解决1.二脉冲对消原理2.重频参差滤波器原理：3.仿真实验Matlab代码：c=3e8;%光速PRF=2000;%脉冲重复频率Tr=1/PRF;%脉冲重复间隔two_pulseH=@(w) 1-exp(-1j.*w.*Tr);%二脉冲对消系统函数Doppler_w=@(v,lambda) 2.*v./lambda;%多普勒频率lambda=c/(3e9);%S 波段雷达 2

原创 matlab imadjust将暗图像分别在RGB与HSV域增加亮度

对于RGB来说，增加亮度是直接对三个颜色维做变换，而对于HSV来说，增加亮度只需要在V这一维做变换即可。这两种方法得出来的结果并不相同，效果图如下所示。代码如下：srcImage=imread('lula.jpg');subplot(221)imshow(srcImage)title("原图")hsvImage=rgb2hsv(srcImage);%转hsvhsvImage(:,:,3)=hsvImage(:,:,3).*0.6;%变暗0.6darkImage=hsv2rgb(hsvIma

原创 matlab彩色图像变暗

左侧是原图，右侧是变暗图像主要是将图片从RGB域转换为HSV域，然后将亮度V分量降低，最后代码如下：

原创 latex 二重积分、三重积分、闭合积分公式

举个例子，散度公式：∯SA⋅dS=∭V∇⋅AdV\oiint \limits_{S} \mathbf{A} \cdot \mathrm{d}\mathbf{S} = \iiint \limits_{V} \nabla\cdot \mathbf{A} \mathrm{d}VS∬​​A⋅dS=V∭​∇⋅AdV代码：\begin{equation}\oiint \limits_{S} \mathbf{A} \cdot \mathrm{d}\mathbf{S} = \iiint \limits_{V}

原创 MATLAB simulink 2FSK调制与解调实验，附仿真文件（西电B测）

simulink仿真文件链接：（待审核通过后发）整个报告是用LaTeX写的，摘要目录总结参考文献什么的就不放了，直接截取核心部分。如有需要可以私信给完整文件。如有不足之处望指出，欢迎交流！温馨提示：文章是用来分享技术实现，提高知识水平的，而不是直接抄袭过B测的。正文开始：要求及说明任务要求任务说明方案设计及论证数字带通传输系统介绍2FSK基本原理总体设计思路结果与分析初始信号的设定调制信号仿真解调的仿真仿真文件由zhb完

原创 MATLAB 相控阵雷达原理及仿真（雷达原理报告）

一项新的雷达原理作业，依旧是拿word写的，懒得转LaTeX了，截图凑合看吧，需要word的私信，有问题欢迎评论区留言。一、 问题提出二、问题解决1.相控阵雷达原理：2.天线图仿真实验：3.实际信号仿真实验：三、心得感悟附代码：theta=45*pi/180;%波束角度，以度算d=10;%阵列间距,mN=16;%阵元总数lambda=10;%信号波长,mc=3e8;%光速f=c/lambda;%信号频率dR=d.*sin(theta);%相邻单元路程

原创 MATLAB LFCM雷达调频法测距

雷达原理老师留的作业，干了一整天，收获蛮大，拿word写的报告，夜已经深了（21点），懒得转成LaTeX了，就直接截图放上来了。一、问题提出LFCM（线性调频连续波）雷达。发射信号带宽1GHz,时宽10us. 3个目标径向距离分别为100m,105m,105.1m。分析中频信号带宽。仿真产生回波信号及去斜处理过程，对中频输出做傅里叶变换，观察距离维回波特性并进行讨论分析。在此基础上，讨论该类雷达性能指标和信号参数之间的关系。二、问题解决1.线性调频信号（LFM）原理：2.线性调频信号（LFM）产

原创 Matlab使用fft画出信号频谱（详细注释）

做雷原作业，需要对信号进行频谱分析，网上一搜太乱了，很多不是我想要的，特此整理。DSP还没学到fft，就不对fft原理进行详细解释了，直接上代码。fs=500;%采样率f1=5;%信号频率f2=10;%信号频率T=1;%时宽1sn=round(T*fs);%采样点个数t=linspace(0,T,n);%时域横坐标x = 3+cos(2*pi*f1*t) + 2.*cos(2*pi*f2*t);%形成三频信号,注意第二个频率信号幅度为2，直流幅度为3figure(1);plot(t,x)

原创 matlab 图片频域去噪实例（详细注释）

I=imread('lula.jpg');%输入lulaimshow(I)I=imnoise(I, 'gaussian', 0, 0.01);%加高斯白噪声，均值为0，方差为0.01figure(2)imshow(I)J=fftshift(fft2(I));%快速傅里叶变换加shiftfigure(3)%归一化% M=real((J))-min(min(min(real(J))));%减去最小的，搬到0% M=M./max(max(max(M)));%除以最大的，压缩到0-1% M=M.

原创 DSP matlab产生正弦数字采样信号（M2.4）

一些基础知识：模拟信号被采样形成数字信号，假设采样频率是fsf_sfs​，则采样间隔是T=1/fsT=1/f_sT=1/fs​，设nnn是采样序号，则连续时间t=nTt=nTt=nT，其中nnn是正整数。一个连续正弦信号的标准形式如下：Acos⁡(ωt+φ)A\cos(\omega t+\varphi)Acos(ωt+φ)代入t=nTt=nTt=nT有Acos⁡(ωnT+φ)A\cos(\omega nT+\varphi)Acos(ωnT+φ)定义数字角频率为Ω=ωT\Omega = \omega

原创 DSP sawtooth锯齿波与square方波matlab产生（M2.2）

x = sawtooth(t,xmax)产生锯齿波序列，有两个参数，其中第二个参数xmax可省略。该函数周期为2π2\pi2π，t是时间刻度序列，xmax是刻度伸缩系数，介于0到1之间，默认为1，默认幅度从-1到+1锯齿上升。x = square(t,duty)与sawtooth函数类似，周期为2π2\pi2π，t是时间刻度序列，duty是占空比。实际使用代码如下：A = input('The peak value =');%峰值7L = input('Length of sequence =')

原创 LaTeX tikz初探——空间矢量旋转示意图，四元数（4）

zhb学了一个晚上，又完成了一项作业，提供了代码。先放图：代码和相关注释笔记如下：在这里插入代码片

原创 LaTeX tikz初探——利用emoji画GPS卫星3D分布图（3）

zhb学了一个晚上，又提供了3D卫星图的绘制。效果图：笔记写在注释里了% HB-吃井不忘挖水人\documentclass[tikz,border=3.14mm]{standalone}\usepackage{tikz-3dplot}\usepackage{emoji}\usetikzlibrary{patterns,backgrounds,math}\begin{document} % 调整视角，绕x轴顺时针70，绕z顺时针110，默认是从z正方向向下看，x轴指向右侧 \tdplot

原创 LaTeX tikz初探——利用emoji画GPS卫星2D分布图（2）

zhb提供的代码，我拿来啃了一段，写了些笔记注释。emoji包需用LaTeX最新版本，实测texlive2018不行，texlive2020可以。编译系统需要选择：LuaLatex效果图：笔记写注释里了\documentclass[tikz,border=2pt]{standalone} \usepackage{ctex}\usepackage{emoji}\begin{document} \begin{tikzpicture}\usetikzlibrary{intersections

原创 python AES对称加密文件、解密文件

咱也不是学计算机的，咱也不是学网络安全的，咱更不是学密码学的，所以东拼西凑一堆代码，能用就行。该加解密的秘钥是自己输入的密码加电脑固定序列号，包括网卡MAC地址、CPU序列号、硬盘序列号、主板序列号、BIOS序列号。这样可以保证加密文件只能在加密时的计算机打开，在其他人的计算机无法打开，不想用这个功能的可以在生成秘钥那里把序列号删掉。序列号的获取参考我写的这篇文章：https://blog.csdn.net/weixin_42845306/article/details/114899172我下面的

原创 LaTeX tikz初探——基本图形绘制（1）

导师让我们写的GPS作业，点名要用LaTeX里的tikz库，又是需要自学的东西。LaTeX和Tikz 的安装就略过了，网上方法一大堆。照着用户手册，直奔主题。首先，试试它给的模板：\documentclass{article} % say \usepackage{tikz} \begin{document} We are working on \begin{tikzpicture} \draw (-1.5,0) -- (1.5,0); \draw (0,-1.5) -- (0,1.5);