亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站! | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? mai_processor.cpp

?? 用matlab程序實現WCDMA系統的仿真
?? CPP
?? 第 1 頁 / 共 2 頁
字號:
	double *NumChanPtr;
	unsigned i,j;
	unsigned SamplesPerInterferer;
	double fee;
	int Offset;
	double OffsetDouble;
	double *RealUplinkSignalPtr,*TempRealUplinkSignalPtr;
	double *ImagUplinkSignalPtr,*TempImagUplinkSignalPtr;
	unsigned N_UplinkSignal;
	unsigned Stop;
	double *RealMpathSignalPtr,*TempRealMpathPtr;
	double *ImagMpathSignalPtr,*TempImagMpathPtr;
	mxArray *mxControlBitsPtr;
	double *InphaseUplinkChipsPtr,*QuadUplinkChipsPtr;
	mxArray *mxBetaCPtr,*mxBetaDPtr;
	double *BetaCPtr,*BetaDPtr;
//	FILE *fp;
//	float sam;

	srand( (unsigned)time( NULL ) );

	SamplesPerInterferer = PulseLength + SamplesPerChip*(CHIPS_PER_FRAME-1);

	//Find the Maximum Delay
	MaxDelay=0;
	TempDelayPtr=DelayPtr;
	for (i = 0; i <NumDelays;i++)
	{
		if (*TempDelayPtr > MaxDelay) MaxDelay = *TempDelayPtr++;
		else TempDelayPtr++;
	}

	RealUplinkSignalPtr = (double *) mxCalloc(NCurrentSignal+MaxDelay,sizeof(double));
	if (RealUplinkSignalPtr == NULL)
		mexErrMsgTxt("\n RealUplinkSignalPtr Array Not Allocated!--exiting \n");
	ImagUplinkSignalPtr = (double *) mxCalloc(NCurrentSignal+MaxDelay,sizeof(double));
	if (ImagUplinkSignalPtr == NULL)
		mexErrMsgTxt("\n ImagUplinkSignalPtr Array Not Allocated!--exiting \n");

	RealMpathSignalPtr = (double *) mxCalloc(NCurrentSignal+MaxDelay,sizeof(double));
	if (RealMpathSignalPtr == NULL)
		mexErrMsgTxt("\n RealMpathSignalPtr Array Not Allocated!--exiting \n");
	ImagMpathSignalPtr = (double *) mxCalloc(NCurrentSignal+MaxDelay,sizeof(double));
	if (ImagMpathSignalPtr == NULL)
		mexErrMsgTxt("\n ImagMpathSignalPtr Array Not Allocated!--exiting \n");

	//Generate Array for Control Bits

	mxNumChanPtr = mxCreateDoubleMatrix(1,1,mxREAL);
	if (mxNumChanPtr == NULL)
		mexErrMsgTxt("\n mxNumChanPtr Matrix Not Allocated!--exiting \n");
	NumChanPtr = mxGetPr(mxNumChanPtr);
	*NumChanPtr = (double) 1.0;

	mxControlBitsPtr = mxCreateDoubleMatrix(1,CONTROL_BITS,mxREAL);
	if (mxControlBitsPtr == NULL)
		mexErrMsgTxt("\n mxControlBitsPtr Matrix Not Allocated!--exiting \n");
	ControlBitsPtr = mxGetPr(mxControlBitsPtr);

	RealCodePtr = (short int *) mxCalloc(CHIPS_PER_FRAME,sizeof(short int));
	if (RealCodePtr == NULL)
		mexErrMsgTxt("\n RealCodePtr Array Not Allocated!--exiting \n");

	ImagCodePtr = (short int *) mxCalloc(CHIPS_PER_FRAME,sizeof(short int));
	if (ImagCodePtr == NULL)
		mexErrMsgTxt("\n ImagCodePtr Array Not Allocated!--exiting \n");

	mxScrambleCodePtr = mxCreateDoubleMatrix(1,CHIPS_PER_FRAME,mxCOMPLEX);
	if (mxScrambleCodePtr == NULL)
		mexErrMsgTxt("\n mxScrambleCodePtr Matrix Not Allocated!--exiting \n");
	RealScrambleCodePtr = mxGetPr(mxScrambleCodePtr);
	ImagScrambleCodePtr = mxGetPi(mxScrambleCodePtr);

	mxBetaCPtr = mxCreateDoubleMatrix(1,1,mxREAL);
	if (mxBetaCPtr == NULL)
		mexErrMsgTxt("\n mxBetaCPtr Matrix Not Allocated!--exiting \n");
	BetaCPtr = mxGetPr(mxBetaCPtr);
	*BetaCPtr = (double) 1.0;

	mxBetaDPtr = mxCreateDoubleMatrix(1,1,mxREAL);
	if (mxBetaDPtr == NULL)
		mexErrMsgTxt("\n mxBetaDPtr Matrix Not Allocated!--exiting \n");
	BetaDPtr = mxGetPr(mxBetaDPtr);
	*BetaDPtr = (double) 1.0;

	for (i = 0;i < NumInterferers; i++)
	{
		//Generate Control Frame
		TempControlBitsPtr=ControlBitsPtr;
		for (j = 0; j < CONTROL_BITS; j++)
		{
			fee = ((double) rand() / (double) RAND_MAX) - 0.5;
			if (fee <= 0.0)	*TempControlBitsPtr++ = (double) -1.0;
			else *TempControlBitsPtr++ = (double) 1.0;
		}
/******************************************************************
		//DEBUG CODE
		TempControlBitsPtr=ControlBitsPtr;
		fp = fopen("C:/MATLABR11/work/LGICtest/testdata.txt","r");
		for (j = 0; j < CONTROL_BITS; j++)
		{
			fscanf(fp, "%e\n", &sam);
			*TempControlBitsPtr++ = (double) sam;
		} 
		fscanf(fp, "%e\n", &sam);
		ScrambleCodeNumber = (unsigned long) sam;
		fclose (fp);


//DEBUG CODE
/***************************************************************************/

		//Generate Scramble code for each interferer


		ScrambleCodeNumber =unsigned long (((double) rand() / (double) RAND_MAX)* (pow( (double) 2.0, (double) 24.0) - (double) 1));

		scramble_long(ScrambleCodeNumber, RealCodePtr, ImagCodePtr);

		//Store Scramble Code in an mxArray
		TmpRealCodePtr = RealCodePtr;
		TmpImagCodePtr = ImagCodePtr;
		TmpRealScrambleCodePtr = RealScrambleCodePtr;
		TmpImagScrambleCodePtr = ImagScrambleCodePtr;

		for (j = 0; j < CHIPS_PER_FRAME; j++)
		{
			*TmpRealScrambleCodePtr++ = (double) *TmpRealCodePtr++;
			*TmpImagScrambleCodePtr++ = (double) *TmpImagCodePtr++;
		}

		*mxRightHandSide = mxNumChanPtr;
		*(mxRightHandSide+1) = mxScrambleCodePtr;
		*(mxRightHandSide+2) = (mxArray *) mxChanCodePtr;
		*(mxRightHandSide+3) = mxControlBitsPtr;
		*(mxRightHandSide+4) = mxBetaCPtr;
		*(mxRightHandSide+5) = mxBetaDPtr;
		N_rhs = 6;
		N_lhs = 2;

		ErrCode = mexCallMATLAB(N_lhs,mxLeftHandSide,N_rhs,mxRightHandSide,"gen_uplink_wcdma_frame");
		if(ErrCode != 0)
			mexErrMsgTxt("\nError Calling gen_uplink_wcdma_frame from MAI_Processr\n");

		//Prepare Parameters for caling GenWCDMAUplinkSignal
//		InphaseUplinkChipsPtr = mxGetPr(mxUplinkFramePtr);
//		QuadUplinkChipsPtr = mxGetPi(mxUplinkFramePtr);
		InphaseUplinkChipsPtr = mxGetPr(mxLeftHandSide[0]);
		QuadUplinkChipsPtr = mxGetPi(mxLeftHandSide[0]);

		TempRealNextPtr = RealNextSignalPtr + i*SamplesPerInterferer;
		TempImagNextPtr = ImagNextSignalPtr + i*SamplesPerInterferer;

		GenWCDMAUplinkSignal(InphaseUplinkChipsPtr,QuadUplinkChipsPtr,
							InphaseUplinkChipsPtr,QuadUplinkChipsPtr,
							InphaseUplinkChipsPtr,QuadUplinkChipsPtr,
							CHIPS_PER_FRAME,PulseShapePtr,PulseLength,SamplesPerChip,
							TempRealNextPtr,TempImagNextPtr,SamplesPerInterferer);

		//Generate Offset

		OffsetDouble = (double) rand() / (double) RAND_MAX -  0.5;
		OffsetDouble *= 2 * MaxOffset * SamplesPerChip;
		if ((OffsetDouble - floor(OffsetDouble)) < 0.5)
			Offset = (int) OffsetDouble;
		else
			Offset = (int) OffsetDouble + 1;

		//************************
		//DEBUG CODE
		//Offset = 10;
		//DEBUG CODE
		//************************

		TempRealUplinkSignalPtr = RealUplinkSignalPtr;
		TempImagUplinkSignalPtr = ImagUplinkSignalPtr;

		if (Offset > (int) MaxDelay)
		{
			TempRealLastPtr = RealLastSignalPtr + (i+1)*NLastSignal - Offset;
			TempImagLastPtr = ImagLastSignalPtr + (i+1)*NLastSignal - Offset;
			TempRealCurrentPtr = RealCurrentSignalPtr + i*NCurrentSignal;
			TempImagCurrentPtr = ImagCurrentSignalPtr + i*NCurrentSignal;

			for (j = 0; j < (unsigned) Offset; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealLastPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagLastPtr++;
			}
			
			Stop = NCurrentSignal + MaxDelay - Offset;
			for (j = 0; j < Stop; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealCurrentPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagCurrentPtr++;
			}
		}
		else if (Offset >= 1)
		{
			TempRealLastPtr = RealLastSignalPtr + (i+1)*NLastSignal - Offset;
			TempImagLastPtr = ImagLastSignalPtr + (i+1)*NLastSignal - Offset;
			TempRealCurrentPtr = RealCurrentSignalPtr + i*NCurrentSignal;
			TempImagCurrentPtr = ImagCurrentSignalPtr + i*NCurrentSignal;
			TempRealNextPtr = RealNextSignalPtr + i*NNextSignal;
			TempImagNextPtr = ImagNextSignalPtr + i*NNextSignal;

			for (j = 0; j < (unsigned) Offset; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealLastPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagLastPtr++;
			}
			
			for (j = 0; j < NCurrentSignal; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealCurrentPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagCurrentPtr++;
			}

			Stop = MaxDelay - Offset;
			for (j = 0; j < Stop; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealNextPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagNextPtr++;
			}
		}
		else if (Offset < 0)
		{
			TempRealCurrentPtr = RealCurrentSignalPtr + i*NCurrentSignal - Offset;
			TempImagCurrentPtr = ImagCurrentSignalPtr + i*NCurrentSignal - Offset;
			TempRealNextPtr = RealNextSignalPtr + i*NNextSignal;
			TempImagNextPtr = ImagNextSignalPtr + i*NNextSignal;

			Stop = NCurrentSignal + Offset;
			for (j = 0; j < Stop; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealCurrentPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagCurrentPtr++;
			}

			Stop = MaxDelay - Offset;
			for (j = 0; j < Stop; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealNextPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagNextPtr++;
			}
		}

		else  //Offset = 0
		{
			TempRealCurrentPtr = RealCurrentSignalPtr + i*NCurrentSignal;
			TempImagCurrentPtr = ImagCurrentSignalPtr + i*NCurrentSignal;
			TempRealNextPtr = RealNextSignalPtr + i*NNextSignal;
			TempImagNextPtr = RealNextSignalPtr + i*NNextSignal;

			for (j = 0; j < NCurrentSignal; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealCurrentPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagCurrentPtr++;
			}

			for (j = 0; j < MaxDelay; j++)
			{
				*TempRealUplinkSignalPtr++ = *TempRealNextPtr++;
				*TempImagUplinkSignalPtr++ = *TempImagNextPtr++;
			}
		}

		if (M_PATH_FLAG)
		{
			//ApplyMultipathChannel

			//First Initialize output vector

			TempRealMpathPtr = RealMpathSignalPtr;
			TempImagMpathPtr = ImagMpathSignalPtr;
			N_UplinkSignal = NCurrentSignal + MaxDelay;
			for (j = 0; j < N_UplinkSignal; j++)
			{
				*TempRealMpathPtr++ = (double) 0.0;
				*TempImagMpathPtr++ = (double) 0.0;
			}

			MAI_Mpath_Channel(RealUplinkSignalPtr,ImagUplinkSignalPtr,N_UplinkSignal,
							  DelayPtr,NumDelays,MPathAmpPtr,fDoppler,SamplesPerChip,
							  RealMpathSignalPtr,ImagMpathSignalPtr);
		}
		else	//Do Not apply Multipath 
		{
			TempRealMpathPtr = RealMpathSignalPtr;
			TempImagMpathPtr = ImagMpathSignalPtr;
			TempRealUplinkSignalPtr = RealUplinkSignalPtr;
			TempImagUplinkSignalPtr = ImagUplinkSignalPtr;

			N_UplinkSignal = NCurrentSignal + MaxDelay;
			for (j = 0; j < N_UplinkSignal; j++)
			{
				*TempRealMpathPtr++ = *TempRealUplinkSignalPtr++;
				*TempImagMpathPtr++ = *TempImagUplinkSignalPtr++;
			}
		}
			
		TempRealMAIPtr = RealMAISignalPtr;
		TempImagMAIPtr = ImagMAISignalPtr;
		TempRealMpathPtr = RealMpathSignalPtr;
		TempImagMpathPtr = ImagMpathSignalPtr;

		for (j = 0; j < N_UplinkSignal; j++)
		{
			*TempRealMAIPtr++ += *TempRealMpathPtr++;
			*TempImagMAIPtr++ += *TempImagMpathPtr++;
		}

	}
	mxFree(RealUplinkSignalPtr);
	mxFree(ImagUplinkSignalPtr);
	mxFree(RealMpathSignalPtr);
	mxFree(ImagMpathSignalPtr);
	mxDestroyArray(mxNumChanPtr);
	mxDestroyArray(mxControlBitsPtr);
	mxFree(RealCodePtr);
	mxFree(ImagCodePtr);
	mxDestroyArray(mxScrambleCodePtr);
	mxDestroyArray(mxBetaCPtr);
	mxDestroyArray(mxBetaDPtr);

}


			


