fminbnd循环求和

replace

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

) ~) c+ N3 g5 k3 [! iclc;
! [$ x7 l7 X* R' S5 w' qclear all;
  I- Z" O# G  S, v: O; e$ Dsyms t x;
M=[2 4 8 16 32 64 128 256];
1 K  N7 k, V7 F* d0 e- W7 v: Qv=10.5;
sigma=10^(-3);
3 `* t8 D8 d1 G  \2 Uz = sym('z');
x = sym('x');
; @$ |9 |9 f2 M1 Z% D) `/ KP0=10^(-3); %目标误码率
; r2 U, e8 _6 a( t8 ?7 Bt=exp(-1/sigma^2);
0 G! u" f3 |1 `9 K- s  ]r=20; %平均SNR
2 a- b) ]% \! I% fR=10.^(r/10);
* b# ^% W+ w/ c3 \: m9 Fw=0;
for i=1:7

        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);
        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))
        w=w+f(z)
end

[z,feval] = fminbnd( w, 0, 1 )

someone

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

I_believe

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

clc;
+ H% s& j, d0 B8 q: Tclear all;
, k# Y" S" L- D6 H9 nsyms t x;
: K0 _. I( m6 j: M8 M9 N; }M=[2 4 8 16 32 64 128 256];
v=10.5;
4 N# D+ y5 F) J, w8 N" e4 ^sigma=10^(-3);
% z = sym('z');
: ~8 F( Q0 e3 n" ~. I% x = sym('x');
5 C5 C- R  M: N0 @& m6 T$ S+ IP0=10^(-3); %目标误码率
( E8 N/ r7 q" Y7 x9 Tt=exp(-1/sigma^2);
* ^) Z/ \7 d7 V( K0 a! \8 S! Vr=20; %平均SNR
R=10.^(r/10);
w=0;
for i=1:7
$ q, P, j5 T$ X. X0 l7 X, ~+ S# Z
        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);
) c% [9 u: e2 z6 J        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))
        if i==1
            w=@(z) f(z);
        else
% z% U1 f; o1 d% K  D& _% C0 a            w=@(z) w(z)+f(z);
        end
        
4 x+ t0 y; w9 e. L8 X+ E: aend

7 d% ?, \9 z* ^( l+ B[z,feval] = fminbnd( w, 0, 1 );

R_myself

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