fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里
4 e+ ?8 d* T; s
clc;
" l0 G' o7 X! L2 `9 Qclear all;
2 a5 p4 D5 m; Lsyms t x;
M=[2 4 8 16 32 64 128 256];
0 I, e$ L" `; r8 \v=10.5;
1 p: ?; ?4 G5 E0 ]2 ^5 y" F4 _1 Y% Y1 Bsigma=10^(-3);
z = sym('z');
x = sym('x');
P0=10^(-3); %目标误码率
. h5 t. V0 w+ P# z) z, \! l. Bt=exp(-1/sigma^2);
r=20; %平均SNR
5 C) x' ?5 H, Y1 x/ KR=10.^(r/10);
w=0;
for i=1:7
) K) `6 d& u$ A2 _! b
$ m$ Q& k6 u6 A        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);
3 ]7 P! j& n2 ^: Q+ d. ?) Q$ N        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 )
. W' I( a+ J! F# }& k/ {- T- d* z

someone

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
/ C* p; u: E" E9 F$ g
* d8 }. Y* ~0 J$ yclc;
clear all;
syms t x;
6 R- B. c8 F8 E8 U5 p; \M=[2 4 8 16 32 64 128 256];
, T( Z; ?; P) v6 n8 O4 ov=10.5;
5 m6 t; N  q; |. Q: ~' T$ g- K% Wsigma=10^(-3);
% z = sym('z');
% x = sym('x');
P0=10^(-3); %目标误码率
# v' G8 ]' V0 ^/ `; D& Y+ zt=exp(-1/sigma^2);
& g$ A0 w  h7 \r=20; %平均SNR
R=10.^(r/10);
w=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);
% [. C9 M& m- E  D" \        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);
  p- M  A& e- p        else
" T1 I" U4 l( q            w=@(z) w(z)+f(z);
/ ~% M6 l- h! y4 f8 S4 E  Y, g        end
        
0 J; _8 L: V( R# Q$ f( N; X6 k; oend
$ c% A% E% `# n/ b) E6 U
[z,feval] = fminbnd( w, 0, 1 );

R_myself

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