先看QPSK的。
用C语言计算出模拟的QPSK解调信号I路和Q路数据,然后以gnuplot作图。C代码演示了如何通过命令行参数输入噪声大小及解调参考相位。同时也给出了标准正态分布随机数的计算函数。
C代码文件名为:QPSKconst.c
编译:
gcc QPSKconst.c
生成a.exe
执行时输入噪声大小系数,如0.1,解调相位角(度数),如15:
a.exe 0.1 15
则从屏幕输出1000行星座图信号坐标点。在gnuplot中作图即可。作图环境可用以下设置。
set size square
set grid
unset key
然后,在gnuplot命令窗中输入作图命令
gnuplot> plot [-2:2] [-2:2] "<a.exe 0.3 0" w p pt 6 lc 3
"<a.exe 0.3 0"
是执行带参数的exe文件并将结果重定向输入到plot命令中。
得:
又,减小噪声,得
gnuplot> plot [-2:2] [-2:2] "<a.exe 0.05 0" w p pt 6 lc 3
)
又,相位偏移-10
度,得
plot [-2:2] [-2:2] "<a.exe 0.05 -10" w p pt 6 lc 3
)
噪声很大的情况,星座图点散开。
plot [-2:2] [-2:2] "<a.exe 1 0" w p pt 6 lc 3
)
附:C代码
//QPSK信号星座图QPSKconst.c
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define PI 3.14159265
double randn()//标准高斯噪声产生(0均值,方差1)
{
double r1,r2;
r1=(double)rand()/RAND_MAX;
r2=(double)rand()/RAND_MAX;
return sqrt(-2*log(r1+1e-100))*cos(2*PI*r2);//1e-100防止溢出
}
main(int argc, char *argv[])
{
double x,y,x1,y1,s=45;//s相位旋转角(度)
int i;
double a=0.05;
srand(1234);//随机数种子
if(argc!=3)
{
//如果输入参数不足或多了则按默认参数计算
printf("#Usage: a.exe att angle\n");
printf("#Default: a.exe 0.05 45\n");
}
else
{
a=atof(argv[1]);
s=atof(argv[2]);
}
s=s/180.0*PI;//角度制转弧度
for(i=0; i<1000; i++)
{
//标准QPSK解调信号
x=((double)rand()/RAND_MAX>0.5)? 1:-1;
y=((double)rand()/RAND_MAX>0.5)? 1:-1;
x=x+a*randn();//加复高斯噪声
y=y+a*randn();
x1=x*cos(s)-y*sin(s);//相位旋转
y1=x*sin(s)+y*cos(s);
printf("%f\t%f\n",x1,y1);//输出星座图数据
}
}
以此类似,可得16QAM、64QAM、BPSK、8PSK的星座图程序。
//QPSK信号星座图QAMconst.c
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define PI 3.14159265
double randn()//标准高斯噪声产生(0均值,方差1)
{
double r1,r2;
r1=(double)rand()/RAND_MAX;
r2=(double)rand()/RAND_MAX;
return sqrt(-2*log(r1+1e-100))*cos(2*PI*r2);
}
main(int argc, char *argv[])
{
double x,y,x1,y1,s=0;//s相位旋转角(度)
int i;
int M=sqrt(64);
double a=0.05;
srand(1234);//随机数种子
if(argc!=4)
{
//如果输入参数不足或多了则按默认参数计算
printf("#Usage: a.exe att angle\n");
printf("#Default: 64QAM a.exe 0.05 0 64\n");
}
else
{
a=atof(argv[1]);
s=atof(argv[2]);
M=sqrt(atoi(argv[3]));
}
s=s/180.0*PI;//角度制转弧度
for(i=0; i<1000; i++)
{
//标准QAM解调信号
x=(rand()%M)*2-M+1;
y=(rand()%M)*2-M+1;
x=x+a*randn();//加复高斯噪声
y=y+a*randn();
x1=x*cos(s)-y*sin(s);//相位旋转
y1=x*sin(s)+y*cos(s);
printf("%f\t%f\n",x1,y1);//输出星座图数据
}
}
编译作图64QAM:
gnuplot> plot [-12:12] [-12:12] "<a.exe 0.1 0 64" w p pt 6 lc 3
gnuplot> plot [-12:12][-12:12] "<a.exe 0.1 0 16" w p pt 6 lc 3
M-PSK
//M-PSK信号星座图MPSKconst.c
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define PI 3.14159265
double randn()//标准高斯噪声产生(0均值,方差1)
{
double r1,r2;
r1=(double)rand()/RAND_MAX;
r2=(double)rand()/RAND_MAX;
return sqrt(-2*log(r1+1e-100))*cos(2*PI*r2);
}
main(int argc, char *argv[])
{
double x,y,x1,y1,s=0;//s相位旋转角(度)
int i;
int M=8;
double a=0.05;
srand(12345678);//随机数种子
if(argc!=4)
{
//如果输入参数不足或多了则按默认参数计算
printf("#Usage: a.exe att angle\n");
printf("#Default: 8PSK a.exe 0.05 0 8\n");
}
else
{
a=atof(argv[1]);
s=atof(argv[2]);
M=atoi(argv[3]);
}
s=s/180.0*PI;//角度制转弧度
for(i=0; i<1000; i++)
{
//标准MPSK解调信号
x=(double)(rand()%M)/M*2*PI;
y=10*sin(x);
x=10*cos(x);
x=x+a*randn();//加复高斯噪声
y=y+a*randn();
x1=x*cos(s)-y*sin(s);//相位旋转
y1=x*sin(s)+y*cos(s);
printf("%f\t%f\n",x1,y1);//输出星座图数据
}
}
编译作图:
plot [-12:12][-12:12] "<a.exe 0.1 0 2" w p pt 6 lc 3
)
8PSK
plot [-12:12][-12:12] "<a.exe 0.1 0 8" w p pt 6 lc 3
文章浏览阅读1k次。1、HTML5 Web SQL 数据库 Web SQL 数据库 API 并不是 HTML5 规范的一部分,但是它是一个独立的规范,引入了一组使用 SQL 操作客户端数据库的 APIs。如果你是一个 Web 后端程序员,应该很容易理解 SQL 的操作。Web SQL 数据库可以在最新版的 Safari, Chrome 和 Opera 浏览器中工作。2、核心方法 以下是规范中定义的三个_方式准则的定义
文章浏览阅读4.1k次,点赞2次,收藏6次。spring Boot 中使用线程池异步执行多个定时任务在启动类中添加注解@EnableScheduling配置自定义线程池在启动类中添加注解@EnableScheduling第一步添加注解,这样才会使定时任务启动配置自定义线程池@Configurationpublic class ScheduleConfiguration implements SchedulingConfigurer..._springboot启动后自动开启多个线程程序
文章浏览阅读1.1k次。在项目的target文件夹下把之前"mvn clean package"生成的压缩包(我的是jar包)删掉重新执行"mvn clean package"再执行"mvn clean install"即可_mvn clean install有errors
文章浏览阅读974次。Navicat连接mysql数据库时,不断报1405错误,下面是针对这个的解决办法:MySQL服务器正在运行,停止它。如果是作为Windows服务运行的服务器,进入计算机管理--->服务和应用程序------>服务。如果服务器不是作为服务而运行的,可能需要使用任务管理器来强制停止它。创建1个文本文件(此处命名为mysql-init.txt),并将下述命令置于单一行中:SET PASSW..._nvarchar链接不上数据库
文章浏览阅读2.2k次。Python的requests模块是一个常用的HTTP库,用于发送HTTP请求和处理响应。_python requests 参数
文章浏览阅读2.7w次,点赞7次,收藏50次。APT攻击APT攻击是近几年来出现的一种高级攻击,具有难检测、持续时间长和攻击目标明确等特征。本文中,整理了近年来比较典型的几个APT攻击,并其攻击过程做了分析(为了加深自己对APT攻击的理解和学习)Google极光攻击2010年的Google Aurora(极光)攻击是一个十分著名的APT攻击。Google的一名雇员点击即时消息中的一条恶意链接,引发了一系列事件导致这个搜_2010谷歌网络被极光黑客攻击
文章浏览阅读1.1k次。JS代码 /** * 生命周期函数--监听页面加载 */ onLoad: function (options) { setTimeout( function(){ wx.showToast({ title: '黄菊华老师', }) },2000 ) },说明该代码只执行一次..._微信小程序 settimeout 向上层传值
文章浏览阅读48次。uploadify2.1.4如何能使按钮显示中文博客分类:uploadify网上关于这段话的搜索恐怕是太多了。方法多也试过了不知怎么,反正不行。最终自己想办法给解决了。当然首先还是要有fla源码。直接去管网就可以下载。[url]http://www.uploadify.com/wp-content/uploads/uploadify-v2.1.4...
文章浏览阅读9.6k次,点赞5次,收藏36次。戴尔服务器安装VMware ESXI6.7.0教程(U盘安装)一、前期准备1、下载镜像下载esxi6.7镜像:VMware-VMvisor-Installer-6.7.0-8169922.x86_64.iso这里推荐到戴尔官网下载,Baidu搜索“戴尔驱动下载”,选择进入官网,根据提示输入服务器型号搜索适用于该型号服务器的所有驱动下一步选择具体类型的驱动选择一项下载即可待下载完成后打开软碟通(UItraISO),在“文件”选项中打开刚才下载好的镜像文件然后选择启动_vmware-vcsa-all-6.7.0-8169922.iso
文章浏览阅读2k次。百度语音技术永久免费的语音自动转字幕介绍基于百度语音技术,识别率97%无时长限制,无文件大小限制永久免费,简单,易用,速度快支持中文,英文,粤语永久免费的语音转字幕网站: http://thinktothings.com视频介绍 https://www.bilibili.com/video/av42750807 ...
文章浏览阅读7.6k次,点赞2次,收藏9次。Instrumentation是一种直接修改程序二进制文件的方法。其可以用于程序的调试,优化,安全等等。对这个词一般的翻译是“插桩”,但这更多使用于软件测试领域。【找一些相关的例子】Dyninst可以动态或静态的修改程序的二进制代码。动态修改是在目标进程运行时插入代码(dynamic binary instrumentation)。静态修改则是直接向二进制文件插入代码(static b_dyninst
文章浏览阅读2.9k次。部署asp网站到云服务器 内容精选换一换通常情况下,需要结合客户的实际业务环境和具体需求进行业务改造评估,建议您进行服务咨询。这里仅描述一些通用的策略供您参考,主要分如下几方面进行考虑:业务迁移不管您的业务是否已经上线华为云,业务迁移的策略是一致的。建议您将时延敏感型,有快速批量就近部署需求的业务迁移至IEC;保留数据量大,且需要长期稳定运行的业务在中心云上。迁移方法请参见如何计算隔离独享计算资源..._nas asp网站