应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助高手帮忙用matlab接下面这个算式，具体见附件
61/1返回列表
查看: 795  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

likuihei

金虫 (正式写手)

[求助] 求助高手帮忙用matlab接下面这个算式，具体见附件

本人在科研过程中找到了这个公式，并做试验验证，公式的计算规则如下，不知道要怎么用matlab计算，请高手帮忙算一下：

AN=((0.02105*xN^0.9208)/(1- xN)^0.9208-aN)*0.0002778                   1式
BN=(1.515*0.2842* xN /(1+0.2842* xN)-bN)*0.0002778                    2式
M=11.8554*A+2.0937*B                                           3式
总式为0.0010652*(3600*0.8/300-3168*cN/298)-M=0                    4式                                 

①N=1时，a1 =0.007，b1=0.048，x1 =0.0903
将a1、b1、x1代入解出4式中的cN值为c1；

②N=2时a2=0.007+ A1
        b2=0.048+ B1
        x2= c1
代入取出值c2

③ N=3时a3=0.007+ A1+ A2
        b3=0.048+ B1+ B2
        x3= c2
代入取出值c3

④ N=4时a4=0.007+ A1+ A2+ A3
          b4=0.048+ B1+ B2+ B3
         x4= c3
代入取出值c4

⑤ N=5时a5=0.007+ A1+ A2+ A3+ A4
          b5=0.048+ B1+ B2+ B3+ B4
         x5= c4
代入取出值c5
以此类推。
请高手给出m文件并将计算的结果绘制成曲线，cN-N
或者给出表格形式的cN，N
若不理解请看附件[ Last edited by likuihei on 2011-10-3 at 22:05 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
每天重复同样的事情，需要很大的耐性
1楼2011-10-03 22:03:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

likuihei(金币+20): 2011-10-04 18:41:28
N = 10;
n = 1:N;
a = zeros(1,N);        b = zeros(1,N);
c = zeros(1,N);        A = zeros(1,N);
B = zeros(1,N);        x = zeros(1,N);

a(1) = 0.007;        b(1) = 0.048;        x(1) =0.0903;

for i=1:N
    A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778;
    B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778;
    M = 11.8554*A(i)+2.0937*B(i);
    c(i) = (3600*0.8/300-M/0.0010652)*298/3168;
    a(i+1) = a(i)+A(i);
    b(i+1) = b(i)+B(i);
    x(i+1) = c(i);
end
plot(n,c,'bo-')
赞一下 回复此楼
2楼2011-10-04 10:17:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likuihei

金虫 (正式写手)
引用回帖:
2楼: Originally posted by change0618 at 2011-10-04 10:17:45:
N = 10;
n = 1:N;
a = zeros(1,N);        b = zeros(1,N);
c = zeros(1,N);        A = zeros(1,N);
B = zeros(1,N);        x = zeros(1,N);

a(1) = 0.007;        b(1) = 0.048;        x(1) =0.0903;

for i=1:N
    A(i) = ((0.02105* ...

高手对不住，我已开始把计算要求搞错了，下面是正确的计算要求，请您再帮忙修改一下，谢谢！再加十个金币吧，可以吗
赞一下 回复此楼
每天重复同样的事情，需要很大的耐性
3楼2011-10-04 18:41:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
N = 50;
a = zeros(1,N);        b = zeros(1,N);
c = zeros(1,N);        A = zeros(1,N);
B = zeros(1,N);        x = zeros(1,N);
M = zeros(1,N);        S = zeros(1,N);
a(1) = 0.007;        b(1) = 0.048;        x(1) =0.0903;

for i=1:N
    A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778;
    B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778;
    M(i) = 11.8554*A(i)+2.0937*B(i);
    S(i) = sum(M(1:i));
    c(i) = (3600*0.8/300-S(i)/0.0010652)*298/3168;
    a(i+1) = a(i)+A(i);
    b(i+1) = b(i)+B(i);
    x(i+1) = c(i);
end
plot(1:N,c,'ro-')
赞一下 回复此楼
4楼2011-10-04 19:34:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likuihei

金虫 (正式写手)
引用回帖:
4楼: Originally posted by change0618 at 2011-10-04 19:34:00:
N = 50;
a = zeros(1,N);        b = zeros(1,N);
c = zeros(1,N);        A = zeros(1,N);
B = zeros(1,N);        x = zeros(1,N);
M = zeros(1,N);        S = zeros(1,N);
a(1) = 0.007;        b(1) =  ...

您好，我要怎么把金币给你呢，刚才20个给过了，再追加追加不上了，你知道怎么弄吗？不知道的话，下次有问题请教你再给你可以吗？
赞一下 回复此楼
每天重复同样的事情，需要很大的耐性
5楼2011-10-04 20:29:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)
6楼2011-10-04 21:17:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 likuihei 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