基于博弈论的频谱分配的matlab程序_第1页
基于博弈论的频谱分配的matlab程序_第2页
基于博弈论的频谱分配的matlab程序_第3页
基于博弈论的频谱分配的matlab程序_第4页
基于博弈论的频谱分配的matlab程序_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、是认知无线电场景下，基于博弈论的频谱分配的matlab程序-Cognitive radio scenario, based on game theory matlab spectrum allocation processclc clear all B=15; BER=10(-4); K=1.5/(logm(0.2/BER); SNR=8,10; k=log2(1+K.*SNR); r1=10; r2=12; x=0; y=1; z=1; % b=0,0; % c=x+y*(sum(b); % p=r.*k.*b-b.*c; a(1)=0.001; for i=1:300 a(2)=0.001

2、; for j=1:300 eigenvalue(1)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)+10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); eigenvalue(2)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)-10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)

3、*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); if (eigenvalue(1)<1)&&(eigenvalue(1)>-1)&&(eigenvalue(2)<1)&&(eigenvalue(2)>-1) a(2)=a(2)+0.002; end end if a(2)>0 a(2)=a(2)-0.002; else a(2)=0; end f2(i)=a(2); f1(i)=a(1); a(1)=a(1)+0.002; end figure(1) plot(f1,f2,'-r

4、9;) axis(0,0.25,0,0.25) hold on B=15; BER=10(-4); K=1.5/(logm(0.2/BER); SNR=10,10; k=log2(1+K.*SNR); r1=10; r2=12; x=0; y=1; z=1; % b=0,0; % c=x+y*(sum(b); % p=r.*k.*b-b.*c; a(1)=0.001; for i=1:300 a(2)=0.001; for j=1:300 eigenvalue(1)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)+10/3*(a(1)*(

5、k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); eigenvalue(2)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)-10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); if (eigenvalue(1)<1)&&(eigenvalue(1)>-1)&a

6、mp;&(eigenvalue(2)<1)&&(eigenvalue(2)>-1) a(2)=a(2)+0.002; end end if a(2)>0 a(2)=a(2)-0.002; else a(2)=0; end f2(i)=a(2); f1(i)=a(1); a(1)=a(1)+0.002; end figure(1) plot(f1,f2,'-.*b') axis(0,0.25,0,0.25) B=15; BER=10(-4); K=1.5/(logm(0.2/BER); SNR=10,8; k=log2(1+K.*SNR)

7、; r1=10; r2=12; x=0; y=1; z=1; % b=0,0; % c=x+y*(sum(b); % p=r.*k.*b-b.*c; a(1)=0.001; for i=1:300 a(2)=0.001; for j=1:300 eigenvalue(1)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)+10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); eigenvalue(2)=1-

8、10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)-10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); if (eigenvalue(1)<1)&&(eigenvalue(1)>-1)&&(eigenvalue(2)<1)&&(eigenvalue(2)>-1) a(2)=a(2)+0.002; end end if a(2)>0 a

9、(2)=a(2)-0.002; else a(2)=0; end f2(i)=a(2); f1(i)=a(1); a(1)=a(1)+0.002; end figure(1) plot(f1,f2,'-g') axis(0,0.25,0,0.25) B=15; BER=10(-4); K=1.5/(logm(0.2/BER); SNR=7,7; k=log2(1+K.*SNR); r1=10; r2=12; x=0; y=1; z=1; % b=0,0; % c=x+y*(sum(b); % p=r.*k.*b-b.*c; a(1)=0.001; for i=1:300 a(2

10、)=0.001; for j=1:300 eigenvalue(1)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)+10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); eigenvalue(2)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)-10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(

11、2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); if (eigenvalue(1)<1)&&(eigenvalue(1)>-1)&&(eigenvalue(2)<1)&&(eigenvalue(2)>-1) a(2)=a(2)+0.002; end end if a(2)>0 a(2)=a(2)-0.002; else a(2)=0; end f2(i)=a(2); f1(i)=a(1); a(1)=a(1)+0.002; end figure(1) plot(f1,f2,

12、9;-.y') axis(0,0.25,0,0.25) B=15; BER=10(-4); K=1.5/(logm(0.2/BER); SNR=8,8; k=log2(1+K.*SNR); r1=10; r2=12; x=0; y=1; z=1; % b=0,0; % c=x+y*(sum(b); % p=r.*k.*b-b.*c; a(1)=0.001; for i=1:300 a(2)=0.001; for j=1:300 eigenvalue(1)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)+10/3*(a(1)*(k(

13、2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); eigenvalue(2)=1-10/3*(2*a(1)*k(1)-a(1)*k(2)+2*a(2)*k(2)-a(2)*k(1)-10/3*(a(1)*(k(2)-2*k(1)*a(2)*(k(1)-2*k(2)+(2*a(2)*k(2)-a(2)*k(1)-2*a(1)*k(1)+a(1)*k(2)2)(0.5); if (eigenvalue(1)<1)&&(eigenvalue(1)>-1)&&(eigenvalue(2)<1)&&(eigenvalue(2)>-1) a(2)=a(2)+0.002; end end if a(2)>0 a(2)=a(2)-0.002; else a(2)=0; end f2(i)=a(2); f1(i)=a(1); a(1)=a(1)+0.002; end figure(1) plo

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论