Scribd
Upload
40 views

Numerical Analysis Lab 6: Syed Muhammad Hussain Rizvi 241736 Bscs-8a

This document summarizes the results of a numerical analysis lab involving curve fitting and interpolation. [1] The document fits a linear curve to the given time and value data using polyfit, and determines the slope and intercept. [2] It uses polyval to estimate y-values at two time points, and interp1 to interpolate those same points, showing the differences. [3] The difference between polyval and interp1 is explained - polyval uses the fitted linear curve while interp1 interpolates based on the original data points.

Uploaded by

Hussain Rizvi
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
40 views

Numerical Analysis Lab 6: Syed Muhammad Hussain Rizvi 241736 Bscs-8a

This document summarizes the results of a numerical analysis lab involving curve fitting and interpolation. [1] The document fits a linear curve to the given time and value data using polyfit, and determines the slope and intercept. [2] It uses polyval to estimate y-values at two time points, and interp1 to interpolate those same points, showing the differences. [3] The difference between polyval and interp1 is explained - polyval uses the fitted linear curve while interp1 interpolates based on the original data points.

Uploaded by

Hussain Rizvi
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

Numerical Analysis Lab 6

Syed Muhammad Hussain Rizvi


241736
Bscs-8a

Task code : (Missing Code and Updated previous qs code)

%Given Data
sprintf('<-----Part 1----->')

t=[0:8];
y=[40.12 66.78 80.17 86.71 80.77 66.78 44.41 10.51 -32.60];
% get polynomial coefficients (pc) of bets fit line
pc = polyfit(t,y,1); % 1 means linear fit
plot(t,y,'ro');
hold on;
plot(t,polyval(pc,t),'b-');
txt = sprintf('Best fit line y=%.2fx + %.2f -> Part-1',pc(1),pc(2));

%1. What is the slope and intercept? Hint: relate to polyfit coefficients
intercept=sprintf('y-intercept : %.2fx',pc(2))
gradient=sprintf('slope : %.2fx',pc(1))

%2. What is the value of y when t=4.5 and t=8.5 Hint: polyval
sprintf('<-----Part 2----->')

plot(4.5,polyval(pc,4.5),'bo');
plot(8.5,polyval(pc,8.5),'bo');
sprintf('Value a t= 4.5 : %.2f',polyval(pc,4.5))
sprintf('Value at t=8.5 : %.2f',polyval(pc,8.5))

%3. Use interp1 for part 2 - why the different values


%x2 = [(4.5) (8.5)];
sprintf('<-----Part 3----->')

interp1(t,y,4.5);
plot(4.5,interp1(t,y,4.5),'go');
plot(8.5,interp1(t,y,8.5),'go');
sprintf('Interpolation at t=8.5 : %.2f',interp1(t,y,4.5))
sprintf('Interpolation at t=8.5 : %.2f',interp1(t,y,8.5))
plot(x2,y2linear,'go');
legend('Data points -> Part-1',txt,'Value at t=4.5 -> Part-2','Value at t=8.5 -> Part-
2','Interpolation at t=4.5 -> Part-3');
1

Part 1 console output:

Part 2 console output:

Part 3 console output:


2

Graph output:

Q) Why the different values in part 2 and part 3

Ans : This is because in part 2 we use polyval which estimates the value according to the line
draw by polyfit where and the interpl function estimates the values at given points using the
imaginary curve drawn with points as shown above

You might also like