加载中…
个人资料
飞天小苗苗
飞天小苗苗
  • 博客等级:
  • 博客积分:0
  • 博客访问:75,673
  • 关注人气:306
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
留言
加载中…
访客
加载中…
博文

参考代码:

clear all;clc;close all;
%% 采样点
x=[0 :0.2:10];
y=erf(x);

%% 多项式曲线拟合
i=2;
while(1)
[p,S]=polyfit(x,y,i);
    if S.normr<=0.05
        break;
    else 
        i=i+1;
    end
end
p_x=polyval(p,x);

%% 绘图
plot(x,y,'ro',x,p_x,'b')
xlabel('x'),ylabel('y=sin(x)')
legend('y(x)','p(x)')
title([num2str(i),'th-degree polynomial'])
gtext(['p(x)=',poly2str(p,'
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 

题目:

给定数据
x=[0 :0.2:10];
y=erf(x);

要求最大残差模不超过0.05,求多项式拟合的结果。

阅读  ┆ 评论  ┆ 转载 ┆ 收藏 

参考答案:

第一题答案:
>> A=[9 0 -4 1 0 0 2];
>> dA=polyder(A);
>> dA_str=poly2str(dA,'n')

dA_str =

   54 n^5 - 16 n^3 3 n^2

 

第二题答案:
>> B=[-3 1 0 -13 2];
>> int_B=polyint(B,9);
>> int_B_sym=poly2sym(int_B,'s')
 
int_B_sym =
 
- (3*s^5)/5 s^4/4 - (13*s^2)/2 2*s 9

阅读  ┆ 评论  ┆ 转载 ┆ 收藏 


阅读  ┆ 评论  ┆ 转载 ┆ 收藏 

参考答案:

第一题答案:
>> p=[1 6 0 2 1 1];
>> p_2=polyval(p,2)

p_2 =

   139

第二题答案:

>> p=[1 6 0 2 1 1];
>> r=roots(p)

r =

  -6.0509 0.0000i
   0.3783 0.6219i
   0.3783 - 0.6219i
  -0.3529 0.4329i
  -0.3529 - 0.4329i

第三题答案:

>> r=[0.5 -1];
>> p=poly(r)

=

  Columns through 3

    1.0000   -1.5000   -1.5000

  Column 4

    1.0000
第四题答案:
>> p1=[2 3 1];
>> p2=[5 0 1];
>> prod=conv(p1,p2)
prod =
    10    15       
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 


阅读  ┆ 评论  ┆ 转载 ┆ 收藏 





参考代码:

clear all;clc;close all;

[X,Y,Z,V]=flow(5);
figure;
slice(X,Y,Z,V,[1 9.5],2,[-1 3]);
shading flat

[Xq,Yq,Zq]=meshgrid(0.1:0.5:10,-3:0.5:3,-3:0.5:3);
Vq=interp3(X,Y,Z,V,Xq,Yq,Zq);
figure;
slice(Xq,Yq,Zq,Vq,[5 9],2,[-2 2])
shading flat

阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
使用flow(5)生成X,Y,Z,V的数据,使用slice命令绘制X=1,9.5 Y=2 Z=-1,3 五处位置的切面,绘制粗糙切面图。使用interp3命令绘制精细切面图(使用meshgrid命令生成新的三变量,增量设定为0.5)。
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
参考代码:
clear all;clc;close all;

service 10:10:30;
years 1950:10:1990;

wage [150.697 199.592 203.625;
179.323 195.072 250.287;
203.212 200.092 322.767;
226.505 237.706 426.730
249.633 257.281 598.243];
w_16_1988 interp2(service,years,wage,16,1988)
 
运行结果:
w_16_1988 =

  250.0226
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
根据工作起始年和工龄,对应相应的工资,表格如下

用二维插值法 找出1988年开始工作16年的工资情况



阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
  

新浪BLOG意见反馈留言板 电话:4000520066 提示音后按1键(按当地市话标准计费) 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 会员注册 | 产品答疑

新浪公司 版权所有