m基于matlab的TDSCDMA系统性能仿真

news2025/2/26 3:43:14

目录

1.算法概述

2.仿真效果预览

3.核心MATLAB代码预览

4.完整MATLAB程序

1.算法概述

        TD-SCDMA的中文含义为时分复用同步码分多址接入，是由中国第一次提出、在无线传输技术(RTT)的基础上完成并已正式成为被ITU接纳的国际移动通信标准。这是中国移动通信界的一次创举和对国际移动通信行业的贡献，也是中国在移动通信领域取得的前所未有的突破 。
        TD-SCDMA中的TD指时分复用，也就是指在TD-SCDMA系统中单用户在同一时刻双向通信(收发)的方式是TDD(时分双工)，在相同的频带内在时域上划分不同的时段(时隙)给上、下行进行双工通信，可以方便地实现上、下行链路间的灵活切换。例如根据不同的业务对上、下行资源需求的不同来确定上、下行链路间的时隙分配转换点，进而实现高效率地承载所有3G对称和非对称业务。与FDD模式相比，TDD可以运行在不成对的射频频谱上，因此在当前复杂的频谱分配情况下它具有非常大的优势。TD-SCDMA通过最佳自适应资源的分配和最佳频谱效率，可支持速率从8kb/s到2Mb/s以及更高速率的语音、视频电话、互联网等各种3G业务 

       TD-SCDMA的发展始于1998年初，当时在国家邮电部的直接领导下，由原电信科学技术研究院组织队伍在 SCDMA技术的基础上，研究和起草符合IMT-2000要求的TDSCDMA建议草案。该标准草案以智能天线、同步码分多址、接力切换、时分双工为主要特点，于ITU征集IMT-2000第三代移动通信无线传输技术候选方案的截止日1998年6月30日提交到ITU，从而成为IMT2000的15个候选方案之一。ITU综合了各评估组的评估结果。在1999年11月举行的赫尔辛基ITU-RTG8/1第18次会议上和2000年5月举行的伊斯坦布尔ITU-R全会上，TD-SCDMA被正式接纳为CDMATDD制式的方案之中国无线通信标准研究组(CWTS)作为代表中国的区域性标准化组织，自1999年5月加入3GPP后，经过4个月的充分准备，与项目协调组(3 GPPPCG)、技术规范组(TSG)进行了大量协调工作，在同年9月向3GPP建议将TD- SCDMA纳入3GPP标准规范的工作内容。1999年12月在法国尼斯举行的3GPP会议上，提案被无线接入网(3 GPPTSGRAN)全会所接受，正式确定将TD- SCDMA纳入 Release200(后拆分为R4和R5)的工作计划中，并将 TD-SCDMA简称为即低码片速率TDD方案(Low Code rate, LCRTDD) [4]  。

在TD-SCDMA系统中，用到了以下几种主要关键技术 ：
(1)时分双工方式(Time Division Duplexing)；
(2)联合检测(Joint Detection)；
(3)智能天线(Smart Antenna)；
(4)上行同步(Uplink Synchronous)；
(5)软件无线电(Soft Radio)；
(6)动态信道分配(Dynamic Channel Allocation)；
(7)功率控制(Power control)；
(8)接力切换(Baton Handover)；
(9)高速下行分组接入技术(High Speed Downlink Packet Access)  。
 

2.仿真效果预览

matlab2022a仿真测试如下：

 

 

3.核心MATLAB代码预览

snr_indb=1:8;
for k=1:length(snr_indb)
snr=10^(snr_indb(k)/10);
sgma=1;
eb=2*(sgma^2)*snr;
LC=31;
echip=eb/LC;
N=1000;%number of bits  transmitted
%creat PN codes
fbconnection=[0 1 0 0 1]; 
mseq=m_sequence(fbconnection); 
fbconnection1=[0 0 1 0 1];
fbconnection2=[0 1 1 1 1]; 
goldseq=gold_seq(fbconnection1,fbconnection2); 
%N=2^length(fbconnection)-1; 
ind1=find(mseq==0);
mseq(ind1)=-1;
ind2=find(goldseq==0);
goldseq(ind2)=-1;            %creat 31*31 gold sequence
temp=goldseq;
pn_seq1=temp(3,:);
pn_seq2=temp(4,:);
pn_seq3=temp(7,:);
pn_seq4=temp(10,:);
pn_seq5=temp(15,:);
pn_seq6=temp(20,:);
pn_seq7=temp(26,:);
pn_seq8=temp(31,:);
% comput  matrix R
pp=[pn_seq1;pn_seq2;pn_seq3;pn_seq4;pn_seq5;pn_seq6;pn_seq7;pn_seq8];
RR=(1/LC)*pp*pp';
RR1=inv(eb*eye(8));
R=inv(RR+RR1);

