| 查看: 761 | 回复: 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); |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有71人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有258人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
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
3楼2020-05-09 22:40:29
