用多行代码编写一个表达式

1次观看(过去30天)

显示旧的注释

Sai maadesh kumar n r 2022年10月9日

0
链接

直接链接到这个问题

https://es.mathworks.com/matlabcentral/answers/1820465-writing-a-single-expression-in-multiple-lines-of-code

评论道: 明星黾 2022年10月9日

我必须定义一个矢量如下，作为我的代码的一部分，在一个LIVE脚本，必须导出为pdf:

Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2)); .．.

K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)) k3 * y (1) * (2) + Kprime3 * y(2);…

K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3)); .．.

K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2)); .．.

k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));

粗体部分是要聚焦的行。当我将它导出为pdf文件时，它的显示如下:

密码要流出边境了。最初，当第二行与其他三行水平对齐时，当该行的一部分没有打印到pdf中时。如果我给出三个点(…)，然后在下一行(给出回车键后)写入粗体行的一部分，如下所示:

Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2)); .．.

K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y(2))…

k3 * y (1) * (2) + Kprime3 * y(2);…

K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3)); .．.

K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2)); .．.

k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));

则输出如下:

如何将整个代码打印为pdf格式的实时脚本代码导出?

8的评论
显示隐藏旧的评论

Sai maadesh kumar n r 2022年10月9日

主程序

本程序计算了尼龙6的反应物转化率和平均分子量混合均匀的等温批处理反应堆

                                此代码使用RK-4方法解决ode
                               
                                tEnd = 12;% HRS时间
                               
                                H = 0.01;%步长;
                               
                                N = tEnd/h;%。iter。
                               
                                %初始条件
                               
                                Y0 =[8.8 0 0 0 0.177]';% y向量的初始条件
                               
                                t = 0 (1, N);%的时间向量
                               
                                N y = 0(5日);%y-array: row- time, column - y1, y2…
                               
                                y (: 1) = y0;初始化y @t=0的第一列
                               
                                导数的%函数
                               
                                为i = 1: N
                               
                                按步长值改变时间
                               
                                t (i + 1) = t(我)+ h;%行向量
                               
                                k1 = f (t (i), y(:,我));
                               
                                k2 = f (t(我)+ (h / 2), y(:,我)+ (h / 2) * k1);
                               
                                k3 = f (t(我)+ (h / 2), y(:,我)+ (h / 2) * k2);
                               
                                k4 = f (t (i) + h, y(:,我)+ h * k3);
                               
                                使用RK-4中斜率的加权平均值更新y
                               
                                Y (:，i+1) = Y (:，i)+(h/6)*(k1+2*k2+2*k3+k4);
                               
                                结束
                               
                                %subplot用于绘制两个图形
                               
                                图(1)
                               
                                情节(t (1 - ((1:) / y (1,1))),“r”）;
                               
                                包含(的时间(人力资源)）;
                               
                                ylabel (“caprolactum转换”）;
                               
                                标题(“Caprolactum转换”）;
                               
                                图(2)
                               
                                情节(t(113 *(4:)。/ y (5,:)));
                               
                                包含(的时间(人力资源)）;
                               
                                ylabel (“Avg.MW”）;
                               
                                标题(“Avg兆瓦。”）;

函数来计算导数

                                函数Ddt = f (t, y)
                               
                                %的常量
                               
                                T = 538;%开尔文
                               
                                %参数，用于计算速率常数
                               
                                A0=[5.99e +05, 1.89e +10, 2.29e +09];
                               
                                E0=[1.99e +04, 2.33e +04, 2.28e +04];
                               
                                Ac=[4.31E+07, 1.21E+04, 1.64E+10];
                               
                                Ec=[1.88E+04, 2.07E+04, 2.01E+04];
                               
                                H=[1.92 2e +03， -5.95e +03， -4.04e +03];
                               
                                S=[-7.89, 0.944， -6.95];
                               
                                y3 = y (3);在计算速率常数时，只使用这个y
                               
                                %速率常数(大写字母与RK法的k斜率区别)
                               
                                K1 = kCalc (y3、E0 A0, Ac, Ec, T, 1);
                               
                                K2 = kCalc (y3、E0 A0, Ac, Ec, T, 2);
                               
                                K3 = kCalc (y3、E0 A0, Ac, Ec, T, 3);
                               
                                Kprime1 = kprimeCalc (K1、S、H T 1);
                               
                                Kprime2 = kprimeCalc (K1、S、H T 2);
                               
                                Kprime3 = kprimeCalc (K1, S、H T, 3);
                               
                                %包含每个dy/dt的向量
                               
                                Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2));.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)) k3 * y (1) * (2) + Kprime3 * y (2);.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3));.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2));.．.
                               
                                k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));
                               
                                %函数来确定k
                               
                                函数k = kCalc (y3 E0、A0 Ac, Ec, T,我)
                               
                                k = A0 (i) * exp (e0(我)/ (1.987 * T)) + Ac (i) * y3 * exp(中(我)/ (1.987 * T));
                               
                                结束
                               
                                确定k'的%函数
                               
                                函数kprime = kprimeCalc (k S H T i)
                               
                                kprime = k / (exp ((S (i) / 1.987) - (H (i) / (1.987 * T))));
                               
                                结束
                               
                                结束

明星黾 2022年10月9日

对于我来说，代码运行没有错误，关于代码本身和格式-

此代码使用RK-4方法解决ode

tEnd = 12;% HRS时间

H = 0.01;%步长;