%generate  codes

for  i=1:N
temp1=rand;
if(temp1<0.5),dsource1(i)=-1;
else   dsource1(i)=1;
end;
temp2=rand;
if(temp2<0.5),dsource2(i)=-1;
else   dsource2(i)=1;
end;
temp3=rand;
if(temp3<0.5),dsource3(i)=-1;
else   dsource3(i)=1;
end;
temp4=rand;
if(temp4<0.5),dsource4(i)=-1;
else   dsource4(i)=1;
end;
temp5=rand;
if(temp5<0.5),dsource5(i)=-1;
else   dsource5(i)=1;
end;
temp6=rand;
if(temp6<0.5),dsource6(i)=-1;
else   dsource6(i)=1;
end;
temp7=rand;
if(temp7<0.5),dsource7(i)=-1;
else   dsource7(i)=1;
end;
temp8=rand;
if(temp8<0.5),dsource8(i)=-1;
else   dsource8(i)=1;
end;
end;

number_of_err=0;
number_of_err1=0;
number_of_err2=0;

for  i=1:N
%将每个bit repeat LC=31 times
for  j=1:LC
    repeatdata1(j)=dsource1(i);
    repeatdata2(j)=dsource2(i);
    repeatdata3(j)=dsource3(i);
    repeatdata4(j)=dsource4(i);
    repeatdata5(j)=dsource5(i);
    repeatdata6(j)=dsource6(i);
    repeatdata7(j)=dsource7(i);
    repeatdata8(j)=dsource8(i);
end;
 for ii=0:2^8-1
            j=1;aa=ii;
            while aa~=0
                bk(j,ii+1)=rem(aa,2);
                aa=floor(aa/2);
                j=j+1;
            end
 end
 bk=bk*2-1;
%tranmit signal is:

trans_sig1=sqrt(echip)* repeatdata1.*pn_seq1;
trans_sig2=sqrt(echip)* repeatdata2.*pn_seq2;
trans_sig3=sqrt(echip)* repeatdata3.*pn_seq3;
trans_sig4=sqrt(echip)* repeatdata4.*pn_seq4;
trans_sig5=sqrt(echip)* repeatdata5.*pn_seq5;
trans_sig6=sqrt(echip)* repeatdata6.*pn_seq6;
trans_sig7=sqrt(echip)* repeatdata7.*pn_seq7;
trans_sig8=sqrt(echip)* repeatdata8.*pn_seq8;

%add AWGN  noise
noise=sgma*randn(1,LC);
 % receive signal
 for  j=1:LC
     rtemp(j)=trans_sig1(j)+trans_sig2(j)+trans_sig3(j)+trans_sig4(j)+trans_sig5(j)+trans_sig6(j)+trans_sig7(j)+trans_sig8(j);
 end;
 r=rtemp+noise;
r1=r;r2=r;r3=r;r4=r;r5=r;r6=r;r7=r;r8=r;

%CD JUDGE
I1=sum(r1.*pn_seq1);
I2=sum(r2.*pn_seq2);
I3=sum(r3.*pn_seq3);
I4=sum(r4.*pn_seq4);
I5=sum(r5.*pn_seq5);
I6=sum(r6.*pn_seq6);
I7=sum(r7.*pn_seq7);
I8=sum(r8.*pn_seq8);

% DD JUDGE
I=R*[I1,I2,I3,I4,I5,I6,I7,I8]';
y=[I1,I2,I3,I4,I5,I6,I7,I8]';

for ii=1:2^8
  
    c(ii)=2*(bk(:,ii))'*y-(bk(:,ii))'*R*bk(:,ii);
       
end
[m,mm]=max(c);
bb=bk(:,mm);
if(bb(1)~=dsource1(i)),number_of_err2=number_of_err2+1;end;
if(bb(2)~=dsource2(i)),number_of_err2=number_of_err2+1;end;
if(bb(3)~=dsource3(i)),number_of_err2=number_of_err2+1;end;
if(bb(4)~=dsource4(i)),number_of_err2=number_of_err2+1;end;
if(bb(5)~=dsource5(i)),number_of_err2=number_of_err2+1;end;
if(bb(6)~=dsource6(i)),number_of_err2=number_of_err2+1;end;
if(bb(7)~=dsource7(i)),number_of_err2=number_of_err2+1;end;
if(bb(8)~=dsource8(i)),number_of_err2=number_of_err2+1;end;

