Plotting a result with 3 different variables

20 Dec 2022
1 Answer
Answer Accepted
Updated 22 Dec 2022
3 Views (30 days)

0 votes

Hi everybody,
I have a equation that has three variables (Ca, Ka and W). I am trying to get the result for each Ca, Ka and W values and plot the results. If I assign 10 variables from 1 to 10 to each of them then I'd like to have 10*10*10=1000 results for my equation. I used linspace command but it gave me total of 10 results where each of my variables started from 1 and increased to 10 all together. Is there a way to create a code to assign those variables individually?
I am truly a newbie in Matlab so I apologize if I'm asking a really dumb question :)
clear all
close all
clc
Ca = linspace(1,10,10);
Ka = linspace(1,10,10);
W = linspace(1,10,10);
fun(Ca,Ka,W)
function Tr = fun(Ca,Ka,W,ma,m,k,Mu,W0,Zeta,Wa,Beta,ra,b,c,d,e,f,g)
ma=2;
m=40;
k=25;
Mu=ma/m;
W0=sqrt(k/m);
Zeta=Ca/(2*ma*W0);
Wa=sqrt(Ka/ma);
Beta= sqrt(Wa/W0);
r=W/W0;
Tr = sqrt(((2.*Zeta.*r).^2+(r.^2-Beta.^2).^2)./((2.*Zeta.*r).^2.*(r.^2-1+Mu.*r.^2).^2+(Mu.*r.^2.*Beta.^2-(r.^2-1).*(r.^2-Beta.^2)).^2));
plot (W,Tr)
end

1 Comment
Show -1 older comments Hide -1 older comments

Rik
Rik on 20 Dec 2022
You probably want ndgrid, although I don't fully understand what you want. Can you explain what results you want to see?

Sign in to comment.

Sign in to answer this question.

 Accepted Answer

Torsten
Torsten on 20 Dec 2022
Edited: Torsten on 20 Dec 2022
Ran in:

1 vote

