2楼: Originally posted by cooooldog at 2014-10-25 18:50:25
1. 第一个是正常的现象; eps是实数双精度表示中的一个最小间隔,
http://www.mathworks.cn/help/matlab/ref/eps.html
数值计算,双精度时, 达到 sqrt(eps)的水平通常就已经很漂亮了;
如果想要更高精度, 你用mathem ...

3楼: Originally posted by cooooldog at 2014-10-25 19:00:32
第二个问题只知道你求一个一元三次方程,找到三个根; y(1:3)
每个根又对于x(1:3), z(1:3)

具体的物理意义是什么? 你想要的图又是什么?

4楼: Originally posted by cooooldog at 2014-10-25 19:08:24
如果要的是空间曲线;
把对应的值分别赋给:

y=;
x=;
z=;
然后
plot3(x,y,z)

8楼: Originally posted by 小楼寒 at 2014-10-25 19:59:33
第1题.同意楼上
第2题.修改后的代码，不知是不是你所要的结果：
clear
format long
syms k1 k01 VV  % z=dφ/dt
V=[];
z=[];
for j=1:1:40
    VVstep=0.0002;        %VV代表电极电势φ
    VVfinal=0.01 ...

9楼: Originally posted by 309065816 at 2014-10-25 22:34:16
嗯嗯，谢谢。我的基础比较差，所以编的程序都比较繁琐，按部就班的来。我明天运行一下试试看先。
1、变量名是之前老师有要求过，所以用这种，方便以后辨别。
2、我自己编程过程中也有虚部产生，可能是在运算过程 ...

close all
clear all
clc
%format long
%syms k1 k01 VV  % z=dφ/dt
maxiter = 40;
V=zeros(maxiter,1);
z=zeros(maxiter,3);

for jloop=1:1:maxiter
    VVstep=0.0002;        %VV代表电极电势φ
    VVfinal=0.0160;
    VV=0.008+VVstep*jloop;
    k1=150*exp(19.462*VV);
    k01=150*exp(-19.462*VV);
    equation_a=k01*200;
    equation_b=-(5*k1+1000+k01*200);
    equation_c=2*5*k1+1000+50*k01;
    equation_d=-5*k1;
    Coefficient=[equation_a equation_b equation_c equation_d];
    temp=roots(Coefficient);
   
    equation_y01=temp(1);
    equation_x01=(200*equation_y01^2-200*equation_y01+50)/5;
    equation_z01=(0.059-VV)/(1.1*10^(-7))-96485*(k1*(1-equation_y01)^2-k01*equation_x01*equation_y01)/(10*10^(-10));
   
    equation_y02=temp(2);
    equation_x02=(200*equation_y02^2-200*equation_y02+50)/5;
    equation_z02=(0.059-VV)/(1.1*10^(-7))-96485*(k1*(1-equation_y02)^2-k01*equation_x02*equation_y02)/(10*10^(-10));
   
    equation_y03=temp(3);
    equation_x03=(200*equation_y03^2-200*equation_y03+50)/5;
    equation_z03=(0.059-VV)/(1.1*10^(-7))-96485*(k1*(1-equation_y03)^2-k01*equation_x03*equation_y03)/(10*10^(-10));
   
   % hold on
    V(jloop)=VV;
    z(jloop,:)=[equation_z01,equation_z02,equation_z03];
end
disp('Display the four coefficients of cubic polynomial:')
for iloop=1:length(Coefficient)
    disp([num2str(iloop),': ',num2str(Coefficient(iloop),'%20.12f')])
end

% Find real values in z
index01 = find(abs(imag(z(:,1)))<sqrt(eps));
index02 = find(abs(imag(z(:,2)))<sqrt(eps));
index03 = find(abs(imag(z(:,3)))<sqrt(eps));
indices = intersect(intersect(index01,index02),index03);

plot(V(indices),z(indices,:))
legend('equation z01','equation z02','equation z03')