% DD make decision
if(I(1)<0), desion1=-1;
else desion1=1;
end;
if(I(2)<0), desion2=-1;
else desion2=1;
end;
if(I(3)<0), desion3=-1;
else desion3=1;
end;
if(I(4)<0), desion4=-1;
else desion4=1;
end;
if(I(5)<0), desion5=-1;
else desion5=1;
end;
if(I(6)<0), desion6=-1;
else desion6=1;
end;
if(I(7)<0), desion7=-1;
else desion7=1;
end;
if(I(8)<0), desion8=-1;
else desion8=1;
end;

if(desion1~=dsource1(i)),number_of_err=number_of_err+1;end;
if(desion2~=dsource2(i)),number_of_err=number_of_err+1;end;
if(desion3~=dsource3(i)),number_of_err=number_of_err+1;end;
if(desion4~=dsource4(i)),number_of_err=number_of_err+1;end;
if(desion5~=dsource5(i)),number_of_err=number_of_err+1;end;
if(desion6~=dsource6(i)),number_of_err=number_of_err+1;end;
if(desion7~=dsource7(i)),number_of_err=number_of_err+1;end;
if(desion8~=dsource8(i)),number_of_err=number_of_err+1;end;



%CD  MAKE DECISION
% CD  make  decision
if(I1<0), desion11=-1;
else desion11=1;
end;
if(I2<0), desion21=-1;
else desion21=1;
end;
if(I3<0), desion31=-1;
else desion31=1;
end;
if(I4<0), desion41=-1;
else desion41=1;
end;
if(I5<0), desion51=-1;
else desion51=1;
end;
if(I6<0), desion61=-1;
else desion61=1;
end;
if(I7<0), desion71=-1;
else desion71=1;
end;
if(I8<0), desion81=-1;
else desion81=1;
end;

if(desion11~=dsource1(i)),number_of_err1=number_of_err1+1;end;
if(desion21~=dsource2(i)),number_of_err1=number_of_err1+1;end;
if(desion31~=dsource3(i)),number_of_err1=number_of_err1+1;end;
if(desion41~=dsource4(i)),number_of_err1=number_of_err1+1;end;
if(desion51~=dsource5(i)),number_of_err1=number_of_err1+1;end;
if(desion61~=dsource6(i)),number_of_err1=number_of_err1+1;end;
if(desion71~=dsource7(i)),number_of_err1=number_of_err1+1;end;
if(desion81~=dsource8(i)),number_of_err1=number_of_err1+1;end;
end
err(k)=number_of_err/(3*N)%MMSE ber
err1(k)=number_of_err1/(3*N)%CD ber
err2(k)=number_of_err2/(3*N)
end
%plot(snr_indb,err);
semilogy(snr_indb,err,'k');%MMSE ber
grid on;
hold on;
semilogy(snr_indb,err1,'r');%CD ber
hold on
semilogy(snr_indb,err2,'^-')

legend('MMSE','CD','不加信道估计');
01_006_m

4.完整MATLAB程序

matlab源码说明_我爱C编程的博客-CSDN博客

V

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/9960.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

解读ESSumm: Extractive Speech Summarization from Untranscribed Meeting

解读ESSumm: Extractive Speech Summarization from Untranscribed Meeting

ESSumm: Extractive Speech Summarization from Untranscribed Meeting 从非转录会议中提取语音摘要&#xff08;非转录&#xff1a;原音频&#xff0c;没有转化为文本或者其他格式&#xff09; 论文地址 https://arxiv.org/abs/2209.06913 摘要 在本文中&#xff0c;我们为…
阅读更多...
JS中的迭代器、可迭代对象、生成器

JS中的迭代器、可迭代对象、生成器

迭代器 迭代器&#xff08;iterator&#xff09;&#xff0c;是使用户在容器对象&#xff08;container&#xff0c;例如链表或数组&#xff09;上可以遍历访问的对象&#xff0c;使用该接口无需关心容器对象的内部实现细节。 从上面定义上来看&#xff0c;迭代器是帮助我们对…
阅读更多...
SpringMVC 源码学习 返回值处理

SpringMVC 源码学习 返回值处理

