fminbnd循环求和

replace

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

clc;
9 c$ d6 |* W/ W" e* k: I6 e# rclear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
9 i9 i. }2 ~8 y( Q0 l7 Mv=10.5;
sigma=10^(-3);
z = sym('z');
3 B2 R6 Y% [( ~& ^& fx = sym('x');
5 v8 n; P6 J. ?( _$ S6 e! Z7 p: j9 UP0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
r=20; %平均SNR
. l" z% X- P4 Y, ZR=10.^(r/10);
w=0;
4 G& j, u) b# N7 G$ d3 v0 ufor i=1:7
6 }4 \; o* M' ~4 b
* c6 S9 |4 Q4 c/ z0 h9 H* f9 D        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
! r2 x4 X: p& e% C8 E! M* N) S        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)
; ^8 o) \2 f, u& z1 f( U. J5 bend

, L% a* R# U+ R2 p1 w[z,feval] = fminbnd( w, 0, 1 )

someone

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

I_believe

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

clc;
clear all;
, f( e$ q4 ]2 ?, K4 xsyms t x;
M=[2 4 8 16 32 64 128 256];
) P, F8 Y7 ]+ l. I! a$ f/ h7 T1 nv=10.5;
% q& X3 p7 S. o, ^sigma=10^(-3);
9 X9 U5 g: k$ o; h% z = sym('z');
! h' I/ C8 g3 [) y( ^% x = sym('x');
" _8 ~6 H0 @; w3 o# x# H$ lP0=10^(-3); %目标误码率
& b, ~9 p" q* V, ht=exp(-1/sigma^2);
r=20; %平均SNR
R=10.^(r/10);
  w1 |; a3 B# ^) {/ Aw=0;
for i=1:7
7 I# R  `$ L- p: L. m3 L# x- L
# Q! Y+ O7 i  Z$ q        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);
6 L# }  u7 U2 J+ ]0 g+ b9 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))
2 i3 B8 D% S. q        if i==1
4 x# K$ p, D- Y            w=@(z) f(z);
        else
            w=@(z) w(z)+f(z);
) `$ d5 ]/ e3 f$ i3 {        end
        
! t( M: |( ~) d& b6 Rend
4 p  z) {7 ~  u3 j
) g: T1 {* e/ ~0 R[z,feval] = fminbnd( w, 0, 1 );

R_myself

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