jian_an_na的博客

先记录代码，代码的具体含义以后再细理解。

在仿真时，经常需要从文件夹里载入变量。在网上搜了一小会，没有结果，特此记录，方便以后翻阅。代码如下：a = importdata('filename\datainsyn_i.txt');很久没有写博客，也很久没有停下来整理知识点了。再过一年就得找工作了，慌。...

在进行FFT、IFFT时分别需除sqrt(N),乘sqrt(N），如下：Data = randi([0,1],100,1)+i*randi([0,1],100,1);ifftData = ifft(Data)*sqrt(100);fftData = fft(Data)./sqrt(100);这样做的目的是使能量归一化，即，使得时域和频域数据的能量一样。测试如下：Data_inf = randi([0,1],100,1)+i*randi([0,1],100,1);ifftData = iff

还是上次的小活（求最高频偏和频偏序列长度的关系）。刚给老师看了运行结果，老师很快就指出了我求频偏的那句代码存在错误，并且给我讲解了求频偏的原理。记录下来。（我真是只小菜鸡qaq，老师威武因此，在估计频偏时，应当对前后两段序列中的一段取共轭，之后在点乘。实现代码如下：% x_noise是列向量cfo_est = angle(x_noise(nn).'*conj(x_noise(nn+Nlength(i))))*360/(2*pi);% 采用下面的代码也一样的cfo_est = angle(x_

最近做一个老师给的小活，有一些代码上的易错点或者小收获，记录下来。MATLAB里共轭转置与转置a = [1+1j,2+2j];a'= [1-1j;2-2j]; % ' 表示求共轭转置a.'= [1+1j;2+2j]; % ' 表示求转置（这个问题已经错了两次了，事不过三啊qaq因此，再求信号功率时，可以用如下代码：sigPow = sig*sig'; % sig是行向量而不是sigPow = sig*sig.'; % sig是行向量lengend 中加入运行产生的数

在将代码由软件仿真转向硬件实现时，有一个特别需要注意的点：硬件里没有浮点数。而在硬件实现中，如果数据是小数，经常需要扩位来使其整数化。因此，有一个很重要的问题：在MATLAB如何确定一个向量或矩阵中是否有小数。下面给出一种方法：%% 使用rem求余函数find(rem(a,1)~=0); % a是待确定矩阵若a无小数，则返回一个空向量；反之，则返回小数的位置。...

randintrandint(m,n) 随机产生取值0/1的size为（m,n)的矩阵randint(m,n,[a b]） 随机产生取值a至b的size为（m,n)的矩阵randint(m,n,b) 随机产生取值0至b的size为（m,n)的矩阵randirandi([a,b],m,n) 随机产生取值a至b的size为（ m,n)的矩阵randi(b,m,n)