# Plotting in MATLAB

1. Sep 13, 2011

### bear_25

Hello, i am new in using MATLAB, and i stuck with these two graphs (attached). Can anyone help me, please. Thank you.

#### Attached Files:

• ###### img001.jpg
File size:
25.4 KB
Views:
164
2. Sep 13, 2011

### Hootenanny

Staff Emeritus
Welcome to Physics Forums.

What exactly are you stuck with? Are you trying to plot them?

3. Sep 13, 2011

### bear_25

Thanks.

no, i am confused how to write this in Command Window in MATLAB, because for now i just start to use MATLAB.

4. Sep 13, 2011

### Hootenanny

Staff Emeritus
I'm afraid that I still don't understand. What exactly are you trying to do?

5. Sep 13, 2011

### bear_25

well, i need to plot this two graph in MATLAB, but i don t know how.

6. Sep 13, 2011

### Hootenanny

Staff Emeritus
So you do want to plot these, like I said earlier.

I'm not aware of any in-built function to plot shaded regions (you could plot the boundaries quite easily). However, the plotregion function from the File Exchange would likely do what you need: http://www.mathworks.com/matlabcentral/fileexchange/9261-plot-2d3d-region

7. Sep 13, 2011

### bear_25

Can you tell me how can i plot the boundaries? I really need help.

8. Sep 13, 2011

### Hootenanny

Staff Emeritus
I can't give you a full tutorial on MATLAB plotting. Instead, have a look at this: http://www.mathworks.co.uk/help/techdoc/ref/plot.html , try somethings and come back with specific questions.

9. Sep 13, 2011

### bear_25

Well, i have done this graph in MATLAB:
rA = 25;
kA = 2;
hAmax = 7;

rB = 40;
kB = 4;
hBmax = 27;

rC = 60;
kC = 5;
hCmax = 55;

R = 200;

i = 0;
%%
for dR = 0:0.001:4
i = i + 1;
hA(1,i) = rA*(1 - cos(asin((sqrt(2*R*dR + dR.^2) + rA*sin(acos((rA - kA)/rA)))/rA))) - kA;
hB(1,i) = rB*(1 - cos(asin((sqrt(2*R*dR + dR.^2) + rB*sin(acos((rB - kB)/rB)))/rB))) - kB;
hC(1,i) = rC*(1 - cos(asin((sqrt(2*R*dR + dR.^2) + rC*sin(acos((rC - kC)/rC)))/rC))) - kC;
end

%%
hA(hA<0) = 0;
hA(abs(hA)>hAmax) = hAmax;

hB(hB<0) = 0;
hB(abs(hB)>hBmax) = hBmax;

hC(hC<0) = 0;
hC(abs(hC)>hCmax) = hCmax;

%%
figure
plot(0:0.001:4,hA,'r')
hold on
plot(0:0.001:4,hB,'b')
plot(0:0.001:4,hC,'k')
xlabel('Bočna deformacija, dR/mm')
ylabel('Visina krojne naslage, h/mm')
grid on

Now i need to get those two graphs but i don t know how to set boundaries inside the graphs which are attached