来自 娱乐 2019-11-02 14:14 的文章
当前位置: 澳门皇冠金沙网站 > 娱乐 > 正文

昨11月国足球赢了

问题:杰士邦:今益求精,头头事到nBlackBerry:赢得优质nSamsung:恭喜中国足球,昨日看小编的n金立科学技术:只要有1%的或是,大家就有99%的信仰n雷蛇:追寻无可比拟,正是肩负14亿期望,寸步不让n安踏体育:不常交恶不比翻身n刚好食物:国“仁”自豪n……n

"增速是GDP的三倍! 高盛:现在华夏的投资时机在这里四大行业n3小时前  nn曹泽熙 小说总量 219篇n 关注n摘要:在华夏,值得投资的不单有BATJ,还会有新兴工业、新花费和医治常规等比相当多家底。n*本文来源华尔街视线(微信ID:wallstreetcn),小编曹泽熙。更加多优秀资源音信请登录wallstreetcn.com,或下载华尔街见闻APP。*nn高盛把最近几年来中夏族民共和国的退换归纳为多个定义:“新的中华”(New China),在此个概念中,新兴工业、网络、新费用和正规行当是值得关切的投资时机。nn高盛猜想,那四大行当的均匀拉长大概为60%,大致是中国GDP年增长速度的三倍。nn高盛将各样行业又细分了2-3个行业,并依附不一样行当的影响力举行了赋权:nnn新兴工业占比21%(此中高级成立业占8%,IT创设业占9%,清洁能源占3%);新开销占比32%(其Hong Kong中华电力有限集团动小车占2%,娱乐行当占14%,教育行业占16%);互连网占21%(个中电商、游戏占15%,网络金融占6%);健康行当占26%(当中医治健康服务占19%,医治安保卫证占7%)。n那四大行当有何特点呢?nn首先,增进将三回九转保持强硬。在过去三年中,那四大行当的年复合增进率为18%,高于GDP年复合增加率8%。我们感觉,在现在,那四大行业的年复合增加率将达19%。nnn第二,不仅局限在网络领域。中夏族民共和国的网络和手艺集团吸引了广大投资人目光。不过,我们只付与了互连网22%的权重。除了互连网,饱含了娱乐、旅游、教育等在内的新成本,以致健康行当和新颖工业行当等等,比重都比表现抢眼的互连网集团大。nnn第三,分裂于“新经济”。互连网和新兴工业的百分比占到了四分之二,而正规行业和新成本的增加率约为15%,那和“新经济”的定义差别。nnn第四,增进不要单纯是量的加强。那四大行业不知足于夺取越来越多的市镇占有率,同时还在商海上有所更加大的话语权,并在股票总市值链中处于越来越高等的未有。nn最终,有谈得来的周期和逻辑。那四大行当的抓好受到经济周期的震慑相当小,有投机的提升周期。nnn从宏观经济的角度来明白那第四次全国代表大会行当,最简易的章程实在把宏观经济数据拆分为第生龙活虎、二、第三行业业,以致投资和说话。之所以这四大行当能成为华夏的辽阳行业,高盛以为,主要基于以下多少个成分:nn首先,2013年来讲,第第三行当业的增长速度超过第第二行业业:nnn第二,中华夏族民共和国耗费者的花费行为发生肯定更动,网购正在吞吃更加多的线下购物市集分占的额数:nnn第三,越来越多的“软基本建设”(水务、教育、医疗、文娱体育)一败涂地,拉动中华新花费:nnn第四,在出口方面,中黄炎子孙民共和国的低附送值产品出口比重下落,高附赠值产品出口比重上涨:nnn此外,中夏族民共和国政府也会有生机勃勃多种安排推动那四大行业的前进。nn包蕴二〇一一年就写入“十九五”规划的七大新兴行业;贰零壹贰年国家总括局引入的五个总计门类:高科技(science and technology)成立业和高科技(science and technology)服务业,器重观测于研究开发投入;中夏族民共和国还在2014年建议了“中华夏族民共和国构建2025”安插,希望拉长创立业水平;今年,人民政坛还发表了推动人工智能演变的有关设计,希望在创制业、种植业、金融、教育以致行政事务等地点发布人工智能的职能。nn宏观数据表现能够、还大概有政策支持,那四大行业发表现状如何?nn高盛从以下多少个方面进行了拆解解析:nn第生龙活虎,智力资本。研发投入是直接体现智力资本投入的多少。二〇一六年,中华夏儿女民共和国在研究开发方面包车型地铁投入大概占到了GDP的2%,同繁多沸腾经济体万分。专利申请方面,二〇一五年,中华夏儿女民共和国申请专利数占到了满世界的29%,已经超(英文名:jīng chāo)越欧洲,并和美利坚合众国一定。nn第二,人力资本。中华夏族民共和国超级重视教育投入,2014年的教育投入大大略攻下GDP的3%。长时间的教化投入让中华有着丰硕的工程技术人力能源。nn第三,基础设备建设。在过去十年中,中华夏儿女民共和国高科学技术行个中央建设年复合增加率达13%。"

