| 查看: 639 | 回复: 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); |
回复此楼» 猜你喜欢
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有161人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有1人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17808.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
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
