fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里

clc;
) }& V5 [  F) e: w& r' E& jclear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
0 A" P0 E- w6 a) @1 H6 ?/ s4 |. D. pv=10.5;
sigma=10^(-3);
z = sym('z');
x = sym('x');
P0=10^(-3); %目标误码率
+ t, }% a* C' g* a/ V, V2 Lt=exp(-1/sigma^2);
r=20; %平均SNR
, e* {3 {, p* A6 M+ d/ `R=10.^(r/10);
5 Z' }5 v/ h0 K9 ow=0;
for i=1:7

        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
" b( E5 r1 K' e' ]        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
        f=@(z) integral(@(x) 0.5.*eRFc(0.5.*x.*v.*z.*t.^(z.^2.*x.^2).*sin(pi./M(i)).*(x./sigma.^2).*(t.^(-x.^2./2))),h(i),h(i+1))
9 e6 U% n! Z. U* D2 K        w=w+f(z)
end
. j1 C6 ^5 Y" W# N7 W' P
[z,feval] = fminbnd( w, 0, 1 )
% N7 ^, F$ Z! ], H% W

someone

仅供参考，匿名函数f在定义之后，调用时f(z)里的z要是具体数值，否则无法计算匿名函数值

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了

2 x6 ~9 n5 @4 l/ G- k& I) a# [* b5 Hclc;
clear all;
/ M3 {$ n! O; j* Vsyms t x;
8 G. k2 z1 b3 g6 K6 Y" I5 ^M=[2 4 8 16 32 64 128 256];
v=10.5;
) i7 }0 ?2 H+ ]1 jsigma=10^(-3);
, J4 ?  O# m: E/ g( v3 V- o% z = sym('z');
% x = sym('x');
P0=10^(-3); %目标误码率
8 Z+ E$ K2 m2 K. T6 \4 F% vt=exp(-1/sigma^2);
' Q+ {  L* U8 z$ l# D, \2 u; Hr=20; %平均SNR
R=10.^(r/10);
  J3 E1 M$ D2 t) Qw=0;
for i=1:7
- @& k- C$ n  m' Q& O' ^
0 M, n5 @2 l5 G! c2 ]" s# q9 T        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
) A- u' E0 Q+ D: W% B        f=@(z) integral(@(x) 0.5.*erfc(0.5.*x.*v.*z.*t.^(z.^2.*x.^2).*sin(pi./M(i)).*(x./sigma.^2).*(t.^(-x.^2./2))),h(i),h(i+1))
% F0 ]! G/ |, \        if i==1
            w=@(z) f(z);
  i# J4 }- O% d& Y, d# B( W; l7 _# U  G        else
8 I, f- {) W5 j. u            w=@(z) w(z)+f(z);
        end
        
( q6 l0 V% y4 A3 Qend
) }* Z+ f1 l: ?5 S* D% {4 k/ B
[z,feval] = fminbnd( w, 0, 1 );

R_myself

检查一下目标函数及参数设置