Contour mapping not applied to all points.As you can see from the figure provided. The result of my code is the Contour mapping of the noise sources using pcolor and fcontour

As you can see from the figure provided. The result of my code is the Contour mapping of the noise sources using pcolor and fcontour command. But, the some of the points on the bottom and the right side are not covered with contours. Additinally, the two perpendicular line was also supposed to display contour color lines. But it does not. Please Help!!

Code:

m=201;

n=101;

x=0:10:1000;

X=repmat(x,m,1);

b=0:10:2000;

y=b';

Y=repmat(y,1,n);

%xy-axis coordinates

%Conveyors

x1=800;

y1=1600;

%First Truck

x2=700;

y2=200;

%Second Truck

x3=600;

y3=1400;

%Third Truck

x4=200;

y4=1400;

%First Grader

x5=600;

y5=1000;

%Second Grader

x6=400;

y6=1000;

%Third Grader

x7=200;

y7=800;

%First Dozer

x8=600;

y8=500;

%Second Dozer

x9=400;

y9=600;

%Third Dozer

x10=300;

y10=1500;

%First Crushing Plant

x11=900;

y11=1500;

%Second Crushing Plant

x12=700;

y12=1700;

%First Washing Plant

x13=900;

y13=200;

%Second Washing Plant

x14=200;

y14=1700;

%Distances away from Origin

%Conveyors

Rc1=sqrt(((0).^2+(y1-Y).^2));

Rc2=sqrt(((x1-X).^2)+((0).^2));

%Trucks

Rt1=sqrt(((x2-X).^2)+((y2-Y).^2));

Rt2=sqrt(((x3-X).^2)+((y3-Y).^2));

Rt3=sqrt(((x4-X).^2)+((y4-Y).^2));

%Grader

Rg1=sqrt(((x5-X).^2)+((y5-Y).^2));

Rg2=sqrt(((x6-X).^2)+((y6-Y).^2));

Rg3=sqrt(((x7-X).^2)+((y7-Y).^2));

%Dozers

Rd1=sqrt(((x8-X).^2)+((y8-Y).^2));

Rd2=sqrt(((x9-X).^2)+((y9-Y).^2));

Rd3=sqrt(((x10-X).^2)+((y10-Y).^2));

%Crushing Plant

Rcp1=sqrt(((x11-X).^2)+((y11-Y).^2));

Rcp2=sqrt(((x12-X).^2)+((y12-Y).^2));

%Washing Plant

Rwp1=sqrt(((x13-X).^2)+((y13-Y).^2));

Rwp2=sqrt(((x14-X).^2)+((y14-Y).^2));

%Sound Power Levels of Mining Conveyors with Standard Steel Idler Rolls

Lw_31_5= 0.0;

Lw_63=82;

Lw_125=93;

Lw_250=95;

Lw_500=105;

Lw_1000=106;

Lw_2000=103;

Lw_4000=96;

Lw_8000=82;

%Crushing Plant Sound Power Levels

LwC_31_5=76;

LwC_63=87;

LwC_125=95;

LwC_250=105;

LwC_500=110;

LwC_1000=111;

LwC_2000=109;

LwC_4000=103;

LwC_8000=91;

%WashingPlant Sound Power Levels

LwW_31_5=91;

LwW_63=95;

LwW_125=99;

LwW_250=107;

LwW_500=110;

LwW_1000=114;

LwW_2000=115;

LwW_4000=110;

LwW_8000=98;

% Horizontal Conveyor Sound Pressure Level

Lp_31_5_1=(Lw_31_5)-(10*log10(pi()*(Rc1)));

Lp_63_1=(Lw_63)-(10*log10(pi()*(Rc1)));

Lp_125_1=(Lw_125)-(10*log10(pi()*(Rc1)));

Lp_250_1=(Lw_250)-(10*log10(pi()*(Rc1)));

Lp_500_1=(Lw_500)-(10*log10(pi()*(Rc1)));

Lp_1000_1=(Lw_1000)-(10*log10(pi()*(Rc1)));

Lp_2000_1=(Lw_2000)-(10*log10(pi()*(Rc1)));

Lp_4000_1=(Lw_4000)-(10*log10(pi()*(Rc1)));

Lp_8000_1=(Lw_8000)-(10*log10(pi()*(Rc1)));

Log_Total_C1=(10.^(Lp_31_5_1/10))+10.^(Lp_63_1/10)...

+10.^(Lp_125_1/10)+10.^(Lp_250_1/10)+ 10.^(Lp_500_1/10)...

+10.^(Lp_1000_1/10)+10.^(Lp_2000_1/10)+10.^(Lp_4000_1/10)...

