clear all;clc;
syms t;
f = sqrt(t);
fun1 = @(t) eval(f);
%%
a = 2.0;
b = 2.4;
n = 5;
for i=1:n
x(i)=a+(b-a)/(n-1)*(i-1);
y(i,1)= fun1(x(i));
end
%%
for j=1:n-1
for k=j:n-1
temp=y(k+1,j)-y(k,j);
y(k+1,j+1)=temp/(x(k+1)-x(k+1-j));
end
d_f(j)=y(j,j);
d_f(j+1)=y(j+1,j+1);
end
table = [x' y];
%%
c = d_f;
p=c(1);
q=1;
syms X
for i=2:n
q=q*(X-x(i-1));
p=p+c(i)*q;
end
p = simple(p)
p = vpa(p)
%%
for i=1:301
tt(i)=a+(b-a)/300*(i-1);
Nn(i) = eval(subs(p,'X','tt(i)'));
end
for i=1:301
h(i)=a+(b-a)/300*(i-1);
yy(i) = fun1(h(i));
end
%%
figure;
subplot(121);
plot(h,yy,'r')
title('f(x)')
subplot(122);
plot(tt,Nn,'b')
title('Newton interpolation polynomial')追问
你这回答的是什么？
...能用MATLAB帮忙画一下,下面函数所构成的曲面,以及曲线,万分感谢!
上面函数所构成的曲面，以及曲线可以用mesh（）函数和plot3（）函数来绘制。>> x=0:50;y=-25:25;>> [X,Y]=meshgrid(x,y);>> Z=25*(X-sqrt(631.25-Y.^2)).\/(16.7-sqrt(631.25-Y.^2));>> mesh(X,Y,Z) %3D曲面图 >> xlabel('x');ylabel('y');zlabel('z')>> x=...
matlab求多条曲线包络线。 见下图绿线该用什么语句画出? 万分感谢!
y=exp(x\/2); %你自己的函数！ 我是举例 y2=20*exp(x\/3);y3=30*exp(x\/5);plot(x,y,x,y2,x,y3)Y=[y;y2;y3];ymin=min(Y);hold on plot(x,ymin,'*')
matlab求二阶微分方程数值解,万分感谢。顺便能给下matlab解微分方程数 ...
matlab求二阶微分方程数值解,万分感谢。顺便能给下matlab解微分方程数值解 我来答 1个回答 #热议# 已婚女性就应该承担家里大部分家务吗?米子701 2016-10-06 知道答主 回答量:13 采纳率:0% 帮助的人:2万 我也去答题访问个人页 关注 展开全部 本回答由提问者推荐 已赞过 已踩过< 你对这...
谁帮我用MATLAB计算一下矩阵最大特征值和最大特征值所对应的特征向量...
>> [d,v]=eig([1 1\/7 1\/3 7 1 1\/5 3 5 1])d = -0.1327 -0.0663 - 0.1149i -0.0663 + 0.1149i -0.4094 -0.2047 + 0.3546i -0.2047 - 0.3546i -0.9026 0.9026 0.9026 v = 3.7089 0 0 0 -0.3545 + 1.5823i ...
用matlab求该函数的最小值,万分感谢!忘记怎么搜索了。。
求导=0，取x ,然后带入比较min(f(x))
你好,有个关于matlab simulink仿真的问题困扰我很久了,麻烦帮我看下...
1、建议用变步长求解器，例如ode45。如果用定步长求解器，最大步长不要设置为auto，设置稍微小一些的数，例如0.1。2、直接使用变步长求解器可能会因代数环的存在导致出错，这是因为PID控制器中的du\/dt模块是直接馈通的（Direct Feedthrough）。考虑使用 PID Controller (with Approximate Derivative) 模块...
如何使用Matlab,对一组数据进行FFT变换,得到频谱分析,万分感谢。
我的调试信号，你自己是电流电压数据的话，最开始通过load指令载入就是 Y = fft(S,N); %做FFT变换 Ayy = abs(Y); %取模 Ayy=Ayy\/(N\/2); %换算成实际的幅度 Ayy(1)=Ayy(1)\/2;F=([1:N]-1)*Fs\/N; %换算成实际的频率值，Fn=(n-1)*Fs\/N stem(F(1:N\/2),Ayy(1:...
matlab中这个方程式用fzero怎么求解,希望大神帮忙一下,万分感谢!急急...
题主给出的方程式，用fzero（）函数求解。可以这样做：1、自定义给出方程的函数 fun=@(x)100-(x*sqrt((x\/10)^2+1)*sqrt((x\/100)^2+1));2、设定x0初值，这里取 x0=1 3、然后用fzero函数求解 4、实现代码与结果 fun=@(x)100-(x*sqrt((x\/10)^2+1)*sqrt((x\/100)^2+1));...
谁能帮我注解一下这段MATLAB程序 万分感谢!!
I=filter(b,a,I);%%得到滤波器的系数b和a，并对数据I进行滤波 for i=1:COUNT%%COUNT应该是个数吧，还有STEP [maxx lag]=max(abs(I((i-1)*STEP+1:i*STEP)));%%在I中取STEP大小的窗口，并计算窗口内的最大值，maxx是最大值，lag是最大值所在的列 if maxx>=0.4 n=n+1;%%如果最...
如何用MATLAB将10个点拟合出一条曲线?求详细步骤、解题思路!万分...
如果是空间中的10个点，就在命令窗口中输入sftool，敲回车，打开空间曲线拟合工具箱；使用拟合工具箱是比较简单的办法，你也可以自己在命令框里使用matlab命令完成曲线拟合，polyfit二项式拟合，lsqcurvefit曲线拟合，regress一元回归，即使没有工具箱，matlab也可以通过各种命令完成任务，不要像楼上一样偷懒哦 ...
