| 查看: 604 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
waitingywy新虫 (小有名气)
|
[求助]
MATLAB出图小程序找错误，求帮忙 已有1人参与
|
|
|
这段代码无法出图，不是很懂，求大神指教 bats_codes_out=cell(5,M); sucess_rates=[]; m=10; sucess_rate_sum=0; fail_rate_sum=0; for i=1:1:m n=10; % round((1+epsilong)*K/M); mu_n=randsample([1:length(mu)],n,true,mu); D=max(mu_n); d_min=ceil(D/K); bats_codes=[]; for j=1:n bats_codes_out=bats_encode(d_min,j,q,codes,mu_n,K,g,L,M); bats_codes_bec=becpacket(M,e,bats_codes_out,g); bats_codes_in=bats_codes_inner(q,M,bats_codes_bec,g); bats_codes_bec_bec=becpacket(M,e,bats_codes_in,g); bats_codes=[bats_codes,bats_codes_bec_bec]; end [fail_rate,de_codes,x]=bats_decodes(bats_codes,M,n,K,L,g); sucess_rate(i)=1-fail_rate X(i)=x; sucess_rate_sum=sucess_rate_sum+sucess_rate(i); fail_rate_sum=fail_rate_sum+fail_rate; end sucess_aver=sucess_rate_sum/m; fail_rate_aver=fail_rate_sum/m; suc_one=0; for i=1:m if sucess_rate(i)==1 suc_one=suc_one+1; end end % for i=1:K % if ~(isempty(de_codes{2,i})) % j=j+1; % end % end % s1=1:length(mu); % subplot(121); % plot(s1,mu); % s2=1:length(mu1); % subplot(122); % plot(s2,mu1); |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有266人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
3楼2020-05-09 22:40:29
waitingywy
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1223
- 散金: 168
- 红花: 3
- 帖子: 247
- 在线: 64.7小时
- 虫号: 2220691
- 注册: 2013-01-04
- 性别: GG
- 专业: 零件成形制造
|
clear all; clc; g=8; q=2^g; K=100; L=2; e=0.1; M=16; codes=randsrc(L,K,[0:q-1]); codes=gf(codes,g); mu=importdata('simDegreeM16K100p0.1.mat'); for i=1:size(mu,1) if i>K mu(K)=mu(K)+mu(i); mu(i)=0; end end sum(mu) bats_codes_out=cell(5,M); sucess_rates=[]; m=10; sucess_rate_sum=0; fail_rate_sum=0; for i=1:1:m n=10; % round((1+epsilong)*K/M); mu_n=randsample([1:length(mu)],n,true,mu); D=max(mu_n); d_min=ceil(D/K); bats_codes=[]; for j=1:n bats_codes_out=bats_encode(d_min,j,q,codes,mu_n,K,g,L,M); bats_codes_bec=becpacket(M,e,bats_codes_out,g); bats_codes_in=bats_codes_inner(q,M,bats_codes_bec,g); bats_codes_bec_bec=becpacket(M,e,bats_codes_in,g); bats_codes=[bats_codes,bats_codes_bec_bec]; end [fail_rate,de_codes,x]=bats_decodes(bats_codes,M,n,K,L,g); sucess_rate(i)=1-fail_rate X(i)=x; sucess_rate_sum=sucess_rate_sum+sucess_rate(i); fail_rate_sum=fail_rate_sum+fail_rate; end sucess_aver=sucess_rate_sum/m; fail_rate_aver=fail_rate_sum/m; suc_one=0; for i=1:m if sucess_rate(i)==1 suc_one=suc_one+1; end end % for i=1:K % if ~(isempty(de_codes{2,i})) % j=j+1; % end % end % s1=1:length(mu); % subplot(121); % plot(s1,mu); % s2=1:length(mu1); % subplot(122); % plot(s2,mu1); |
2楼2020-05-09 18:20:14