?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
国产一区啦啦啦在线观看| 日韩美女视频一区| 日韩成人午夜电影| 在线综合亚洲欧美在线视频| 亚洲第一成年网| 欧美一级在线观看| 久久97超碰国产精品超碰| 精品成人a区在线观看| 国产一区二区三区久久悠悠色av | 亚洲人成精品久久久久久 | 国产精品麻豆视频| 91浏览器在线视频| 亚洲成人7777| 久久新电视剧免费观看| 成人一级片在线观看| 国产精品久久久久aaaa樱花| 色妹子一区二区| 爽好多水快深点欧美视频| 日韩欧美中文字幕公布| 国产成人精品亚洲日本在线桃色| 日韩一区中文字幕| 欧美日韩精品一区二区三区四区 | 欧美va亚洲va国产综合| 国产传媒日韩欧美成人| 最新日韩av在线| 欧美日韩精品一区二区三区| 精品夜夜嗨av一区二区三区| 国产女人aaa级久久久级| av成人老司机| 日本成人超碰在线观看| 国产精品欧美综合在线| 欧美日韩国产在线播放网站| 国产精品自在欧美一区| 亚洲伊人伊色伊影伊综合网| 精品国产人成亚洲区| 色哟哟欧美精品| 国产在线看一区| 一区二区三区在线不卡| 久久嫩草精品久久久精品一| 欧美性猛交xxxx黑人交| 成人午夜视频网站| 奇米777欧美一区二区| 国产精品家庭影院| 精品国产免费人成电影在线观看四季| 9l国产精品久久久久麻豆| 毛片不卡一区二区| 亚洲人精品午夜| 国产性做久久久久久| 制服.丝袜.亚洲.中文.综合| av一区二区不卡| 国产在线播放一区二区三区| 亚洲自拍另类综合| 欧美高清在线一区二区| 欧美一区二区免费| 在线免费观看日本一区| 国产精品一区二区黑丝| 蜜臀久久99精品久久久久久9 | 亚洲一区在线免费观看| 日本一区二区在线不卡| 日韩欧美123| 91精品欧美久久久久久动漫| 欧美最新大片在线看| 9色porny自拍视频一区二区| 国产剧情一区二区| 麻豆高清免费国产一区| 午夜在线电影亚洲一区| 一区二区在线看| 亚洲色图都市小说| 国产精品动漫网站| 中文字幕欧美日本乱码一线二线 | 欧美天堂一区二区三区| 91啪亚洲精品| 成人午夜电影小说| 国内精品在线播放| 激情五月婷婷综合| 另类综合日韩欧美亚洲| 麻豆久久久久久| 久久精品国产99国产精品| 蜜桃久久av一区| 激情小说欧美图片| 国产精品77777| 国产精品一品二品| 国产成人无遮挡在线视频| 国产精品一区二区91| 国产a久久麻豆| 成人h版在线观看| 波多野结衣中文字幕一区| av亚洲精华国产精华| 99国产精品久| 欧美性大战xxxxx久久久| 欧美午夜精品久久久| 制服丝袜在线91| wwwwww.欧美系列| 中文字幕乱码亚洲精品一区| 自拍偷拍国产亚洲| 亚洲地区一二三色| 奇米精品一区二区三区在线观看一| 久久国产尿小便嘘嘘| 国产激情视频一区二区三区欧美 | 美女mm1313爽爽久久久蜜臀| 久99久精品视频免费观看| 国产激情偷乱视频一区二区三区| 成人激情免费网站| 欧美吻胸吃奶大尺度电影| 欧美一级高清大全免费观看| 久久久不卡影院| 一区二区三区色| 久久精品国产网站| 成a人片亚洲日本久久| 欧美午夜宅男影院| 精品av综合导航| 亚洲视频免费观看| 美国欧美日韩国产在线播放| 成人免费毛片aaaaa**| 在线观看日韩毛片| 精品成人a区在线观看| 亚洲女与黑人做爰| 激情综合色综合久久综合| 91亚洲大成网污www| 91精品国产一区二区人妖| 国产午夜三级一区二区三| 亚洲精品视频自拍| 狠狠色狠狠色合久久伊人| 日本久久电影网| 精品国产91九色蝌蚪| 一区二区三区中文在线| 久久国产欧美日韩精品| 色婷婷国产精品久久包臀| 精品久久久久一区| 亚洲精品久久7777| 国产精品一区专区| 欧美日韩大陆一区二区| 国产女人aaa级久久久级| 午夜精品福利一区二区蜜股av| 国产成人久久精品77777最新版本| 欧美性色欧美a在线播放| 久久九九全国免费| 免费在线看成人av| 色久综合一二码| 欧美国产激情一区二区三区蜜月| 天天操天天色综合| 色8久久人人97超碰香蕉987| 国产亚洲欧美一区在线观看| 日韩电影在线一区二区三区| 91视频免费看| 欧美激情综合五月色丁香小说| 午夜一区二区三区视频| 成人精品视频网站| 精品国产一二三区| 五月激情综合网| 欧洲色大大久久| 亚洲人成7777| 成人av电影免费在线播放| 久久久久久久综合日本| 青青青爽久久午夜综合久久午夜| 欧美午夜一区二区| 亚洲丝袜制服诱惑| 波多野结衣在线aⅴ中文字幕不卡| 337p日本欧洲亚洲大胆色噜噜| 青娱乐精品在线视频| 欧美日韩aaaaa| 亚洲丰满少妇videoshd| 欧美最新大片在线看 | 国产三级精品三级在线专区| 久久99精品国产.久久久久久 | 色婷婷狠狠综合| 亚洲免费在线观看| 色综合一区二区| 中文字幕在线不卡一区 | 污片在线观看一区二区| 欧美亚洲禁片免费| 亚洲一区二区在线免费观看视频 | 欧美天堂一区二区三区| 亚洲欧美另类图片小说| 色综合中文字幕国产 | 国产福利一区二区三区视频在线| 久久综合色8888| 国产精品一区在线观看你懂的| 国产色产综合色产在线视频| 成人精品一区二区三区四区| 专区另类欧美日韩| 日本伦理一区二区| 偷偷要91色婷婷| 精品日韩欧美在线| 国产精品白丝jk白祙喷水网站 | 欧美一区二区三区思思人| 麻豆成人在线观看| 国产欧美日韩麻豆91| 99久久精品免费| 亚洲一区二区av在线| 日韩一区二区视频| 极品少妇xxxx精品少妇偷拍| 中文字幕免费不卡| 欧美色综合影院| 久久99精品网久久| 亚洲人成网站色在线观看| 777色狠狠一区二区三区| 国模冰冰炮一区二区| 国产精品久久一卡二卡| 欧美午夜宅男影院|