求解难题

Write a program to solve the following ordinary differential equation by

  • basic Euler method
  • improved Euler method
  • four-order Runge-Kutta method

澳门皇冠金沙网站 1=3 end{cases} xin[0, 1.5])

and calculate y(1.5) with stepsize=0.1, 0.1/2, 0.1/4, 0.1/8

Compare it with analytic solution (in figure)

澳门皇冠金沙网站 2 = dfrac{3}{1 x^3})

回答:

主程序

program main
    use ODE
    implicit none
    real :: x0=0.0,y0=3.0,a=0.0,b=1.5,h=0.1
    integer :: n,i,j
    character(len=512) :: filename


    print *,'The analytic solution of y(1.5)= ',analytic_solution(1.5)
    print *

    print *,'Basic Euler Method'
    print '(5x,a,f3.1,3x,a,f3.1,3x,a,f3.1,3x,a,f3.1)','Set initial  x0= ',x0,'y0= ',y0,'a= ',a,'b= ',b
    do j=0,3
        h=0.1/2**j
        print '(5x,a,f6.4)','Set h= ',h
        call init(x0,y0,a,b,h)
        n=size(x)-1
        call Basic_Euler
        print '(8x,a,f10.8)','y(1.5)= ',y(n)
        write(filename, *) j
        filename='BEM_'//Trim(AdjustL(filename))//'.txt'
        open(101,file=filename)
        write(101,'(f3.1,f10.8)') (x(i),y(i),i=0,n)
        close(101)
        deallocate(x,y)
    end do
    print *

    print *,'Improved Euler Method'
    print '(5x,a,f3.1,3x,a,f3.1,3x,a,f3.1,3x,a,f3.1)','Set initial  x0= ',x0,'y0= ',y0,'a= ',a,'b= ',b
    do j=0,3
        h=0.1/2**j
        print '(5x,a,f6.4)','Set h= ',h
        call init(x0,y0,a,b,h)
        n=size(x)-1
        call Improved_Euler
        print '(8x,a,f10.8)','y(1.5)= ',y(n)
        write(filename, *) j
        filename='IEM_'//Trim(AdjustL(filename))//'.txt'
        open(101,file=filename)
        write(101,'(f3.1,f10.8)') (x(i),y(i),i=0,n)
        close(101)
        deallocate(x,y)
    end do
    print *

    print *,'Four Order Runge-Kutta Method'
    print '(5x,a,f3.1,3x,a,f3.1,3x,a,f3.1,3x,a,f3.1)','Set initial  x0= ',x0,'y0= ',y0,'a= ',a,'b= ',b
    do j=0,3
        h=0.1/2**j
        print '(5x,a,f6.4)','Set h= ',h
        call init(x0,y0,a,b,h)
        n=size(x)-1
        call Four_Order_Runge_Kutta
        print '(8x,a,f10.8)','y(1.5)= ',y(n)
        write(filename, *) j
        filename='RKM_'//Trim(AdjustL(filename))//'.txt'
        open(101,file=filename)
        write(101,'(f3.1,f10.8)') (x(i),y(i),i=0,n)
        close(101)
        deallocate(x,y)
    end do
    print *

end program main

求解ODE方程的要害措施写在ODE模块中:

module ODE
    implicit none
    private
    public :: x,y,init,Basic_Euler,Improved_Euler,Four_Order_Runge_Kutta,analytic_solution

    real,allocatable :: x(:),y(:)
    real :: x0,y0,a,b,h
    integer :: n,i

contains

    function f(x,y)
        implicit none
        real :: f,x,y
        f=-x*x*y*y
    end function f

    function analytic_solution(x) result(f)
        implicit none
        real :: f,x
        f=3.0/(1 x*x*x)
    end function

    subroutine init(x0_,y0_,a_,b_,h_)
        implicit none
        real :: x0_,y0_,a_,b_,h_
        x0=x0_
        y0=y0_
        a=a_
        b=b_
        h=h_
        n=int((b-a)/h)
        allocate(x(0:n),y(0:n))
        x=(/ (a i*h,i=0,n) /)
        y=0
        y(0)=y0
    end subroutine init

    subroutine Basic_Euler()
        implicit none
        do i=1,n
            y(i)=y(i-1) h*f(x(i-1),y(i-1))
        end do
    end subroutine Basic_Euler

    subroutine Improved_Euler()
        implicit none
        real :: y_
        do i=1,n
            y_=y(i-1) h*f(x(i-1),y(i-1))
            y(i)=y(i-1) h/2*(f(x(i-1),y(i-1)) f(x(i),y_))
        end do
    end subroutine Improved_Euler

    subroutine Four_Order_Runge_Kutta()
        implicit none
        real :: k1,k2,k3,k4
        do i=1,n
            k1=f(x(i-1),y(i-1))
            k2=f(x(i-1) h/2,y(i-1) h/2*k1)
            k3=f(x(i-1) h/2,y(i-1) h/2*k2)
            k4=f(x(i-1) h,y(i-1) h*k3)
            y(i)=y(i-1) h/6*(k1 2*k2 2*k3 k4)
        end do
    end subroutine Four_Order_Runge_Kutta
