function [x,f,n]=shuzhiniudun(x0,escep,nita,N)
%x0给定初始值,escep根的容许误差,nita是函数的容许误差,N最大迭代次数
syms X
fun=X^3+X^2-3*(X+1);
g=diff(fun,X);
gx=double(subs(g,X,x0));
f=double(subs(fun,X,x0));
x1=x0-f/gx;
n=1;
while(abs(f)>nita||abs(x1-x0)>escep)
x0=x1;
gx=double(subs(g,X,x0));
f=double(subs(fun,X,x0));
x1=x0-f/gx;
n=n+1;
if(gx==0||n>N)
display('算法失败');
break;
end
end
x=x1;

f=double(subs(fun,X,x));

end

1个回答

10 个月之前 回复

C语言牛顿迭代法，正常运行输出错误结果，求助

python新人，牛顿迭代中循环判断">"和"<"有什么分别

#include<iostream> #include<conio.h> using namespace std; //x(n + 1) = (xn+a/xn)/2 void NewtonSqrt(double a,double &x){ a = x - (x + a / x) / 2; if (a<= 0.00001) return; x = (x + a / x) / 2; NewtonSqrt(a, x); } void main(){ double a = 3; double x = a / 2; cout << x <<endl; _getch(); }

![图片说明](https://img-ask.csdn.net/upload/201811/28/1543408731_531502.jpg) ![图片说明](https://img-ask.csdn.net/upload/201811/28/1543408623_611407.png) ![图片说明](https://img-ask.csdn.net/upload/201811/28/1543408639_341222.jpg) 伪代码如最后一张图黑色框内，请问应该怎么编matlab程序呢？ 或者有谁知道论文 Optimizing the Gaussian kernel function with the formulated kernel target alignment criterion for two-class pattern classification 的代码吗？

matlab实现牛顿法逻辑回归

# include <stdio.h> # include <math.h> # include <stdlib.h> # define e 2.714 int main () { double a,b; b=1; a=b-(log(4)*4); while (fabs(b-a)>1*e-6) { b=a; a=b-log(b+3)/(1/(b+3)); } printf("a=%ld");

#include <iostream> #include <cmath> #include <iomanip> using namespace std; int main() { double a, x0, x1, t; int n = 1; cin >> a; x0 = a / 2; x1= 1.0/ 2 * (x0 + a / x0); while (x1 - x0 >= pow(10, -5)) { t = x1; x1 = 1.0/ 2 * (x0 + a / x0); x0 = x1; n++; } cout << fixed<<setprecision(6)<<x1; return 0; }
y=ln(x+3),x0=1,求方程的根，用牛顿法。我的答案总乱码啊，为什么
#include<stdio.h> #include<math.h> int main( ) { int M=1000,k; double x0,x1=1,eps=0.000001; k=0; for(;;) { x0=x1; k++; printf("%.6f\n",x1); x1=x0-log(x0+3)*(x0+3); if(k>=M||fabs(x1-x0)<=eps) break; } if(fabs(x1-x0)<=eps) printf("%.6f\n",x1); }

![图片说明](https://img-ask.csdn.net/upload/201905/03/1556864456_931596.png) 优化问题 类似如图的问题 需要用 梯度下降法和牛顿法求解 ``` function [ outcome ] = f( w,c,A,b ) n=size(A,1); m=size(A,2); outcome=0; for i=1:m outcome=outcome+log(1+exp(-b(i)*(w'*A(:,i)+c))); end outcome=(outcome+0.01*(w'*w+c*c))/m; end ``` 这样的函数如何去对他进行求导或者求HESSIAN 是否函数不能这样定义？

<p><span style="font-size: small;">最近在面试中，碰到了一道关于数值分析的题，由于本人水平有限，没解答出来。还请各位大哥大姐帮忙瞧瞧~~</span></p> <div><span style="font-family: Verdana;"><span><span style="font-size: small;"><span>已知单调连续函数y=f(x)</span><span>的如下数据：</span></span></span></span></div> <div><span style="font-family: Verdana;"><span style="font-size: small;">x(i)     -0.11   0.00  1.50 1.80</span></span></div> <div><span style="font-size: small;">f(x(i))    -1.23    -0.10  1.17   1.58</span></div> <div><span style="font-size: small;"> </span></div> <div><span style="font-family: Verdana;"><span> <p class="MsoNormal"><span style="font-size: small;"><span>求若用插值法计算</span><span style="font-size: 12pt;" lang="EN-US"><span>，<img src="http://courseware.lzu.edu.cn/upload/shuzhifenxi/stjda2.files/image022.gif" border="0" alt="" width="13" height="15"></span></span><span>约为多少时</span><span style="font-size: 12pt;" lang="EN-US"><span><img src="http://courseware.lzu.edu.cn/upload/shuzhifenxi/stjda2.files/image024.gif" border="0" alt="" width="57" height="21"></span></span><span>（小数点后保留</span><span style="font-size: 12pt;" lang="EN-US">5</span><span>位）。</span></span></p> <p class="MsoNormal"><span style="font-size: x-small;"> </span></p> <p class="MsoNormal"><span style="font-size: x-small;"> </span></p> <p class="MsoNormal"><span style="font-size: x-small;"> </span></p> </span></span></div><br /><strong>问题补充：</strong><br />注：用牛顿插值法

win10系统安装教程（U盘PE+UEFI安装）

《奇巧淫技》系列-python！！每天早上八点自动发送天气预报邮件到QQ邮箱

Nginx 软件层面加强Nginx性能优化的面试问答和解决方案
Nginx 软件层面加强Nginx性能优化的面试问答和解决方案 去年我去爱卡汽车面试PHP，一轮和二轮面的都不错，在三轮面到Nginx的时候很多问题当时不知道怎么回答，确实没有深入学习过，花了一段时间的学习，终于能解答Nginx高性能优化的问题了，10月24号为了获得程序员勋章，发布了半个优化笔记，浏览到了1000+，受到这个鼓舞，我抽时间在仔细整理下关于Nginx性能优化的问题，我们从软件说起。...
【管理系统课程设计】美少女手把手教你后台管理
【文章后台管理系统】URL设计与建模分析+项目源码+运行界面 栏目管理、文章列表、用户管理、角色管理、权限管理模块（文章最后附有源码） 1. 这是一个什么系统? 1.1 学习后台管理系统的原因 随着时代的变迁，现如今各大云服务平台横空出世，市面上有许多如学生信息系统、图书阅读系统、停车场管理系统等的管理系统，而本人家里就有人在用烟草销售系统，直接在网上完成挑选、购买与提交收货点，方便又快捷。 试想，若没有烟草销售系统，本人家人想要购买烟草，还要独自前往药...
11月19日科技资讯|华为明日发布鸿蒙整体战略；京东宣告全面向技术转型；Kotlin 1.3.60 发布
「极客头条」—— 技术人员的新闻圈！ CSDN 的读者朋友们早上好哇，「极客头条」来啦，快来看今天都有哪些值得我们技术人关注的重要新闻吧。扫描上方二维码进入 CSDN App 可以收听御姐萌妹 Style 的人工版音频哟。 一分钟速览新闻点！ 6G 专家组成员：速率是 5G 的 10 至 100 倍，预计 2030 年商用 雷军：很多人多次劝我放弃WPS，能坚持下来并不是纯粹的商业决定 ...
C语言魔塔游戏

Python爬虫爬取淘宝，京东商品信息

Java工作4年来应聘要16K最后没要,细节如下。。。

2020年，冯唐49岁：我给20、30岁IT职场年轻人的建议

##1、骇客帝国(1999) 概念：在线/离线，递归，循环，矩阵等 剧情简介： 不久的将来，网络黑客尼奥对这个看似正常的现实世界产生了怀疑。 他结识了黑客崔妮蒂，并见到了黑客组织的首领墨菲斯。 墨菲斯告诉他，现实世界其实是由一个名叫“母体”的计算机人工智能系统控制，人们就像他们饲养的动物，没有自由和思想，而尼奥就是能够拯救人类的救世主。 可是，救赎之路从来都不会一帆风顺，到底哪里才是真实的世界？

CPU对每个程序员来说，是个既熟悉又陌生的东西？ 如果你只知道CPU是中央处理器的话，那可能对你并没有什么用，那么作为程序员的我们，必须要搞懂的就是CPU这家伙是如何运行的，尤其要搞懂它里面的寄存器是怎么一回事，因为这将让你从底层明白程序的运行机制。 随我一起，来好好认识下CPU这货吧 把CPU掰开来看 对于CPU来说，我们首先就要搞明白它是怎么回事，也就是它的内部构造，当然，CPU那么牛的一个东
@程序员，如何花式构建线程？

2020年1月17日，国家统计局发布了2019年国民经济报告，报告中指出我国人口突破14亿。 猪哥的朋友圈被14亿人口刷屏，但是很多人并没有看到我国复杂的人口问题：老龄化、男女比例失衡、生育率下降、人口红利下降等。 今天我们就来分析一下我们国家的人口数据吧！ 一、背景 1.人口突破14亿 2020年1月17日，国家统计局发布了 2019年国民经济报告 ，报告中指出：年末中国大陆总人口（包括31个

Linux自学篇——linux命令英文全称及解释
man: Manual 意思是手册，可以用这个命令查询其他命令的用法。 pwd：Print working directory 显示当前目录 su：Swith user 切换用户，切换到root用户 cd：Change directory 切换目录 ls：List files 列出目录下的文件 ps：Process Status 进程状态 mk
Python实战：抓肺炎疫情实时数据，画2019-nCoV疫情地图

NO.1 　有20瓶药丸，其中19瓶装有1克/粒的药丸，余下一瓶装有1.1克/粒的药丸。给你一台称重精准的天平，怎么找出比较重的那瓶药丸？天平只能用一次。 解法 有时候，严格的限制条件有可能反倒是解题的线索。在这个问题中，限制条件是天平只能用一次。 因为天平只能用一次，我们也得以知道一个有趣的事实：一次必须同时称很多药丸，其实更准确地说，是必须从19瓶拿出药丸进行称重。否则，如果跳过两瓶或更多瓶药

2020年的1月，我辞掉了我的第一份工作

Java坑人面试题系列: 包装类（中级难度）
Java Magazine上面有一个专门坑人的面试题系列: https://blogs.oracle.com/javamagazine/quiz-2。 这些问题的设计宗旨，主要是测试面试者对Java语言的了解程度，而不是为了用弯弯绕绕的手段把面试者搞蒙。 如果你看过往期的问题，就会发现每一个都不简单。 这些试题模拟了认证考试中的一些难题。 而 “中级(intermediate)” 和 “高级(ad
Spring Boot 2.x基础教程：使用JdbcTemplate访问MySQL数据库

By 超神经场景描述：昨天 2 月 3 日，是大部分城市号召远程工作的第一天，全国有接近 2 亿人在家开始远程办公，钉钉上也有超过 1000 万家企业活跃起来。关键词：十一出行 人脸...