SpringMVC中对返回值的数据基本分为两类&#xff1a; 1、响应数据 2、响应页面 一、响应数据 响应数据大多数都是将返回值的格式转换为JSON格式然后展示在页面或者保存i起来。 第一步&#xff1a;在SpringBoot中需要引入json场景 <dependency><groupId>org.sprin…
阅读更多...
如何设计存储架构

如何设计存储架构

步骤 步骤1&#xff1a;估算性能需求 任务 基于具体的业务场景来估算性能需求&#xff0c;包括存储量、读写性能等 挑战 不知道如何估算担心估算不准 步骤2&#xff1a;选择存储系统 任务 根据技术储备、方案优缺点选择合适的存储系统 挑战 不知道有哪些存储系统知道…
阅读更多...
UG NX二次开发（C#）-UI Styler-批量选择点

UG NX二次开发（C#）-UI Styler-批量选择点

1、前言 在设计UG NX的对话框界面时，有时需要选择点，在UI Styler编辑器中已经设置好了可以指定点选择按钮，这个点选择对话框是单选，不是多选的，如果想选择多个点，那么可以采用对象选择按钮，本文介绍下选择点的操作。 2、指定点按钮 2.1 设计UI Styler UG NX的UI Sty…
阅读更多...
vscode中Emmet语法的使用

vscode中Emmet语法的使用

每篇博文的浪漫主义 【镰仓旅拍|落日绝景 极致画质|沉浸式旅行 FX3】 https://www.bilibili.com/video/BV1jg411Y7vC/?share_sourcecopy_web&vd_source385ba0043075be7c24c4aeb4aaa73352 镰仓旅拍|落日绝景 极致画质|沉浸式旅行 FX31.1快速生成HTML结构语法 生成标签直接…
阅读更多...
想要彻底卸载Mac应用程序，还得要用这些方法才行

想要彻底卸载Mac应用程序，还得要用这些方法才行

Mac电脑如果有太多无用的应用程序&#xff0c;很有可能会拖垮Mac系统的运行速度。因此&#xff0c;卸载电脑中无用的软件是优化Mac系统运行速度的最佳方式之一。Mac删除应用程序特别简单&#xff0c;长点击应用点击x&#xff0c;或是直接将应用拖进废纸篓。但是有一些应用长按没…
阅读更多...
Qt——（详细）“项目在Debug构建环境下能运行而在Release构建环境下不能运行”解决方案之一，以及 禁用（黄色）警告

Qt——（详细）“项目在Debug构建环境下能运行而在Release构建环境下不能运行”解决方案之一，以及 禁用（黄色）警告

系列文章目录 提示&#xff1a; 文章目录系列文章目录前言环境一、问题准备工作——为了在Release环境下可以进行断点调试分析二、解决1、根据需求&#xff0c;对函数类型进行更改2、根据需求&#xff0c;在函数内添加“return [int]”延伸——“禁用警告”消除 变量 的“黄色感…
阅读更多...
现代修谱，如何看待支系单飞的现象？

现代修谱，如何看待支系单飞的现象？

族谱与支谱、房谱的区别 现代修谱&#xff0c;修的是什么谱&#xff0c;你知道吗&#xff1f;其实现代修谱的种类有很多种&#xff0c;有支谱、房谱、族谱、宗谱、统谱、通谱等等&#xff0c;而这些在生活中都被我们简称为家谱。 不过在现代修谱里&#xff0c;宗谱、统谱、通谱…
阅读更多...
极智编程 | 谈谈 C++ 中容器 map 和 unordered_map 的区别

极智编程 | 谈谈 C++ 中容器 map 和 unordered_map 的区别

欢迎关注我的公众号 [极智视界]&#xff0c;获取我的更多笔记分享 大家好&#xff0c;我是极智视界&#xff0c;本文来 谈谈 C 中 map 和 unordered_map 的区别。 map 和 unordered_map 都可以看做是一种 key-value 的映射关系&#xff0c;unordered_map 可以理解为 无序版的ma…
阅读更多...
C语言 数组

C语言 数组

C语言 数组一、一维数组1. 数组的创建方式程序清单1程序清单22. 计算数组的元素的个数3. 数组在内存中的存储方式二、二维数组1. 二维数组的创建方式2. 计算二维数组的行和列3. 二维数组在内存中的存储方式三、数组名的含义总结数组名的应用场景数组名作为函数参数一、一维数组…
阅读更多...
22.11.16打卡 mysql学习笔记