+10.^(Lp_8000_1/10);

Lp_Total_C1=10*log10(Log_Total_C1);

%Vertical Conveyor Sound Pressure Level

Lp_31_5_2=(Lw_31_5)-(10*log10(pi()*(Rc2)));

Lp_63_2=(Lw_63)-(10*log10(pi()*(Rc2)));

Lp_125_2=(Lw_125)-(10*log10(pi()*(Rc2)));

Lp_250_2=(Lw_250)-(10*log10(pi()*(Rc2)));

Lp_500_2=(Lw_500)-(10*log10(pi()*(Rc2)));

Lp_1000_2=(Lw_1000)-(10*log10(pi()*(Rc2)));

Lp_2000_2=(Lw_2000)-(10*log10(pi()*(Rc2)));

Lp_4000_2=(Lw_4000)-(10*log10(pi()*(Rc2)));

Lp_8000_2=(Lw_8000)-(10*log10(pi()*(Rc2)));

Log_Total_C2=(10.^(Lp_31_5_2/10))+10.^(Lp_63_2/10)...

+10.^(Lp_125_2/10)+10.^(Lp_250_2/10)+10.^(Lp_500_2/10)...

+10.^(Lp_1000_2/10)+10.^(Lp_2000_2/10)+10.^(Lp_4000_2/10)...

+10.^(Lp_8000_2/10);

Lp_Total_C2=10*log10(Log_Total_C2);

%First Truck 777C

Lp_63_T1=122.44-(10*log10(2*pi()*(Rt1.^2)));

Log_Total_T1=(10.^(Lp_63_T1));

Lp_Total_T1=10*log10(Log_Total_T1);

%Second Truck 777D

Lp_63_T2=122.44-(10*log10(2*pi()*(Rt2.^2)));

Log_Total_T2=(10.^(Lp_63_T2));

Lp_Total_T2=10*log10(Log_Total_T2);

%Third Truck Hitachi

Lp_63_T3=125.46-(10*log10(2*pi()*(Rt3.^2)));

Log_Total_T3=(10.^(Lp_63_T3));

Lp_Total_T3=10*log10(Log_Total_T3);

% First GraderCaterpillar 16H

Lp_63_G1=117.12-(10*log10(2*pi()*(Rg1.^2)));

Log_Total_G1=(10.^(Lp_63_G1));

Lp_Total_G1=10*log10(Log_Total_G1);

%Second Grader Caterpillar 16H

Lp_63_G2=117.12-(10*log10(2*pi()*(Rg2.^2)));

Log_Total_G2=(10.^(Lp_63_G2));

Lp_Total_G2=10*log10(Log_Total_G2);

%Third Grader Caterpillar 16H

Lp_63_G3=117.12-(10*log10(2*pi()*(Rg3.^2)));

Log_Total_G3=(10.^(Lp_63_G3));

Lp_Total_G3=10*log10(Log_Total_G3);

%First Dozer Cat D10T

Lp_63_D1=120.53-(10*log10(2*pi()*(Rd1.^2)));

Log_Total_D1=(10.^(Lp_63_D1));

Lp_Total_D1=10*log10(Log_Total_D1);

%Second Dozer Cat D10T

Lp_63_D2=120.53-(10*log10(2*pi()*(Rd2.^2)));

Log_Total_D2=(10.^(Lp_63_D2));

Lp_Total_D2=10*log10(Log_Total_D2);

%Third Dozer Cat D10T

Lp_63_D3=120.53-(10*log10(2*pi()*(Rd3.^2)));

Log_Total_D3=(10.^(Lp_63_D3));

Lp_Total_D3=10*log10(Log_Total_D3);

%Crushing Plant 1

Lp_31_5_Cp1=(LwC_31_5)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_63_Cp1=(LwC_63)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_125_Cp1=(LwC_125)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_250_Cp1=(LwC_250)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_500_Cp1=(LwC_500)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_1000_Cp1=(LwC_1000)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_2000_Cp1=(LwC_2000)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_4000_Cp1=(LwC_4000)-(10*log10(2*pi()*(Rcp1.^2)));

Lp_8000_Cp1=(LwC_8000)-(10*log10(2*pi()*(Rcp1.^2)));

Log_Total_Cp1=(10.^(Lp_31_5_Cp1/10))+10.^(Lp_63_Cp1/10)...

+10.^(Lp_125_Cp1/10)+10.^(Lp_250_Cp1/10)+10.^(Lp_500_Cp1/10)...

+10.^(Lp_1000_Cp1/10)+10.^(Lp_2000_Cp1/10)+10.^(Lp_4000_Cp1/10)...

+10.^(Lp_8000_Cp1/10);