end module ODE

其中init措施是用来初阶化:

    subroutine init(x0_,y0_,a_,b_,h_)
        implicit none
        real :: x0_,y0_,a_,b_,h_
        x0=x0_
        y0=y0_
        a=a_
        b=b_
        h=h_
        n=int((b-a)/h)
        allocate(x(0:n),y(0:n))
        x=(/ (a i*h,i=0,n) /)
        y=0
        y(0)=y0
    end subroutine init

analytic_solution是数值解:

    function analytic_solution(x) result(f)
        implicit none
        real :: f,x
        f=3.0/(1 x*x*x)
    end function

不是广告,小编只是以为这一次大家合营社的文案不错
澳门皇冠金沙网站 3

Basic Euler Method

流程图:

澳门皇冠金沙网站 4

原理:
![][3]
[3]: http://latex.codecogs.com/gif.latex?begin{cases} y(x_{n 1})=y(x_n) hf(x_n,y(x_n)) O(h^2) y(x_0)=y_0 end{cases}

代码:

    subroutine Basic_Euler()
        implicit none
        do i=1,n
            y(i)=y(i-1) h*f(x(i-1),y(i-1))
        end do
    end subroutine Basic_Euler

出口结果:

澳门皇冠金沙网站 5

Basic Euler 方法

图例:

澳门皇冠金沙网站 6

不一样h每步迭代结果

Improved Euler Method

流程图:

澳门皇冠金沙网站 7

原理:
![][4]
[4]: http://latex.codecogs.com/gif.latex?begin{cases} bar y_{n 1}=y_n hf(x_n,y_n) y(x_{n 1})=y(x_澳门皇冠金沙网站,n) dfrac{h}{2}[f(x_n,y(x_n)) f(x_{n 1},bar y_{n 1})] O(h^2) y(x_0)=y_0 end{cases}

代码:

    subroutine Improved_Euler()
        implicit none
        real :: y_
        do i=1,n
            y_=y(i-1) h*f(x(i-1),y(i-1))
            y(i)=y(i-1) h/2*(f(x(i-1),y(i-1)) f(x(i),y_))
        end do
    end subroutine Improved_Euler

出口结果:

澳门皇冠金沙网站 8

Improved Euler 方法

图例:

澳门皇冠金沙网站 9

差别h每步迭代结果

Four Order Runge-Kutta Method

流程图:

澳门皇冠金沙网站 10

原理:
![][5]
[5]: http://latex.codecogs.com/gif.latex?begin{cases} y_{n 1}=y_n dfrac{h}{6}(K_1 2K_2 2K_3 K_4) K_1=f(x_n,y_n) K_2=f(x_n dfrac{h}{2},y_n dfrac{h}{2}K_1) K_3=f(x_n dfrac{h}{2},y_n dfrac{h}{2}K_2) K_4=f(x_n h,y_n hK_3) end{cases}

代码:

    subroutine Four_Order_Runge_Kutta()
        implicit none
        real :: k1,k2,k3,k4
        do i=1,n
            k1=f(x(i-1),y(i-1))
            k2=f(x(i-1) h/2,y(i-1) h/2*k1)
            k3=f(x(i-1) h/2,y(i-1) h/2*k2)
            k4=f(x(i-1) h,y(i-1) h*k3)
            y(i)=y(i-1) h/6*(k1 2*k2 2*k3 k4)
        end do
    end subroutine Four_Order_Runge_Kutta

出口结果:

澳门皇冠金沙网站 11

4 Order Runge-Kutta 方法

图例:

澳门皇冠金沙网站 12

区别h每步迭代结果

三种情势结果比较

澳门皇冠金沙网站 13

h=0.1时

澳门皇冠金沙网站 14

h=0.1/8时

能够看看,当h非常的大时,二种艺术的间隔依然超大的,当h渐渐减小时,两种方法的结果已基本相近。

本文由澳门皇冠金沙网站发布于娱乐,转载请注明出处:昨11月国足球赢了

关键词: