简单的MATLAB程序，我不会matlab

Question

程序如下，我不太会，或者可以说不会。我自己依葫芦画瓢编的运行不起。希望可以改改：
%直流电动机转速特性曲线
%以n=f（Ia）关系画出曲线
%--------------------------------
%当改变
%下面输入电机的基本数据：
Ф=1;Ce=1;Ra=0.0161;
%规定变量变化范围：
Ia=0:100:2000;U1=750;U2=800;U3=700;
%带入公式，计算出直流电动机转速特性曲线：
n1=(U1/(Ce*Ф))-(Ra.*Ia)./(Ce*Ф);
n2=(U2/(Ce*Ф))-(Ra.*Ia)./(Ce*Ф);
n3=(U3/(Ce*Ф))-(Ra.*Ia)./(Ce*Ф);
plot(Ia,n1,n2，n3，'k')
xlabel('Ia')
ylabel('n')
grid on

Answer 1

虽然我不懂你的程序，但是你里边有几个问题。首先，我复制粘贴到MATLAB之后字母Ф都是现实红色的（红色是有问题咯）。我随便用个字母Q代替了下。还有就是输出图像plot(Ia,n1,n2，n3，'k')这一行中n2，n3，这两个，都是在中文输入法下打出来的。当然也不行，这些符号仅能在英文状态下输入。然后运行有输出。图像在附件中。希望对你有帮助。

>> %直流电动机转速特性曲线

%以n=f（Ia）关系画出曲线

%--------------------------------

%当改变

%下面输入电机的基本数据：

q=1;Ce=1;Ra=0.0161;

%规定变量变化范围：

Ia=0:100:2000;U1=750;U2=800;U3=700;

%带入公式，计算出直流电动机转速特性曲线：

n1=(U1/(Ce*q))-(Ra.*Ia)./(Ce*q);

n2=(U2/(Ce*q))-(Ra.*Ia)./(Ce*q);

n3=(U3/(Ce*q))-(Ra.*Ia)./(Ce*q);

plot(Ia,n1,n2,n3,'k')

xlabel('Ia')

ylabel('n') 

grid on

Answer 2

把 Ф 改为 theta 没有问题了：
theta=1;Ce=1;Ra=0.0161;
%规定变量变化范围：
Ia=0:100:2000;U1=750;U2=800;U3=700;
%带入公式，计算出直流电动机转速特性曲线：
n1=(U1/(Ce*theta))-(Ra.*Ia)./(Ce*theta);
n2=(U2/(Ce*theta))-(Ra.*Ia)./(Ce*theta);
n3=(U3/(Ce*theta))-(Ra.*Ia)./(Ce*theta);
plot(Ia,n1,n2,n3,'k')
xlabel('Ia')
ylabel('n')
grid on