fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里
. E# z) x! k# D5 i+ p
clc;
/ ?8 `5 R$ E# V% w$ O5 L; Aclear all;
8 A1 q# }6 h( Qsyms t x;
M=[2 4 8 16 32 64 128 256];
* B; F$ V# x$ `; V; c- vv=10.5;
2 o# s- e. r5 f: _sigma=10^(-3);
z = sym('z');
x = sym('x');
3 i: r, V  S, b2 _9 CP0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
# r# B0 a( Q" R" E# Ur=20; %平均SNR
R=10.^(r/10);
  X7 Z" _; e6 A2 i* q9 X" M! p2 Yw=0;
for i=1:7

% }0 Q3 b% `, @3 @4 v( ]# B        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))
) o8 Y& H& o% T        w=w+f(z)
end

& I* r" b! f- ?( [) x7 m' i2 F[z,feval] = fminbnd( w, 0, 1 )

someone

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
) o8 G$ ^/ J+ c4 x& i5 {1 Z4 b; r
clc;
clear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
v=10.5;
sigma=10^(-3);
% z = sym('z');
% x = sym('x');
4 }" T( X) K3 f+ H2 ]P0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
9 ~# ~- }0 K6 S" G& Er=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);
7 m. L( G% x# @+ q        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
6 |( N& i$ q; j4 g, d* @            w=@(z) f(z);
1 H, ~6 s7 o8 S* f, K        else
            w=@(z) w(z)+f(z);
        end
        
end
4 H9 u! r3 o+ c' o% X
[z,feval] = fminbnd( w, 0, 1 );

R_myself

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