Ran in:
To plot, you will have to keep one parameter constant and plot "result" over the two remaining parameters.
Ca = linspace(1,10,10);
Ka = linspace(1,10,10);
W = linspace(1,10,10);
[ca,ka,w]=ndgrid(Ca,Ka,W);
result = fun(ca,ka,w)
result =
result(:,:,1) = 1.6465 1.5682 1.4523 1.3668 1.3246 1.3097 1.3081 1.3121 1.3183 1.3251 1.5293 1.4962 1.4659 1.4413 1.4227 1.4094 1.4003 1.3942 1.3904 1.3881 1.4984 1.4821 1.4685 1.4573 1.4480 1.4405 1.4344 1.4296 1.4257 1.4227 1.4866 1.4771 1.4694 1.4630 1.4576 1.4530 1.4491 1.4457 1.4428 1.4404 1.4810 1.4747 1.4698 1.4657 1.4622 1.4591 1.4564 1.4541 1.4520 1.4501 1.4778 1.4735 1.4701 1.4672 1.4647 1.4625 1.4606 1.4589 1.4573 1.4559 1.4759 1.4727 1.4702 1.4681 1.4663 1.4646 1.4632 1.4619 1.4607 1.4596 1.4747 1.4722 1.4703 1.4687 1.4673 1.4660 1.4649 1.4639 1.4629 1.4620 1.4738 1.4719 1.4704 1.4691 1.4680 1.4670 1.4661 1.4652 1.4645 1.4638 1.4732 1.4716 1.4704 1.4694 1.4685 1.4677 1.4669 1.4662 1.4656 1.4650 result(:,:,2) = 0.1859 0.1866 0.1872 0.1878 0.1883 0.1888 0.1892 0.1897 0.1902 0.1906 0.1837 0.1840 0.1842 0.1844 0.1845 0.1846 0.1847 0.1848 0.1848 0.1849 0.1814 0.1815 0.1815 0.1814 0.1814 0.1814 0.1813 0.1812 0.1812 0.1811 0.1797 0.1796 0.1795 0.1795 0.1794 0.1793 0.1792 0.1791 0.1790 0.1789 0.1785 0.1784 0.1783 0.1782 0.1781 0.1780 0.1779 0.1778 0.1777 0.1777 0.1776 0.1775 0.1774 0.1773 0.1772 0.1772 0.1771 0.1770 0.1770 0.1769 0.1770 0.1769 0.1768 0.1767 0.1767 0.1766 0.1766 0.1765 0.1764 0.1764 0.1766 0.1765 0.1764 0.1763 0.1763 0.1762 0.1762 0.1761 0.1761 0.1760 0.1762 0.1762 0.1761 0.1761 0.1760 0.1760 0.1759 0.1759 0.1758 0.1758 0.1760 0.1759 0.1759 0.1758 0.1758 0.1758 0.1757 0.1757 0.1757 0.1756 result(:,:,3) = 0.0748 0.0749 0.0750 0.0750 0.0751 0.0752 0.0752 0.0753 0.0753 0.0754 0.0744 0.0745 0.0745 0.0746 0.0747 0.0747 0.0747 0.0748 0.0748 0.0749 0.0739 0.0740 0.0740 0.0740 0.0741 0.0741 0.0741 0.0741 0.0742 0.0742 0.0734 0.0735 0.0735 0.0735 0.0735 0.0735 0.0735 0.0735 0.0735 0.0735 0.0730 0.0730 0.0730 0.0730 0.0730 0.0730 0.0730 0.0730 0.0730 0.0730 0.0726 0.0726 0.0726 0.0726 0.0726 0.0726 0.0726 0.0726 0.0726 0.0726 0.0723 0.0723 0.0723 0.0723 0.0723 0.0723 0.0723 0.0723 0.0723 0.0723 0.0721 0.0721 0.0721 0.0721 0.0721 0.0720 0.0720 0.0720 0.0720 0.0720 0.0719 0.0719 0.0719 0.0719 0.0719 0.0718 0.0718 0.0718 0.0718 0.0718 0.0717 0.0717 0.0717 0.0717 0.0717 0.0717 0.0717 0.0717 0.0717 0.0717 result(:,:,4) = 0.0407 0.0407 0.0407 0.0408 0.0408 0.0408 0.0408 0.0408 0.0409 0.0409 0.0406 0.0406 0.0406 0.0407 0.0407 0.0407 0.0407 0.0407 0.0407 0.0407 0.0404 0.0405 0.0405 0.0405 0.0405 0.0405 0.0405 0.0405 0.0405 0.0405 0.0403 0.0403 0.0403 0.0403 0.0403 0.0403 0.0403 0.0403 0.0403 0.0403 0.0401 0.0401 0.0401 0.0401 0.0401 0.0401 0.0401 0.0401 0.0401 0.0401 0.0399 0.0399 0.0399 0.0399 0.0399 0.0399 0.0399 0.0399 0.0399 0.0399 0.0397 0.0397 0.0397 0.0397 0.0397 0.0397 0.0397 0.0397 0.0397 0.0397 0.0396 0.0396 0.0396 0.0396 0.0396 0.0396 0.0396 0.0396 0.0396 0.0396 0.0395 0.0395 0.0395 0.0395 0.0395 0.0395 0.0395 0.0395 0.0395 0.0395 0.0394 0.0394 0.0394 0.0394 0.0394 0.0394 0.0394 0.0394 0.0394 0.0394 result(:,:,5) = 0.0257 0.0257 0.0257 0.0257 0.0257 0.0257 0.0257 0.0257 0.0257 0.0257 0.0256 0.0256 0.0256 0.0256 0.0256 0.0257 0.0257 0.0257 0.0257 0.0257 0.0256 0.0256 0.0256 0.0256 0.0256 0.0256 0.0256 0.0256 0.0256 0.0256 0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 0.0254 0.0254 0.0254 0.0254 0.0254 0.0254 0.0254 0.0254 0.0254 0.0254 0.0253 0.0253 0.0253 0.0253 0.0253 0.0253 0.0253 0.0253 0.0253 0.0253 0.0252 0.0252 0.0252 0.0252 0.0252 0.0252 0.0252 0.0252 0.0252 0.0252 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0251 0.0250 0.0250 0.0250 0.0250 0.0250 0.0250 0.0250 0.0250 0.0250 0.0250 result(:,:,6) = 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0177 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0177 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0175 0.0175 0.0175 0.0175 0.0176 0.0176 0.0176 0.0176 0.0176 0.0176 0.0175 0.0175 0.0175 0.0175 0.0175 0.0175 0.0175 0.0175 0.0175 0.0175 0.0174 0.0174 0.0174 0.0174 0.0174 0.0175 0.0175 0.0175 0.0175 0.0175 0.0174 0.0174 0.0174 0.0174 0.0174 0.0174 0.0174 0.0174 0.0174 0.0174 0.0173 0.0173 0.0173 0.0173 0.0173 0.0174 0.0174 0.0174 0.0174 0.0174 0.0173 0.0173 0.0173 0.0173 0.0173 0.0173 0.0173 0.0173 0.0173 0.0173 result(:,:,7) = 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0128 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0129 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0128 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 0.0127 result(:,:,8) = 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0098 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0099 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0098 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 0.0097 result(:,:,9) = 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0077 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0078 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 0.0077 result(:,:,10) = 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0063 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062 0.0062
function Tr = fun(Ca,Ka,W,ma,m,k,Mu,W0,Zeta,Wa,Beta,ra,b,c,d,e,f,g)
ma=2;
m=40;
k=25;
Mu=ma/m;
W0=sqrt(k/m);
Zeta=Ca/(2*ma*W0);
Wa=sqrt(Ka/ma);
Beta= sqrt(Wa/W0);
r=W/W0;
Tr = sqrt(((2.*Zeta.*r).^2+(r.^2-Beta.^2).^2)./((2.*Zeta.*r).^2.*(r.^2-1+Mu.*r.^2).^2+(Mu.*r.^2.*Beta.^2-(r.^2-1).*(r.^2-Beta.^2)).^2));
%plot (W,Tr)
end

2 Comments
Show None Hide None

Tuncer
Tuncer on 22 Dec 2022
That seems to be solving my problem. Thanks a lot!
Rik
Rik on 22 Dec 2022
If this answer solved your question, please consider marking it as accepted answer. If not, feel free to comment with your remaining issues.

Sign in to comment.

More Answers (0)

Categories

Find more on Annotations in Help Center and File Exchange

Products

Release

R2021b

Tags

Asked:

Tuncer
Tuncer
on 20 Dec 2022

Commented:

Rik
Rik
on 22 Dec 2022

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!