Lp_Total_Cp1=10*log10(Log_Total_Cp1);

%Crushing Plant 2

Lp_31_5_Cp2=(LwC_31_5)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_63_Cp2=(LwC_63)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_125_Cp2=(LwC_125)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_250_Cp2=(LwC_250)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_500_Cp2=(LwC_500)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_1000_Cp2=(LwC_1000)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_2000_Cp2=(LwC_2000)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_4000_Cp2=(LwC_4000)-(10*log10(2*pi()*(Rcp2.^2)));

Lp_8000_Cp2=(LwC_8000)-(10*log10(2*pi()*(Rcp2.^2)));

Log_Total_Cp2=(10.^(Lp_31_5_Cp2/10))+(10.^(Lp_63_Cp2/10))...

+10.^(Lp_125_Cp2/10)+10.^(Lp_250_Cp2/10)+10.^(Lp_500_Cp2/10)...

+10.^(Lp_1000_Cp2/10)+10.^(Lp_2000_Cp2/10)+10.^(Lp_4000_Cp2/10)...

+10.^(Lp_8000_Cp2/10);

Lp_Total_Cp2=10*log10(Log_Total_Cp2);

%Washing Plant 1

Lp_31_5_Wp1=(LwW_31_5)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_63_Wp1=(LwW_63)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_125_Wp1=(LwW_125)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_250_Wp1=(LwW_250)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_500_Wp1=(LwW_500)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_1000_Wp1=(LwW_1000)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_2000_Wp1=(LwW_2000)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_4000_Wp1=(LwW_4000)-(10*log10(2*pi()*(Rwp1.^2)));

Lp_8000_Wp1=(LwW_8000)-(10*log10(2*pi()*(Rwp1.^2)));

Log_Total_Wp1=(10.^(Lp_31_5_Wp1/10))+10.^(Lp_63_Wp1/10)...

+10.^(Lp_125_Wp1/10)+10.^(Lp_250_Wp1/10)+10.^(Lp_500_Wp1/10)...

+10.^(Lp_1000_Wp1/10)+10.^(Lp_2000_Wp1/10)+10.^(Lp_4000_Wp1/10)...

+10.^(Lp_8000_Wp1/10);

Lp_Total_Wp1=10*log10(Log_Total_Wp1);

%Washinng Plant 2

Lp_31_5_Wp2=(LwW_31_5)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_63_Wp2=(LwW_63)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_125_Wp2=(LwW_125)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_250_Wp2=(LwW_250)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_500_Wp2=(LwW_500)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_1000_Wp2=(LwW_1000)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_2000_Wp2=(LwW_2000)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_4000_Wp2=(LwW_4000)-(10*log10(2*pi()*(Rwp2.^2)));

Lp_8000_Wp2=(LwW_8000)-(10*log10(2*pi()*(Rwp2.^2)));

Log_Total_Wp2=(10.^(Lp_31_5_Wp2/10))+10.^(Lp_63_Wp2/10)...

+10.^(Lp_125_Wp2/10)+10.^(Lp_250_Wp2/10)+10.^(Lp_500_Wp2/10)...

+10.^(Lp_1000_Wp2/10)+10.^(Lp_2000_Wp2/10)+10.^(Lp_4000_Wp2/10)...

+10.^(Lp_8000_Wp2/10);

Lp_Total_Wp2=10*log10(Log_Total_Wp2);

%Overall Total Noise Combined

Log_Total_Overall=(10.^(Lp_Total_C1/10))+(10.^(Lp_Total_C2/10))...

+(10.^(Lp_Total_T1/10))+(10.^(Lp_Total_T2/10))+(10.^(Lp_Total_T3/10))...

+(10.^(Lp_Total_G1/10))+(10.^(Lp_Total_G2/10))+(10.^(Lp_Total_G3/10))...

+(10.^(Lp_Total_D1/10))+(10.^(Lp_Total_D2/10))+(10.^(Lp_Total_D3/10))...

+(10.^(Lp_Total_Cp1/10))+(10.^(Lp_Total_Cp2/10))+(10.^(Lp_Total_Wp1/10))...

+(10.^(Lp_Total_Wp2/10));

Lp_Total_Overall=10*log10(Log_Total_Overall);

Z=Lp_Total_Overall;

%Colour Contour Plot for the Mining Site

pcolor(X,Y,Z);

shading flat

shading interp

contourf(X,Y,Z,10);

colormap(jet)

colorbar

set(gca,'ydir','reverse')

title('Mining Site Sound Pressure Levels')

xlabel('x-Coordinate (m)')

ylabel('y-Coordiante (m)')

Answer 1

[Danish Teo]

This is the Solution