Matlab code - probably a simple mistake

2011 9 2
1 个回答
3 次查看（30 天）

0 个投票

I am trying to calculate kx, and ky. The conductivity in the x and y direction. The graph is supposed to look like a nice upwards curve, however mine fluxuates up and down(like a wave). I have checked my formulas, by inserting the values into the equation for L, the values of kx and ky appear to be correct. It seems as though the graphing part does not work tho.
I am going from L is 30 to 140, increments of 10. This is the only variable that changes. Here is my code:
clear
clc
% Constants
mfp = 31; % mean free path (nm)
k = 316.72; % conductivity (W / m*k)
% Equations
Lmax = 140;
dL=10;
L=30;
while L<=Lmax
kx(L)=k*(1-((2*mfp)/(3*3.14*L)));
ky(L)=k*(1-(mfp/(3*L)));
L= L+dL;
end
plot(kx,'-'),hold on, plot(ky,'--'), hold on, plot(k);
title(['Conductivity ']);
xlabel('Thickness (nm)','FontSize',14);
ylabel('Conductivity','FontSize',14);
leg1=legend('kx','ky');

请先登录，再回答此问题。

回答（1 个）

Paulo Silva
Paulo Silva 2011-9-2

0 个投票

You are just looking for plot(kx,ky) not all those plots
Also the calculations generate several zero values, that's why plot(kx,ky) shows many lines, you can plot just their values plot(kx,ky,'*') or ignore all the zeros and connect the points with one line:
clear
clc
% Constants
mfp = 31; % mean free path (nm)
k = 316.72; % conductivity (W / m*k)
% Equations
Lmax = 140;
dL=10;
L=30;
while L<=Lmax
kx(L)=k*(1-((2*mfp)/(3*3.14*L)));
ky(L)=k*(1-(mfp/(3*L)));
L= L+dL;
end
kx(~kx)=[];
ky(~ky)=[];
plot(kx,ky);
title(['Conductivity ']);
xlabel('Thickness (nm)','FontSize',14);
ylabel('Conductivity','FontSize',14);

5 个评论
显示 3更早的评论 隐藏 3更早的评论

Chris
Chris 2011-9-2
I am letting the while loop, go from 30 to 140. While computing values for kx, and ky . I am able to plot kx and ky, however, the plots are identical. However, when you manually solve the equations, they are not the same ( like it should be) . I am trying to plot kx and ky, dont worry about the extra plot for k.
Chris
Chris 2011-9-2
I am not by matlab right now, but what exactly is this doing now,
kx(~kx)=[];
ky(~ky)=[];
Paulo Silva
Paulo Silva 2011-9-2
that removes the zeros from the kx and ky vectors so you can get one line connecting all non zero points, I don't know if the results are correct that way, please check if they are correct, the line is there and it goes up like you said.
Chris
Chris 2011-9-2
when you say removes all the zeros, are you referring to all the zeros from 0 to about 30? This just makes the graph from 30 to 140 now? is that what you mean?
It should look like the graph of square root of x to a degree
Paulo Silva
Paulo Silva 2011-9-2
I don't know where the zeros are, I just noticed that they are messing up the result of the plot function and yes the result does look like plot(x,sqrt(x))

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

尚未输入任何标签。

提问:

Chris
Chris
2011-9-2

Community Treasure Hunt

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

Start Hunting!

Translated by