找回密码
 注册
关于网站域名变更的通知
查看: 284|回复: 3
打印 上一主题 下一主题

fminbnd循环求和

[复制链接]
  • TA的每日心情
    开心
    2022-1-21 15:21
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2023-2-6 11:23 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    代码一直在报错,不知道问题在哪里+ k1 d7 q' Q$ T2 R
    $ J5 {, v0 K, m
    clc;
    ) }& V5 [  F) e: w& r' E& jclear all;0 y1 Y! c" z! {1 G- f* P
    syms t x;" V, a* W0 ^) Q: W) f
    M=[2 4 8 16 32 64 128 256];
    0 A" P0 E- w6 a) @1 H6 ?/ s4 |. D. pv=10.5;0 X, E4 [! z+ J# p) |
    sigma=10^(-3);* }8 Q* F" g1 v
    z = sym('z');% Q0 @4 |7 t8 ^% r9 {
    x = sym('x');. Y7 O1 B( ?) m  M
    P0=10^(-3); %目标误码率
    + t, }% a* C' g* a/ V, V2 Lt=exp(-1/sigma^2);  k" Y0 M7 f9 n* m; z5 t$ k
    r=20; %平均SNR
    , e* {3 {, p* A6 M+ d/ `R=10.^(r/10);
    5 Z' }5 v/ h0 K9 ow=0;2 o2 J. V$ L- u
    for i=1:72 r: L8 c6 F) F: a" L2 |
    / ^8 Y5 g- q: ]
            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);: Y3 D/ \2 E' p8 q' k1 _
            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)# _0 d; q9 N& s  |
    end
    . j1 C6 ^5 Y" W# N7 W' P' [% b( `% }) @/ M9 i
    [z,feval] = fminbnd( w, 0, 1 )
    % N7 ^, F$ Z! ], H% W
  • TA的每日心情
    开心
    2022-1-29 15:05
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2023-2-6 13:11 | 只看该作者
    仅供参考,匿名函数f在定义之后,调用时f(z)里的z要是具体数值,否则无法计算匿名函数值
  • TA的每日心情
    开心
    2022-1-29 15:03
  • 签到天数: 2 天

    [LV.1]初来乍到

    3#
    发表于 2023-2-6 13:21 | 只看该作者
    解决如下,可以运行,不过结果看上去有些问题,去查查你的目标函数以及参数是不是都写对了/ Y$ N" P$ C2 V9 `$ a* Y

    2 x6 ~9 n5 @4 l/ G- k& I) a# [* b5 Hclc;% A1 d6 }: a  L. ~4 Q
    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];: S4 Y) D6 M4 ~
    v=10.5;
    ) i7 }0 ?2 H+ ]1 jsigma=10^(-3);
    , J4 ?  O# m: E/ g( v3 V- o% z = sym('z');  p' B# h, i; O# x
    % x = sym('x');, X2 Q4 D* L0 `. N% `: H, A/ s& T
    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! P% `7 p( Z, @- v' ~
    R=10.^(r/10);
      J3 E1 M$ D2 t) Qw=0;2 |; \" e4 t# A  {, a  G" F
    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);- p. {1 @& t& ]0 a
            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- H- `6 T. A* r7 H( ^; k& c" x+ M& L
                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);& H! t: J* S- h3 B5 w" ^. q/ J9 a
            end- d- g) u( g1 \# {) G
            
    ( q6 l0 V% y4 A3 Qend
    ) }* Z+ f1 l: ?5 S* D% {4 k/ B3 K( `7 X- n$ V
    [z,feval] = fminbnd( w, 0, 1 );
  • TA的每日心情
    开心
    2022-1-21 15:08
  • 签到天数: 1 天

    [LV.1]初来乍到

    4#
    发表于 2023-2-6 13:26 | 只看该作者
    检查一下目标函数及参数设置
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-28 00:10 , Processed in 0.125000 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表