CURVE FIT USING MATLAB

-

 curve fit:

It is the process of constructing a curve or a mathematical function,and it has the relationship b/w orginal data set and predicted data set.

Basically,the curve fit depends upon the polynomial equ...

The polynomial equ are show below:

  1. linear equ/first degree polynomial:y=ax+b;
  2. quadratic equ/second degree polynomial:y=ax^2+bx+c;
  3. cubic equ/third degree polynomial:y=ax^3+bx^2+cx+d;
  4. the degree of polynomial can go 4,5 and so on;

coading for linear and quadratic equ:

clear all
close all
clc
%loading data
cp_data=load('data')
temp=cp_data(:,1)
cp=cp_data(:,2)

%curve fit
co_eff1=polyfit(temp,cp,3);
predicted_cp=polyval(co_eff1,temp)
co_eff2=polyfit(temp,cp,1);
predicted_cp2=polyval(co_eff2,temp)
%compare org fit,predicted fit
plot(temp,cp,LineWidth=3);
hold on
plot(temp,predicted_cp,LineWidth=3)
plot(temp,predicted_cp2,LineWidth=3)
xlabel('temperature')
ylabel('specific heat')
legend('original c_p','cubic polynomial','linear polynomial')

process:

First of all,we load the data into the matlab workspace....

In that data,the first column is temprature and 2nd column is specific heat;all that values assign to the variable called temp and cp...

Then,we will enter to the curve fit process;we wants to find co-effi for that data with the help of polyfit command,which is inbuild in matlab

And,we wants to find predicted value of y;so we will use the command polyval;then we will plot that org data and curve fit,then compare itself...

output:

it will also done by curve fitting tools

process:

  • Go to apps and click curve fitting tools,then simply select x data and y data,
  • It will shown projected line and goodness of fit like results also shown as it is..
  • And then click file and click generate code,the code automatically generates...

split the curve fit:

for clearing the error, we assign data to be splitted...

the code formation as it is..we just split the data...

clear all
close all 
clc
%loading a data
cp_data=load('data');
temp1=cp_data(:,1);
cp1=cp_data(:,2);
%spliting the curve fit
co_efficient1=polyfit(temp1(1:1600),cp1(1:1600),3);
co_efficient2=polyfit(temp1(1600:3200),cp1(1600:3200),3);
predicted_cp1=polyval(co_efficient1,temp1(1:1600))
predicted_cp2=polyval(co_efficient2,temp1(1600:3200))
%plotting org valuve and spliting curve fit value
plot(temp1,cp1,'linewidth',3,'color','g')
hold on
plot(temp1(1:1600),predicted_cp1,'LineWidth',3,'color','b')
plot(temp1(1600:3200),predicted_cp2,'LineWidth',3,'color','r')
xlabel('temp')
ylabel('specific heat')
legend('org data','split data 1','split data 2')

output:

HERE,WE COMPARE CURVE FIT AND SPLIT CURVE FIT GRAPH HAVE DIFFRENCES,THE LINE MUCH LIKE A SAME LINE....


Comments

Post a Comment

Popular posts from this blog

simple pendulum using matlab

-
Image
  ORDERLY DEFFRENTIAL EQUATION: In applied engineering and science many physical initial and boundary value problems can be described by ODE. The behavior of many physical problem,particularly those system undergoing time dependent changes can be formulated as ODE. The diff.equ involves one or more physical quantities such as displacement,velocity,temperature etc.. An nth order ODE may be represented as F(x,y,y’,y’’…)=0 SIMPLE PENDULUM: A simple pendulum is a mechanical device,the thin string or massless string is fixed in a supported beam and the other end is fixed to a bob.actually,its act like a newton’s 2 nd law,we give a force to the bob,the pendulum swig due to net force and mass.according to the net force and mass , the bob travels high or low. THE EQU OF A SIMPLE PENDULUM: b-damping co-efficient m-mass of the object g-gravity l-length og the string   First of all,we change the 2 nd order ODE to 1 ST order ODE  According to this formulla,create a function in ma
Read more

DRAG FORCE

-
Image
DRAG FORCE   DRAG FORCE: Drag force is a force,which resistance offered by the fluid like air or water.if we consider a car travelling a point A to B,we put a hand outside thorugh the window,then we feel the opposive of air.its known as drag force.   DRAG FORCE EQUATION: When the drag takes place in the presence of air, it is known as aerodynamic drag. When it takes place in water, it is known as hydrodynamic drag.   F D  = ½c_dρAv 2   Where, C_d-Drag co-efficient,which is dimensionless with value,depending upon the shape of object.(c_d=0.33-0.41 for mini cars and c_d=0.31-0.40 for sport cars etc...). ρ  -rho is the density of fluid. V – velocity of the object. A – cross sectional area of object( In automonile 'A' consider as frontal area). UNIT OF DRAG FORCE:NEWTON or N Drag force is propotional to the square of the velocity.if the velocity of speed increases two times,drag force increases four times. EXAMPLE: A car travelling on a road. A bicycle pedalling
Read more

2r robotic arm

-
Image
2R ROBOTIC ARM ABOUT: It is a type of mechanical arm,usually it’s perform a similar function like a human arm.in the 2r robotics arm,the two links and manipulator were connected by the joints and the joints allow to perform either rotational motion or linear displacement.the links of the manipulator can be consider to form a kinematic chain To create a animation using matlab: First of all,we should understand the concept of forward kinematics 2r robotics arm like how it’s works?,what are the material has been used?,how many linkages and joints has been used? Then,we should do analyse mathematically, ¨     Calculate length and position of link1.. ¨     Calculate length and position of link2.. For plotting purpose;we pick some co-ordinates,which refers to forward kinematics robotics arm. Starting orgin x0=0;y0=0 At link 1: x1=l1cos⍬;y1=l1sin⍬ At link 2:x2=x1+l2cos⍬;y2=y2+l2sin⍬ [x0 y0][x1 y1][x2 y2] write a code to create animation link 2r robotic arm: clear all close all clc %
Read more