22.11.16打卡 mysql学习笔记

马上要考试了, 越到考试越想玩, 烦躁烦躁烦躁, 没学多少, 争取明天把mysql基础篇学完 DCL介绍 2022年11月16日 14:54 DCL主要用来处理数据库有哪些用户可以访问, 每个用户具有什么样的权限 用户管理 2022年11月16日 15:20 所有用户的数据都存放在系统数据库mysql中的user表…
阅读更多...
Python是什么？要如何学习？

Python是什么？要如何学习？

Python 是荷兰人 Guido van Rossum &#xff08;吉多范罗苏姆&#xff0c;中国程序员称其为“龟叔”&#xff09;在 1990 年初开发的一种解释型编程语言。 Python 的诞生是极具戏曲性的&#xff0c;据 Guido 自述记载&#xff0c;Python 语言是在圣诞节期间为了打发无聊的时间而…
阅读更多...
基于全景相机的视觉里程计算法研究

基于全景相机的视觉里程计算法研究

一、视觉里程计 视觉里程计技术首先建立相机的成像模型&#xff0c;接着通过标定算法计算相机参数&#xff0c;最后建立相邻图像的关联并估计相机运动轨迹。 1.1相机在空间中运动的描述 描述相机在三维空间中的运动状态&#xff0c;即求解相机在空间中不同时刻下的位姿关系。相…
阅读更多...
假冒网站引发多重安全风险 | 官方严正声明：切勿在非官方渠道购买或下载 Navicat 软件

假冒网站引发多重安全风险 | 官方严正声明：切勿在非官方渠道购买或下载 Navicat 软件

Navicat 严正声明 近期&#xff0c;有关于 Navicat 假冒网站的事件&#xff0c;不法分子通过仿制官方网站&#xff0c;诱导用户下载盗版软件。Navicat 官方已正式向国内监管部门举报&#xff0c;提请将该不法网站下架。目前&#xff0c;监管部门已介入调查中。 我司呼吁广大用…
阅读更多...
这次，听人大教授讲讲分布式数据库的多级一致性｜TDSQL 关键技术突破

这次，听人大教授讲讲分布式数据库的多级一致性｜TDSQL 关键技术突破

近年来&#xff0c;凭借高可扩展、高可用等技术特性&#xff0c;分布式数据库正在成为金融行业数字化转型的重要支撑。分布式数据库如何在不同的金融级应用场景下&#xff0c;在确保数据一致性的前提下&#xff0c;同时保障系统的高性能和高可扩展性&#xff0c;是分布式数据库…
阅读更多...
C语言之详解内存操作函数

C语言之详解内存操作函数

个人主页&#xff1a;平行线也会相交 欢迎 点赞&#x1f44d; 收藏✨ 留言✉ 加关注&#x1f493;本文由 平行线也会相交 原创 收录于专栏【C/C】 目录前言memcpy模拟实现memmove模拟实现memcmpmemset前言 memcpy叫做内存拷贝&#xff0c;memmove叫做内存移动&#xff0c;memc…
阅读更多...
【数据结构】模拟实现双向链表

【数据结构】模拟实现双向链表

你必须非常努力&#xff0c;才能显得毫不费劲 目录 1.模拟实现双向链表 1.1 DLinkedList的内部类 1.2 DLinkedList的成员属性 1.3 DLinkedList的成员方法 1.3.1 在链表开头插入一个新结点 1.3.2 在链表结尾插入一个新的结点 1.3.3 计算结点个数 1.3.4 在链表任意位置…
阅读更多...
4.构造器，this，修饰符详解

4.构造器，this，修饰符详解

构造器&#xff1a; 构造器也叫构造方法&#xff0c;无返回值。非构造方法必须要有返回类型 主要作用&#xff1a;完成对象的初始化&#xff0c;创造对象时&#xff0c;自动调用构造器初始化对象 即使没有显示地使用static关键字&#xff0c;构造器实际上也是静态方法 JAVA…
阅读更多...
HTML---基础入门知识详解

HTML---基础入门知识详解

1&#xff1a;标签的概念 在别人写的网页中我们会看到许多文字&#xff0c;图片排版整齐&#xff0c;让人看的赏心悦目&#xff0c;这就是用到了标签&#xff0c;或者说标签就是帮我们实现某种作用的工具&#xff0c;比如制作段落&#xff0c;换行&#xff0c;导入图片&#x…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023