同心圆的问题
用matlab做一个100个圆的同心圆,各个圆的颜色由给定的行列式的值来赋予,怎么做?现在有行列式A,根据A的各个值从外圆开始给各个圆付上颜色,用来观察同心圆从内到外的光变化。请问有什么好方法吗?
A=[ 1.33E-41 2.14E-41 1.83E-41 8.16E-42 7.66E-42 2.26E-41 3.80E-41 3.78E-41 2.45E-41 1.14E-41 3.73E-42 1.95E-43 3.72E-42 1.54E-41 2.35E-41 1.72E-41 9.64E-42 2.60E-41 6.45E-41 9.07E-41 8.02E-41 4.80E-41 2.22E-41 9.34E-42 2.88E-42 6.31E-42 2.17E-41 3.08E-41 2.24E-41 2.84E-41 9.03E-41 1.86E-40 2.35E-40 1.99E-40 1.27E-40 8.37E-41 7.17E-41 6.28E-41 6.24E-41 8.37E-41 9.63E-41 8.39E-41 1.31E-40 3.45E-40 6.67E-40 8.79E-40 8.65E-40 7.65E-40 7.76E-40 9.01E-40 1.00E-39 1.05E-39 1.13E-39 1.18E-39 1.11E-39 1.07E-39 1.54E-39 2.72E-39 4.17E-39 5.21E-39 5.78E-39 6.51E-39 7.88E-39 9.63E-39 1.13E-38 1.27E-38 1.40E-38 1.48E-38 1.45E-38 1.28E-38 1.07E-38 9.57E-39 9.60E-39 9.91E-39 9.51E-39 8.27E-39 6.78E-39 5.60E-39 4.84E-39 4.23E-39 3.62E-39 3.19E-39 3.47E-39 4.74E-39 6.47E-39 7.45E-39 6.80E-39 4.86E-39 2.80E-39 1.55E-39 1.10E-39 1.03E-39 1.11E-39 1.60E-39 3.08E-39 6.45E-39 1.24E-38 2.03E-38 2.73E-38 3.02E-38];
如图,横轴代表圆的个数,纵轴代表颜色的变化。x=100的时候是圆心,请大家帮帮忙,谢谢
做个循环画,不就可以了吗
不同的圆传入不同的半径和颜色不就可以了吗?
另外,不明白A这个行列式是什么意思,行列式就是一个数,不知道该怎么用? A行列里的每一个数都代表颜色的变化。再问一下,怎么把不同的颜色由明到暗带入各个圆里??
本帖最后由 leeking30 于 2010-12-9 20:19 编辑
抛砖引玉 最简单的办法
clear;clc
A=[ 1.33E-41 2.14E-41 1.83E-41 8.16E-42 7.66E-42 2.26E-41 3.80E-41 3.78E-41 2.45E-41 1.14E-41 3.73E-42 1.95E-43 3.72E-42 1.54E-41 2.35E-41 1.72E-41 9.64E-42 2.60E-41 6.45E-41 9.07E-41 8.02E-41 4.80E-41 2.22E-41 9.34E-42 2.88E-42 6.31E-42 2.17E-41 3.08E-41 2.24E-41 2.84E-41 9.03E-41 1.86E-40 2.35E-40 1.99E-40 1.27E-40 8.37E-41 7.17E-41 6.28E-41 6.24E-41 8.37E-41 9.63E-41 8.39E-41 1.31E-40 3.45E-40 6.67E-40 8.79E-40 8.65E-40 7.65E-40 7.76E-40 9.01E-40 1.00E-39 1.05E-39 1.13E-39 1.18E-39 1.11E-39 1.07E-39 1.54E-39 2.72E-39 4.17E-39 5.21E-39 5.78E-39 6.51E-39 7.88E-39 9.63E-39 1.13E-38 1.27E-38 1.40E-38 1.48E-38 1.45E-38 1.28E-38 1.07E-38 9.57E-39 9.60E-39 9.91E-39 9.51E-39 8.27E-39 6.78E-39 5.60E-39 4.84E-39 4.23E-39 3.62E-39 3.19E-39 3.47E-39 4.74E-39 6.47E-39 7.45E-39 6.80E-39 4.86E-39 2.80E-39 1.55E-39 1.10E-39 1.03E-39 1.11E-39 1.60E-39 3.08E-39 6.45E-39 1.24E-38 2.03E-38 2.73E-38 3.02E-38];
A=A*1E38;
E=A/max(A);
for i=1:1:100
radius=i;
THETA=linspace(0,2*pi,1000);
RHO=ones(1,1000)*radius;
= pol2cart(THETA,RHO);
Z=ones(1,1000)*A(i);
plot3(X,Y,Z,'linewidth',1.5,'color',.* E(i) );hold on
end
view()
代码贴出来 不用下载费体能了
谢谢leeking30了,非常感谢
页:
[1]