当前位置: 首页 > news >正文

三叶草fz5722周口seo公司

三叶草fz5722,周口seo公司,专业的网站建设电话,室内设计软件哪个比较好1.简述 若向量记号为X,方程组就可以写成F(X)0的形式。 我们知道,对于一元函数的牛顿迭代法求根公式 类似的,对于多元函数求根公式 其中X是向量,是非线性方程组对应的雅可比矩阵。 具体求解的时候,我们可以先通过绘图命令绘制图形…

1.简述

      

若向量记号为X,方程组就可以写成F(X)=0的形式。

我们知道,对于一元函数的牛顿迭代法求根公式

类似的,对于多元函数求根公式

其中X是向量,是非线性方程组对应的雅可比矩阵。

具体求解的时候,我们可以先通过绘图命令绘制图形,看交点。然后将交点附近的值带入迭代

矩阵。最后求出小于误差的收链解。

2.代码

主程序:

function newton%牛顿迭代法解非线性方程组
    syms a x;
    beta=0.15;
    F=[1-beta+sqrt((1-beta)^2+4*a)-2*x;
       beta+a+2*x-(beta+x)^3];
    eps=10e-6;%精度
    num=1;%记步数
    tol=1;%给定误差初始值
    x0=[1;1];%参数赋初始值
    v=[a,x];
        while tol>eps
            Fx=subs(F,v,transpose(x0));
            dF=jacobian(F,v);
            c=subs(dF,v,transpose(x0));
            x=x0-inv(c)*Fx;
            tol=norm(x-x0);
            x0=x;
            num=num+1;
                if (num>10^8)
                    disp('迭代次数大于最大值,可能不收敛')
                    return
                end
        end
%Output
fprintf('解出 a = %g.\n',x0(1))
fprintf('解出 x = %g.\n',x0(2))
fprintf('迭代次数 n = %g次.\n',num)

子程序:

function   [y,n]=newton_fun(F,x0)
        if nargin==2
            eps=1.0e-6;
        end
        num = 0;
        tol = 1;
        v=findsym(F);
        while tol>eps
            Fx = subs(F,v,transpose(x0));
            dF=jacobian(F,v);
            c=subs(dF,v,transpose(x0));
            x=x0-inv(c)*Fx;
            tol=norm(x-x0);
            x0=x;
            num=num+1;
                if (num>10^8)
                    disp('迭代次数大于最大值,可能不收敛')
                    return
                end
        end
        y = x0;
        n = num;
end
 

子程序:

function solve
  syms x
    gamma=0.5;
    m=5;
    theta=12;
    a=1-gamma;
    F=sin((m+1)*x)-a*sin(m*x);
   
    for i=0:m+1
        w0=i*pi/(m+1);
        [x,n]=newton_fun(F,w0);
        w(i+1)=x;
%       beta(i+1)=theta*gama^2/(sqrt(1+a^2-2*a*cos(x))-gama*(1-theta));
        b(i+1)=cos((m+1)*x)-a*cos(m*x);%先算出b
        beta(i+1)=theta*gamma^2/(b(i+1)-gamma*(1-theta));%再算出beta
        num(i+1)=n;
    end
    w
    beta
    num

http://www.khdw.cn/news/25765.html

相关文章:

  • 睢宁县建设局网站视频app推广
  • 网站开发 浏览器最简单的营销方案
  • 做天猫网站价格seo搜索引擎优化兴盛优选
  • 广州公司的网页怎么做的珠海百度seo
  • 济南网站建设推荐企优互联不错淘宝运营培训
  • 做一视频网站多少钱新闻稿代写平台
  • 做的比较简约的网站好网站
  • 广东网站制作多少钱seo排名点击工具
  • 怎么做游戏网站的宣传图片网站优化排名的方法
  • 怎样做个做外贸的网站最近一两天的新闻有哪些
  • 万网网站发布提高工作效率的软件
  • 如何seo网站挣钱螺蛳粉的软文推广
  • 旅游b2b网站建设百度推广销售员好做吗
  • 汕尾网站建设百度免费seo
  • 央企做的好的网站seo排名优化怎么样
  • 文章网站的一级二级怎么做免费发广告的软件
  • 合肥做企业建网站那家好软文推广代写代发
  • 老域名怎么做新网站山东网站建设
  • 专门做情侣装的网站杭州seo联盟
  • 做个网站要多少钱 一般经典软文案例100例
  • 河北建网站做网店自己怎么去推广
  • 卖鞋的网站建设思路新东方雅思培训机构官网
  • 关于网站平台建设调研的函媒体发稿网
  • 网站域名申请步骤线上广告投放方式
  • 建设工程施工合同示范文本2021谷歌seo价格
  • 南京整站优化平台推广计划
  • 邯郸做网站的公司推广网站排名优化seo教程
  • 建网站找那家企业好域名注册后怎么使用
  • 想接网站自己做百度云引擎搜索
  • 手机网站建设服务器什么是网络营销平台