- 博客(30)
- 收藏
- 关注
RSS订阅原创 合成孔径雷达成像——算法与实现图【2.11】与【2.12】仿真——不同β值下的kaiser窗，以及3dB宽度和峰值旁瓣比
先初始化各个参数A = figure();N = 250;len = {};%用来储存标签pslr = [];%用来储存各个pslr值wide_3dB = [];%用来储存各个3dB带宽值其中，{}与[]的区别是：[]用来储存矩阵和向量，每一个[]是一个矩阵或一个向量，{}用来储存cell型数组等，比如字符、矩阵等。若用[]，则最后在输入标签时，会只有一个标签，因为只有一个矩阵，若用{}，则其中会有多个字符，可以正常显示标签。for i = 0:1:6 idx=num2str(i
2021-01-14 04:09:37
2778
1
原创 合成孔径雷达成像——算法与实现图【2.10】仿真——混叠方程可观测频率
书中理论如下：设定fs为400hz，然后利用关系绘制出图。代码如下：fs = 400;F_continue = 0:0.1:2*fs;F_apparent_complex = F_continue-(round(F_continue./fs)).*fs;F_apparent_real = abs(F_apparent_complex);subplot(2,1,1);plot(F_continue,F_apparent_complex);title('(a)复信号');xlabel('可观测
2021-01-08 20:59:51
581
原创 合成孔径雷达成像——算法与实现图【2.9】仿真——复信号采样
T = 1;BW = 300;k = BW/T;% 斜率h = figure();设置初始参数，带宽设置书中给出，h=figure()为了让所有图在一张图中。dt_0 = 1/1200;t_0 = 0:dt_0:T;fs_0 = 1/dt_0;N_0 = length(t_0);p0_0 = t_0+2;p_max_0 = max(p0_0);p_0 = p0_0./p_max_0;f_0 = (0:N_0-1)*fs_0/N_0 - fs_0/2;设置连续的信号，取间
2020-11-25 21:04:15
723
原创 合成孔径雷达成像——算法与实现图【2.8】仿真
T = 1; BW = 150; k = BW/T; h=figure(); N_0 = length(t_0);设置初始参数，带宽在书中给出，为150Hz，h=figure()是为了后面将所有的图在一张图中显示。 p0_0 = t_0+2; p0_max_0 = max(max(p0_0)); p_0 = p0_0./p0_max_0; 信号幅度修正，这里加的2这里的数越大，产生的图像斜率越大。 dt = 1/400;...
2020-11-24 22:06:33
1096
原创 合成孔径雷达成像——算法与实现图【2.7】仿真——奈奎斯特采样率
f0 = 300; fs = 800; Ts = 1/fs;N = 8;先设定参数，这里先做以800Hz为抽样频率的信号，在原图中存在8个抽样点，所以设置抽样点为8个。t1 = linspace( 0, (N-1)*Ts, 501 );ft1 = sin( 2*pi*f0*t1 );设置原信号。t2 = (0:N-1)*Ts;ft2 = sin( 2*pi*f0*t2 );抽取N个点，用800Hz抽样频率选出抽样点。400Hz抽样频率同理可以求得。这里说明，400Hz抽样
2020-11-24 19:17:32
653
原创 合成孔径雷达成像——算法与实现图【2.6】仿真
T = 1;dt = 1/400;t = 0:dt:1;fs = 1/dt;f0 = 15;ft = 100;BW = 75;%带宽75Hzk = BW/T;%斜率为带宽时间比N = length(t);先对确定各个参数，其中带宽书中给出，图的形状为cos(at+bt^2)然后再乘一个削减系数，复数情况为exp(1i*)函数。s1 = cos(2*pi*(f0.*t+0.5*k*t.^2)).*p;s2 = exp(1i*2*pi*((f0-52.5).*t+0.5*k*t.^2)
2020-11-24 10:27:54
1536
原创 合成孔径雷达成像——算法与实现图【2.3】仿真——rect函数、sinc函数的傅里叶变换matlab仿真
sinc函数傅里叶变换后得到rect函数首先先进行初始条件的设定：dt = -200:0.01:200-0.01;fs = 1/dt;dt = dt+(dt==0)*eps;%若有0则选取一个极小值sincx = sin(pi*dt)./(pi*dt);这里确定函数取值范围，取样间隔，取样频率。这里说明一下第三行，它的意义在于若dt=0，则sincx函数没有定义，算出来在矩阵中得出的值是NAN，即计算错误了，所以若dt=0，这行的作用是若值为0，则我们给他加上一个相对于其他值为无穷小量的一个值
2020-11-20 22:53:50
4841
原创 合成孔径雷达成像——算法与实现图【2.2】仿真与解析
代码来源：https://blog.csdn.net/Destiny0321/article/details/53002789?biz_id=102&utm_term=%E5%90%88%E6%88%90%E5%AD%94%E5%BE%84%E9%9B%B7%E8%BE%BE2.2&utm_medium=distribute.pc_search_result.none-task-blog-2allsobaiduweb~default-4-53002789&spm=1018.211
2020-11-19 20:33:56
1891
1
原创 利用python将图片手绘化
首先，我们在这里面会用到的库有PIL的Image与numpy库：from PIL import Imageimport numpy as np然后，我们将要用的图片调出，然后将其转化为灰度图im=Image.open("im.jpg")im=im.convert('L')im.save('灰度.jpg')然后我们要对整个图像取其梯度，然后对图形进行重构a=np.asarray(...
2019-07-06 22:48:43
4096
3
原创 matlab中int函数的用法
我们假设我们要对函数f(x)=x(3+7x+51+x)f(x) = x\left( 3+\frac{7x+5}{1+x} \right)f(x)=x(3+1+x7x+5)从1到2进行积分：syms x y; % 定义符号f(x)=x*(3+(7*x+5)/(1+x)); % 定义被积函数k1=int(f,x,1,2); % 求函数f对符号变量x从1到2的定积分y=double(k1); ...
2019-07-04 19:08:20
88888
2
原创 python学习笔记玫瑰花绘制
画一个玫瑰花：import turtle as t首先定义一个曲线绘制的函数：def DegreeCurve(n,r,d=1): for i in range(n): t.left(d) t.circle(r,abs(d))然后确定初始位置和参数，确定颜色之类的变量：s=0.2# sizet.setup(450*5*s,750*5*s)t.p...
2019-06-04 13:23:23
1744
原创 python学习笔记霍兰德人格分析与错误解决
雷达图霍兰德认为人格兴趣与职业之间有一种内在的对应的关系。我们需要绘制雷达图，并且展示第三方数据。（编译错误AttributeError: Unknown property frac在文章结尾）matplotlib.rcParams[‘font.family’]='SimHei’意为将图片中的字体设置为SimHei，matplotlib.rcParams[‘font.family’]为设置字体...
2019-06-03 21:03:36
3400
6
原创 python学习笔记第八章（体育竞技分析）
我们在这个实例中，我们需要根据对两个人A、B的能力值来预测两人的胜率各是多少。我们假定以下比赛规则：双人击球比赛：1、A、B两个人，回合制的五局三胜。2、开始时一方先发球，直至判分，接下去胜者发球。3、球员只能在发球局得分，15分胜一局。我们处理一个问题时的两种思路有自顶而下和自底而上：自顶向下：（执行）将一个很大的问题表述为许许多多个小问题的组合，然后利用相同的方法来分解小问题，...
2019-05-30 13:11:52
2964
原创 python在获取夏令营信息上的应用小实例
我在填报夏令营申请表的时候，发现其中有一个要统计其中的各个科目成绩分数段的一个选项，而我又想起来这些天学的许多实例，想利用python编出一段可以处理这部分数据的一个程序，应用到的所有的模块均为我之前博客实例中的源代码。填报信息需要的数据如下：如图我们发现我们需要统计所有的成绩在各个分数段有多少科目，我也不知道一共有多少科目，如下是我的思路：1、获取一段不定长的数据我们可以使用计算统计量...
2019-05-29 17:49:37
315
原创 python学习笔记第七章自动绘制图形
自动绘制图形：需求：根据脚本来绘制图形，而不是写代码，是写数据来及进行绘制轨迹。思路：1、定义数据文件的格式2、编写程序，根据文件解析参数绘制图形3、编制数据文件定义格式为自行定义，即数据的顺序，每个数据的作用。在这里我们定义如下：第一个数据：行进距离第二个数据：转向判断，0：左转，1：右转第三个数据：转向角度第四五六个数据：RGB的数值，为在0—1之间的浮点数1、首先是初始...
2019-05-29 15:01:33
1107
原创 python学习笔记第七章文件的使用
一、文件读取整个文件的遍历方法：1、直接使用read函数#遍历文本的方法一fname = input("请输入要打开文件的名称：")fo = open(fname,"r")txt = fo.read()#对全文txt进行处理fo.close()优点为简单、容易，便于理解；缺点是若是读取十分大的文件，比如几个G的文件，会占用十分大的内存去读取。2、按数量处理，逐步读入#遍历文...
2019-05-28 20:58:56
903
原创 python学习笔记第六章文本词频统计
我们这次需要解决的问题是在一篇文章中，哪一些词汇出现的最多，如何去做，我们考虑英文文本和中文的文本。首先，我们先对哈姆雷特的英文文本进行统计词频。1、获取文本并进行归一化def gettext(): txt=open("hamlet.txt","r").read() txt=txt.lower() for ch in '!"#$%&()*+,-./:;<...
2019-05-28 19:57:31
1619
原创 python学习笔记 第六章1
基本统计值的计算1、需要获得多个数据：不知道具体由多少数据，需要定义函数获取。def getnum():#获取不定长的数字 nums=[] inumstr=input("请输入数字（回车退出）：") while inumstr !="": nums.append(eval(inumstr)) inumstr=input("请输入数字（回...
2019-05-25 22:31:39
924
1
原创 python学习笔记 第五章（科赫雪花与科赫曲线）
科赫曲线，也叫做雪花曲线。是一种分形几何，分形几何是一种迭代的几何图形，广泛存在于自然界中。科赫曲线是一种迭代的图形，所以我们是可以利用python中的迭代法来绘制。我们先考虑如何实现迭代。首先我们确定递归的基例，考虑最基本的情况，就是直接绘制一条曲线，没有任何的改变，即0阶的科赫曲线。if n==0:turtle.fd(size)这里的size指的是长度。然后我们来确定递归的...
2019-05-23 17:16:01
13229
7
原创 python学习笔记 第五章2
经典的汉诺塔问题：这里我们可以利用递归的思想去做，递归中重要的三步，我们逐条来实现：1、函数+分支结构2、递归链条3、递归基例函数+分支结构：def hanoi(n,start,end,mid):global countif:else:这里我们可以定义一个函数，里面的参数有：一共有n个圆盘，从start柱子移到end柱子，中间柱子为mid。 这里定义一个全局变量来计算移动...
2019-05-22 22:35:04
133
原创 学python中遇到的问题
1、TypeError: ‘type’ object is not subscriptable类型的定义错误，当时出现的问题如下：文本中写的为for i in range [‘c’,‘C’]，很明显，这里不应该加入range，range后跟数字类型，而[‘c’,‘C’]为列表类型，所以是类型的错误，若出现这类的错误需查找以下情况(a)把类型认错，使用格式的时候使用错误的类型格式。(b)定...
2019-05-19 17:12:14
1661
原创 python学习笔记 第五章1
七段数码管的绘制与拓展：先给出代码：#七段数码管绘制import turtledef drawline(draw):turtle.pendown() if draw else turtle.penup()turtle.fd(40)turtle.right(90)def drawdigit(digit):drawline(True) if digit in [2,3,4,5,6,...
2019-05-19 16:40:36
2577
原创 python学习笔记 第四章2（用蒙特卡罗法计算圆周率）
圆周率的计算：首先用基本公式求解：程序为：#精确圆周率pi=0N=100for k in range(N):pi=pi+1/pow(16,k)(4/(8k+1)-2/(8k+4)-1/(8k+5)-1/(8*k+6))print(“圆周率的值为：{}”.format(pi))运行结果为3.141592653589793下面我们考虑蒙特卡罗方法，思路为随机向一个边长为1...
2019-05-14 20:34:25
1106
1
原创 python学习笔记 第四章1
BMI问题：BMI是对身体质量的刻画，是国际上常用的衡量人体肥胖和健康程度的重要标准，主要用于统计分析。定义为体重（kg）/身高的平方（m）国际和国内的标准如下：代码如下：#BMIheight,weight=eval(input(“请输入身高（米）和体重\（公斤）[逗号隔开]:”))bmi=weight/pow(height,2)print(“BMI数值为：{:.2f}”.fo...
2019-05-12 22:08:06
974
2
原创 python学习笔记 第三章2
time库是python中处理时间的标准库。可以表达计算机时间、提供获取系统时间并格式化输出功能、提供系统级精确计时功能及格式化输出、ti提供精确计时功能用于程序性能分析。时间获取：time() ctime() gmtime()时间格式化：strftime() striptime()程序计时：perf_counter() sleep()进度条实例：（文本进度条）需要根据程序运行不断变化，...
2019-05-11 22:23:01
326
原创 python学习笔记 第三章1
天天向上问题：基本问题：一年365天，每天进步1‰，累计进步多少？一年365天，每天退步1‰，累计剩下多少？其他情况：如果是“双休日退步”呢？双休日需要多努力才能和持续努力一样呢？问题一、二：结果：问题三（双休日）：结果：问题四：考虑A君一年365天每天几步1%，不停歇，B君一年365天，每周工作5天休息2天，休息日下降1%，他需要多努力才和A君一样？大致思路图如下...
2019-05-09 21:00:29
365
原创 python学习笔记——不确定尾数
不确定尾数例子：会发现0.1+0.2不等于0.3，这即为不确定尾数，这不是运行或者编译错误，而是与计算机对数字运算的内部实现原理有关。原理如下：python语言中使用53位二进制来表示浮点类型的小数部分，0.1由0.00011001100110011001100110011001100110011001100110011010来表示，但它不等于0.1，转化为十进制小数时为0.100000...
2019-05-09 01:51:44
2628
1
原创 python学习笔记 第二章
蟒蛇绘制程序：#PythonDraw.pyimport turtleturtle.setup(650,350,200,200)turtle.penup()turtle.fd(-250)turtle.pendown()turtle.pensize(25)turtle.pencolor(“purple”)turtle.seth(-40)for i in range(4):tur...
2019-05-06 22:04:37
248
原创 python学习第一章笔记
温度转化程序：运用软件：Thonny#TempConvert.pyTempStr=input(“请输入带有符号的温度值：”)if TempStr[-1] in [‘F’,‘f’]:C=(eval(TempStr[0:-1])-32)/1.8print(“转换后的温度是{:.2f}C”.format©)elif TempStr[-1] in [‘C’,‘c’]:F=1.8*eval...
2019-05-06 16:37:43
294
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人