N = tEnd/h;%。iter。

%初始条件

Y0 =[8.8 0 0 0 0.177]';% y向量的初始条件

t = 0 (1, N);%的时间向量

N y = 0(5日);%y-array: row- time, column - y1, y2…

y (: 1) = y0;初始化y @t=0的第一列

导数的%函数

为i = 1: N

按步长值改变时间

t (i + 1) = t(我)+ h;%行向量

k1 = f (t (i), y(:,我));

k2 = f (t(我)+ (h / 2), y(:,我)+ (h / 2) * k1);

k3 = f (t(我)+ (h / 2), y(:,我)+ (h / 2) * k2);

k4 = f (t (i) + h, y(:,我)+ h * k3);

使用RK-4中斜率的加权平均值更新y

Y (:，i+1) = Y (:，i)+(h/6)*(k1+2*k2+2*k3+k4);

结束

%subplot用于绘制两个图形

图(1)

情节(t (1 - ((1:) / y (1,1))),“r”）;

包含(的时间(人力资源)）;

ylabel (“caprolactum转换”）;

标题(“Caprolactum转换”）;

图(2)

情节(t(113 *(4:)。/ y (5,:)));

包含(的时间(人力资源)）;

ylabel (“Avg.MW”）;

标题(“Avg兆瓦。”）;

函数Ddt = f (t, y)

%的常量

T = 538;%开尔文

%参数，用于计算速率常数

A0=[5.99e +05, 1.89e +10, 2.29e +09];

E0=[1.99e +04, 2.33e +04, 2.28e +04];

Ac=[4.31E+07, 1.21E+04, 1.64E+10];

Ec=[1.88E+04, 2.07E+04, 2.01E+04];

H=[1.92 2e +03， -5.95e +03， -4.04e +03];

S=[-7.89, 0.944， -6.95];

y3 = y (3);在计算速率常数时，只使用这个y

%速率常数(大写字母与RK法的k斜率区别)

K1 = kCalc (y3、E0 A0, Ac, Ec, T, 1);

K2 = kCalc (y3、E0 A0, Ac, Ec, T, 2);

K3 = kCalc (y3、E0 A0, Ac, Ec, T, 3);

Kprime1 = kprimeCalc (K1、S、H T 1);

Kprime2 = kprimeCalc (K1、S、H T 2);

Kprime3 = kprimeCalc (K1, S、H T, 3);

%包含每个dy/dt的向量

Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2));

K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)) k3 * y (1) * (2) + Kprime3 * y (2);

K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3));

K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2));

k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));

%函数来确定k

函数k = kCalc (y3 E0、A0 Ac, Ec, T,我)

k = A0 (i) * exp (e0(我)/ (1.987 * T)) + Ac (i) * y3 * exp(中(我)/ (1.987 * T));

结束

确定k'的%函数

函数kprime = kprimeCalc (k S H T i)

kprime = k / (exp ((S (i) / 1.987) - (H (i) / (1.987 * T))));

结束

我消除了延拓省略号( .．. )，因为没有理由拥有它们。

有什么问题吗?

．

Torsten 2022年10月9日

编辑:Torsten 2022年10月9日

@sai maadesh kumar n r

使用

                                Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2));
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)).．.
                               
                                K3 * y (1) * (2) + Kprime3 * y (2);
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3));
                               
                                K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2));
                               
                                k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));

而不是

                                Ddt = (k1 * y (1) * (5) + Kprime1 * y y (1) (2) k3 * * y (3) + Kprime3 * (y (3) - y (2));.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)).．.
                               
                                k3 * y (1) * (2) + Kprime3 * y (2);.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y (2) k2 * y (3) ^ 2 + Kprime2 * y (5) * (y (4) - y (3));.．.
                               
                                K1 * y (1) * (5) -Kprime1 * y y (1) (2) + K3 * * y (3) -Kprime3 * (y (3) - y (2));.．.
                               
                                k1 * y (1) * (5) + Kprime1 * y (2) + K2 * y (3) ^ 2-Kprime2 * y (5) * (y (4) - y (3)));

最有可能的是线条

                                K1 * y (1) * (5) -Kprime1 * y (2) 2 * K2 * y y (3) + 2 (2) * * Kprime2 * y (5) * (y (3) - y (2)).．.
                               
                                k3 * y (1) * (2) + Kprime3 * y (2);.．.

将从MATLAB中解释为两列-因此Ddt将是一个大小为(5x2)而不是(5x1)的矩阵。

明星黾 2022年10月9日

我的荣幸!

登录评论。

登录回答这个问题。

答案(0)

登录回答这个问题。

类别

MATLAB 图形图形对象图形对象编程

了解更多图形对象编程在帮助中心而且文件交换

2022世界杯八强谁会赢？

MATLAB

释放

R2022a

世界杯预选赛小组名单社区寻宝

在MATLAB Central中找到宝藏，并发现社区如何帮助您!世界杯预选赛小组名单

开始狩猎!

翻译的

用多行代码编写一个表达式

8的评论
显示隐藏旧的评论

答案(0)

另请参阅

类别

标签

2022世界杯八强谁会赢？

释放

世界杯预选赛小组名单社区寻宝

用多行代码编写一个表达式

8的评论 显示隐藏旧的评论

答案(0)

另请参阅

类别

标签

2022世界杯八强谁会赢？

释放

世界杯预选赛小组名单社区寻宝

8的评论
显示隐藏旧的评论