同心圆的问题

没落贵族 · 发表于 2010-12-9 13:38

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

用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的时候是圆心，请大家帮帮忙，谢谢光强度的变化.bmp

qibbxxt · 发表于 2010-12-9 13:57

做个循环画，不就可以了吗
不同的圆传入不同的半径和颜色不就可以了吗？
另外，不明白A这个行列式是什么意思，行列式就是一个数，不知道该怎么用？

没落贵族 · 发表于 2010-12-9 14:15

A行列里的每一个数都代表颜色的变化。再问一下，怎么把不同的颜色由明到暗带入各个圆里？？

leeking30 · 发表于 2010-12-9 20:18

本帖最后由 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;
[X,Y] = pol2cart(THETA,RHO);
Z=ones(1,1000)*A(i);
plot3(X,Y,Z,'linewidth',1.5,'color',[0 1 0].* E(i) );hold on
end
view([0 0 1])

代码贴出来不用下载费体能了

没落贵族 · 发表于 2010-12-10 15:46

谢谢leeking30了，非常感谢

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 同心圆的问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分